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PREFACE 


This  report  details  the  technical  effort  put  forth  by  CSD  and  its 
subcontractors  during  the  four  phases  of  this  program  for  propellant  nonlinear 
constitutive  theory  extension.  This  includes  the  Phase  I  preliminary  study  in 
which  the  Quinlan  theory  was  critiqued,  alternate  approaches  were  studied  and 
detailed  research  planning  accomplished.  Also  included  are  the  detailed  experi¬ 
mental  evaluations  of  propellant  during  Phases  II  and  III,  the  uniaxial/ 
isothermal  investigation  and  the  two-dimensional  variable  temperature 
investigation.  Detailed  subcontractor  theoretical  development  and  predictions 
are  presented  along  with  the  three-dimensional  investigations. 

This  final  technical  report  on  Contract  No.  F0^6l 1-80-C-0052  consists  of  a 
summary  of  the  phase  II  and  III  laboratory  propellant  evaluation  for  UTP-3001 
and  UTP-19,360B  propellants.  Phase  IV  was  a  validation  of  the  accuracy  of  the 
constitutive  theory  in  a  three-dimensional  state  of  stress  and  strain.  The 
detailed  experimental  results  were  distributed  in  Data  Packages  A  through  G  to 
all  project  personnel.  Additional  details  are  given  in  section  1.0. 

In  performing  the  work  required  for  the  program,  Chemical  Systems  Division 
(CSD)  employed  under  subcontract  the  services  of  five  scientists  of  national  rep¬ 
utation:  Drs.  M.  Quinlan,  M.  E.  Curtin,  W.  L.  Hufferd,  R.  Wool,  and  R.  A. 
Schapery.  The  program  effort  combined  the  comprehensive  experience  and  special¬ 
ized  test  capabilities  of  CSD  in  solid  propellant  mechanical  properties  with  the 
theoretical  expertise  of  these  scientists.  In  addition,  Dr.  J.  E.  Fitzgerald 
was  retained  a3  a  consultant  to  participate  in  the  periodic  technical  reviews  of 
the  program  status. 


1.0  PROGRAM  OBJECTIVES  AND  OVERVIEW 


The  objective  of  the  Propellant  Nonlinear  Constitutive  Theory  (NLCT) 

Program  was  to  develop  and  demonstrate  a  nonlinear  thermomechanical  constitutive 
law  for  solid  rocket  propellants.  The  goal  was  to  achieve  errors  of  less  than 
10%  in  stress  predictions  for  undamaged  propellant  and  less  than  15 1  for 
previously  damaged  propellant.  The  program  was  conducted  in  four  phases: 

•  Phase  I  -  Preliminary  Study 

•  Phase  II  -  Uniaxial  Isothermal  Investigation 

•  Phase  III  -  Two-Dimensional  and  Variable  Temperature  Investigation 

e  Phase  IV  -  Three-Dimensional  Investigation. 

Each  phase  had  a  number  of  tasks  and  subtasks.  Overall  program  logic  is  shown 
in  Figure  1 . 

A  number  of  subcontractors  were  employed  to  evaluate  the  several  NLCT 
approaches  considered,  and  Dr.  J.  E.  Fitzgerald,  Director  of  the  School  of  Civil 
Engineering,  Georgia  Institute  of  Technology,  was  retained  as  a  technical 
advisor  and  consultant  to  the  NLCT  program. 

The  objectives  of  the  four  phases  are  summarized  in  the  following 
paragraphs: 

1.1  PHASE  I  -  PRELIMINARY  STUDY 

The  objective  of  phase  I  was  to  make  detailed  research  plans  for  evaluating 
and  modifying  candidate  nonlinear  constitutive  laws.  Quinlan's  nonlinear 
constitutive  theory11  was  critically  reviewed  and  modifications  to  the 
differential  equation  for  the  bonding  parameter  were  proposed.  In  addition, 
five  other  candidate  constitutive  approaches  were  evaluated  on  a  preliminary 
basis  by  subcontractors  and  CSD  personnel: 

%  Farris  Theory  -  W.  Briggs 

•  Modified  Swanson  Theory  -  D.  Gutierrez-Lemini 

•  Russian  Theories  -  W.  L.  Hufferd 
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•  Time-Dependent  Bonding  -  R.  A.  Schapery 

•  Softening  Approach  -  M.  E.  Gurtin. 

Dr.  S.  M.  Breitling  and  Dr.  R.  P.  Wool  also  provided  information  of  the  probability 
of  successfully  incorporating  a  physically  based  molecular  model  into  any  of  the 
NLCT's  considered  for  evaluation. 

1.2  PHASE  II  -  UNIAXIAL/ISOTHERMAL  INVESTIGATION 

The  objectives  of  phase  II  were  to  continue  theoretical  development  w or*  or. 
the  candidate  nonlinear  constitutive  theories  and  to  evaluate  their  predictive 
capabilities  for  the  uniaxial  isothermal  response  to  two  propellants:  an  HTPB 
(UTP-19,360B)  and  a  PBAN  (UTP-3001).  An  extensive  test  program  was  conducted  on 
these  two  propellants.  Following  the  comparative  evaluation  of  the  constitutive 
theories,  successful  candidates  were  recommended  for  continuation  into 
phase  III. 

1.3  PHASE  III  -  TWO-DIMENSIONAL  A> ^  "ARIABLE  TEMPERATURE  INVESTIGATION 

The  objective  of  phase  III  was  tc  extend  the  theoretical  development  of  the 
candidate  constitutive  theories  to  two-dimensional  and  variable  temperature 
histories.  Stress  time  predictions  were  compared  with  laboratory  tests  of 
the  two  propellants  and  a  decision  was  reached  as  to  which  theories  would 
be  inherent  within  phase  IV. 

1 .4  PHASE  IV  -  THREE-DIMENSIONAL  INVESTIGATION 

The  objective  of  the  final  phase  was  to  conduct  a  critical  validation 
experiment  and  assess  the  accuracy  of  the  constitutive  theory  developed  under 
isothermal  and  transient  temperature  conditions  for  a  three-dimensional  state  of 
stress  and  strain. 


Figure  1 .  Program  Logic 


2.0  SUMMARY,  CONCLUSIONS,  AND  RECOMMENDATIONS 


2.1  SUMMARY  AND  CONCLUSIONS 

The  objective  of  the  Propellant  Nonlinear  Constitutive  Theory  Extension 
Program  was  to  develop  and  demonstrate  an  accurate,  usable,  three-dimensional, 
thermomechanical  constitutive  theory  for  solid  rocket  propellants.  The  target 
criteria  for  selection  of  a  candidate  constitutive  relation  were  two-fold: 

•  Less  than  10$  error  in  stress  predictions  for  undamaged  propellant 

•  Less  than  15$  error  in  stress  predictions  for  previously  damaged 
propellant . 

The  program  was  conducted  in  four  phases. 

The  basis  for  the  program  was  the  limited  success  achieved  by  Quinlan  in 
AFRPL-TR-78-371*  in  modeling  some  of  the  trends  typically  observed  of  solid 
propellants  using  a  constitutive  theory  incorporating  damage.  Accordingly, 
a  significant  effort  was  spent  by  Dr.  M.  Quinlan  at  Cork  University,  Ireland, 
and  by  CSD  in  an  attempt  to  extend  and  improve  the  predictive  capabilities 
of  Quinlan's  fundamental  model.  Though  some  advances  were  made  at  the  theo¬ 
retical  level  to  increase  our  knowledge  of  the  model,  the  model  did  not  achieve 
the  desired  success  and  was  discontinued  before  the  transition  to  combined 
straining-cooling  and  three-dimensional  stress  states. 

Similarly,  Gurtin's  maximum  strain/stress  softening  approach  met  with  only 
marginally  greater  success.  The  model  evolution  was  capable  of  describing  some 
trends  of  solid  propellant  behavior,  but  numerical  difficulties  were  encountered 
in  the  curve  fitting  procedures.  This  method  was  also  discontinued  after 
working  with  complex  uniaxial  test  histories. 

The  remaining  approaches  (namely  those  of  Schapery,  Hufferd,  and  CSD’s 
modification  of  Swanson's  constitutive  theory)  all  achieved  greater  success. 

Due  to  the  extensive  effort  (in  time  and  money)  expended  during  phases  I  and  II 
with  five  analytical  models,  only  CSD's  modification  of  Swanson's  theory  was 
carried  through  to  actual  motor  stress  predictions. 
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Schapery's  approach  is  based  on  a  micro-cracking  model  for  damage.  The 
theory  predicted  stresses  accurately  for  uniaxial  and  biaxial  isothermal  tests. 
The  approach  was  not  continued  to  transient  thermal  nor  motor  stress  predictions 
due  to  limitations  in  time  and  funding. 

Hufferd's  modification  of  Il'yushin's  approach  to  thermovisco-plasticity 
provided  good  predictive  capabilities  for  uniaxial  isothermal  test  histories  and 
was  comparable  to  the  modified  Swanson  theory  for  isothermal  biaxial  and  tran¬ 
sient  cooling/straining  tests.  The  approach  was  not  continued  to  motor  stress 
predictions  due  to  a  lack  of  funds. 

The  Swanson  approach  modified  by  CSD  achieved  the  greatest  success  of  all 
the  nonlinear  analytical  methods.  This  approach  established  a  major  milestone 
in  the  solid  propellant  industry.  Using  the  characterization  data  fiom  three 
simple  laboratory  test  modes,  a  wide  variety  of  very  complex  uniaxial  test 
histories  were  predicted  to  the  desired  accuracy  levels.  These  successful  pre¬ 
dictions  also  included  stress  and  t  '«tt*  response  ranges  outside  of  the 
characterization  ranges.  This  was  the  first  known  accurate  solid  propellant 
stress  prediction  which  extended  beyond  the  nonlinear  characterized  data  base. 
This  analytical  model  also  accurately  predicted  biaxial  and  shear  test  histories 
which  were  not  included  in  the  characterization. 

Because  of  budget  and  schedule  limitations,  this  method  was  not  completely 
adapted  for  thermomechanical  coupling  and  complete  three-dimensional  analysis. 

2.2  RECOMMENDATIONS 

While  predictive  success  was  achieved  with  the  modified  Swanson  theory 
predicting  complex  uniaxial,  biaxial  and  shear  test  histories,  there  was  only 
preliminary  experience  with  full  three-dimensional  stress  axiality  with  combined 
thermal  mechanical  loads.  It  is  recommended  that  additional  work  with  the 
modified  Swanson  theory  be  conducted  with  emphasis  on  combined  thermal  and 
meohanical  loads  and  realistic  three-dimensional  loading  axialities. 

All  four  nonlinear  analysis  methods  could  be  further  developed  and  could 
provide  alternative  methods  for  analysis  of  solid  propellant  structural 
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reapor.se .  The  modified  Swanson  method  is  potentially  compatible  with  inclusion 
in  a  finite  element  analysis,  but  further  work  mentioned  above  is  required 
before  this  transition  can  be  successful. 


3.0  EXPERIMENTAL  PROGRAM 


Laboratory  testing  was  divided  into  categories  of  uniaxial/isothermal, 
two-dimensional ,  and  variable  temperature  and  three-dimensional  investigations. 

The  two  propellants  selected  for  the  program  were  (1)  a  PBAN  used  in  the 
first  stage  of  the  Titan  missile  system  (UTP-3001-750/7768)  and  (2)  a  HTPB 
propellant  developed  for  the  IUS  motor  (UTP-19 ,360B-M00/1777) .  The  first 
numbers  are  the  propellant  designation;  the  next,  the  mixer  size;  and  the 
last,  a  batch  number. 

In  each  of  the  laboratory  test  groups  a  specific  test  of  each  type  has 
been  selected  to  show  the  test  details. 

3.1  UNIAXIAL/ISOTHERMAL  INVESTIGATION 

Testing  uniaxial  specimens  of  UTP-3001  and  UTP-19, 360B  propellants  in 
phase  II  of  the  contract  was  done  for  the  non-damaged  material  (as  indicated  in 
Figure  2)  and  for  damaged  material  (per  Figures  3  and  4).  Most  of  the  tests 
were  run  with  1/2-  x  1/2-  x  6-in.  bars  with  redwood  end  tabs.  The  exceptions  were 
the  stress  endurance  (test  2)  and  constant  rate  (comparison  to  test  4)  data 
which  were  obtained  with  JANNAF  Class  B  specimens.^  Details  of  most  of  the 
individual  test  types  are  discussed  in  subsections  below. 

Details  of  the  other  tests  (which  were  considered  in  earlier  work  but  not 
utilized  in  this  report)  may  be  obtained  from  AFRPL-TR-83-031l32.  All  of 
the  tests  are  shown  in  Figures  2  to  *1  with  sketches  of  the  strain-time  histories. 

The  attachment  linkages  on  both  testers  were  such  that  the  specimens  could 
not  be  put  into  compression  when  the  crosshead  was  returned  to  an  equivalent 
zero  strain  position.  The  strain  measurement  was  done  with  a  linear  potentiome¬ 
ter  attached  to  the  crosshead;  consequently,  the  data  had  to  be  modified  to 
reflect  the  propellant  strain  relaxation  behavior  after  the  stress  had  returned 
to  zero  (free  hanging  specimen). 


Figure  2.  Uniaxial/Isothermal  Non-damaged  Tests 

24^406  R 1 


The  uniaxial  bars  were  machined  from  redwood  boxes  of  propellant.  Tne 
redwood  was  sealed  then  lined  in  the  same  manner  as  a  rocket  motor.  After 
a  partial  cure  of  the  liner,  propellant  was  cast  into  the  box  and  the  system 
cured  to  provide  a  good  bond  to  the  redwood  end  tabs.  The  redwood  box  assembly 
and  finished  specimen  are  shown  in  Figure  5.  After  the  specimen  is  mill 
finished,  a  1/8-in.  hole  is  drilled  on  the  centerline  of  each  end  tab  for 
attachment  to  the  testing  machine. 


All  tests  were  conducted  on  a  Chemical  Systems  Division  (CSD)  manufactured 
six-channel  tester  and  a  modified  Instron.  A  Hewlett-Packard  computer  war,  used 
to  collect  digitized  data  from  the  tests  (see  Appendix  A  details  on  the 
automated  data  reduction  system). 


Figure  3.  Uniaxial/Isothermal  Damaged  Tests 


Nominal  tests  were  run  with  three  samples  per  set 
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Figure  5.  Uniaxial  Bar  Specimen 
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Strain  relaxation  was  measured  on  samples  in  some  of  the  tests  during  final 
unload  cycle.  Cathetometer  measurements  were  made  periodically  and  strain  ver¬ 
sus  time  data  were  plotted.  These  data  were  used  to  estimate  the  relaxation 
behavior  on  cyclic  tests  where  there  was  insufficient  time  for  measurements. 


A  data  modification  was  made  to  estimate  the  peak  or  minimum  stress  and 
strain  points  which  were  not  recorded  by  the  digitized  data  acquisition  system. 
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The  sampling  rate  limited  the  crosshead  rate  that  could  be  used  and  still  obtain 
enough  points  to  adequately  define  a  ramp.  The  available  computer  memory  also 
influenced  the  sampling  rate  in  some  of  the  longer  tests. 

3.1.1  Constant  Rate  Test  No.  1 

Uniaxial  constant  rate  to  failure  tests  wore  conducted  on  6-in.  bars  of 
UTP-3001  and  UTP-19,360B.  The  70  F  tests  were  at  crosshead  rates  of 
10,  1,  0.1,  0.01  and  0.001  in. /min,  while  the  M0  and  120  F  tests  were  at 
10,  1,  and  0.1  in. /min. 

Typical  stress-strain  curves  for  the  3  rates  at  120  F  are  shown  in 
Figure  6  for  UTP-3001  .  To  ensure  an  equivalent  time  data  base,  three  sampler,  of 
each  propellant  were  tested  at  the  same  time  in  the  six  channel  testers.  The 
exception  was  the  10  in. /min  tests  which  were  tested  one  propellant  at  a  time 
due  to  data  sampling  rate  limitations. 

The  6-in.  bar  specimens  always  '  led  at  strain  levels  below  the  level  that 
would  be  obtained  from  JANNAF  specimens  and  always  at  the  propellant  to  wooden 
tab  interface.  Since  the  objective  of  the  tests  was  to  obtain  response  data  and 
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Figure  6.  Constant-Rate  Tests  of  UTP-3001-750/7768  at  124  F  with 

6-in.  Bar  Specimens 
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Figure  8.  Test  No.  3  -  Low-High  Constant-Rate  Tests  of 
UTP-19 , 360B-400/1777  at  70  F 


propellants.  The  samples  were  loaded  to  3^  strain  at  a  crosshead  rate  of 
1  in. /min  for  temperatures  of  20,  43,  73,  and  122  F.  The  load  was  monitored 
with  the  time  while  3train  was  determined  from  cathetometer  measurements 

on  the  samples. 


The  master  stress  relaxation  modulus  data  for  UTP-3001  are  presented  in 
Figure  9  as  typical.  Actual  data  are  available  from  the  data  bank. 

3.1.4  Constant  Rate  Secant  Modulus 

Constant  rate  modulus  tests  were  run  on  the  6-in.  bar  samples  described 
above  and  on  JANNAF  class  B  specimens.  The  6-in.  bar  samples  used  the 


Figure  9.  Master  Modulus  Data  for  UTP-3001 -750-/7768  with 

Experimental  Shift 
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wood-to-wood  distance  as  the  effective  gage  length  while  the  JANNAF's  were 
analyzed  using  2. 70- in.  effective  gage  length  through  the  test  even  though 
it  is  varying.  Strain  was  determined  by  the  crosshead  travel. 

Constant  rate  modulus  data  for  UTP-3001,  with  the  curve  drawn  through 
the  small  strain  portion  of  the  results,  are  compared  to  the  relaxation  modulus 
in  Figure  10.  Similar  data  for  UTP-19,360B  are  available  in  the  data  bank. 

3.1.5  Cyclic  Loading  Test  No.  5 

The  multiple  loading  tests  on  6-in.  bars  of  UTF-3001  and  UTP-19.360B  were 
run  in  five  cycles  with  increasing  strain  levels  for  each  cycle  and  a  rest 
period  between  cycles.  All  tests  were  at  70  F  and  at  crosshead  rates  of 
5,  1,  and  0.1  in. /min.  An  attempt  was  made  with  UTP-3001  to  run  at  10  in. /min 
(i.e.,  planned  rate  instead  of  five),  but  the  data  sampling  rate  did  not  provide 
sufficient  data  points  to  clearly  define  the  load-time  curve,  particularly 
on  the  first  low  strain  cycle.  As  previously  mentioned,  this  data  had  to 
be  reworked . 


modulus 
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Figure  '0.  Master  Constant-Rate  Secant  Modulus  Data  for  UTP-3001 -750/77C8 

JANNAF  Specimens 
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The  rework  of  the  data  consisted  of  estimating  maximum  and  minimum  stress 
and  strain  points  that  were  not  detected  by  the  digitized  data  acquisition 
system.  The  data  reduction  system  computed  strain  from  crosshead  travel; 
this  was  satisfactory  except  at  and  below  zero  stress.  The  sample  linkage 
attachment  was  such  that  specimens  would  not  be  put  into  compression  with  the 
exception  of  test  No.  11,  which  was  run  in  an  Instron  with  rigid  clamp  jaws.  Th 
actual  propellar.t  strain  decay  was  estimated  from  other  tests  where  strain 
recovery  was  monitored  by  cathetometer  measurements  for  the  part  of  the  tests  at 
zero  stress  (i.e.,  no  load  on  the  samples).  The  5-in. /min  crosshead  test  for 
UTP-19,360B  was  selected  as  typical  and  is  shown  in  Figure  11. 

3.1.6  Cyclic  Loading  Test  No.  7 

Cyclic  loading  tests  were  run  on  6-in.  bars  of  UTP-3001  and  UTP-19,360B 
propellants  at  ambient  temperature.  The  cycling  was  for  20  cycles  at  nominal 
strain  levels  of  H ,  8,  and  12%  for  UTP-19,360B  with  UTP-3001  limited  to  l.’l.  At 
the  end  of  the  test  (after  unloading  to  zero  stress.!  the  strain  was  monitored  01 
the  samples  with  a  cathetometer.  The  test  at  a  nominal  12*  strain  is  slow.  0 


temperature 


Stress  While  Cycling  for  UTP-19 , 360B-400/1777 


Figure  12  for  UTP-19, 360B  as  typical.  These  data  were  modified  to  insert  the 
estimated  maximum  stress  points  and  the  propellant  strain  decay  while  at  zero 

stress . 

3.1.7  Relaxation  Test  No.  8 

Stress  relaxation  tests  are  run  on  6-in.  bars  of  UTP-3001  and  UTP-19, 360B 
propellant  for  a  24-hr  period  and  then  monitored  for  strain  decay  after  being 
unloaded.  The  tests  at  ambient  temperature  were  repeated  for  it,  3,  and  12% 
nominal  strain  levels.  They  were  loaded  at  a  crosshead  rate  of  20  in. /min  and 
unloaded  at  1  in. /min  after  the  24-hr  relaxation.  Typical  test  data  are  shown 
in  Figure  13. 

3.1.8  Predaraaged-Relaxation  Test  No.  9 

The  predamaged-relaxation  tests  were  run  with  6-in.  bars  on  UTP-3001  and 
UTP-19.360B  propellants  at  ambient  temperature.  They  were  preloaded  to  12? 
and  unloaded  at  a  crosshead  rate  of  0.1  in. /min  allowed  to  rest,  then  reloaded 
to  8  or  4?  strain  at  20  in. /min.  Af'ei  relaxing  1  hr,  the  samples  were  unloaded 
at  1.0  in. /min  and  strain  was  monitored  after  unloading.  These  tests  were 
repeated  for  a  6%  predamage  strain  followed  by  a  reload  to  4  or  2%  as  above . 

These  data  were  modified  to  obtain  the  peak  stress  and  strain  relaxation 
after  the  samples  were  unloaded.  Strain  was  monitored  with  a  cathetometer 
after  the  relaxation  part  of  the  test.  Typical  test  data  are  shown  in  Figure  14 
for  UTP-19 ,360B. 

3.1.9  Complex  Multiple  Load  Test  No.  10 

The  complex  multiple  load  tests  were  run  with  6-in.  bars  on  UTP-3001  and 
UTP-19, 360B  propellants.  Tests  were  run  at  crosshead  rates  of  5,  1,  and 
0.1  in. /min.  The  test  sequence  was  12  to  8  to  12  to  4?  strain,  then  unloaded, 
reloaded  to  4$  strain,  and  unloaded  (four  cycles)  with  cathetometer  monitoring 
of  strain  decay  on  the  last  unload.  The  same  type  of  sequence  was  repeated 
with  maximum  strains  of  8  and  then  4?  where  the  4?  maximum  strain  was  shortened 
by  one  cycle.  The  5  in. /min,  12?  maximum  strain  test  with  UTP-19, 360B  propellant 
is  typical  and  is  shown  in  Figure  15.  The  data  were  reworked  to  obtain  maximum 
and  minimum  stress  values  as  well  as  strain  decay  for  unloaded  specimens. 

The  cathetometer  strain  after  the  final  unload  was  incorporated  into  the  data. 
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Test  No.  8  -  Stress  While  Step  Straining  for  UTP-800 1-7^0/7768 
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Figure  HI.  Test  No.  9  -  Stress  While  Step  Straining  for  UTP-19 , 760B-400/1777 


temperature 


The  6-in. /min,  12$  strain  test  for  UTP-3001  was  deleted  because  the  modulus 
was  outside  the  range  of  the  remainder  of  the  data  (i.e.,  carton-to-carton 
difference) . 

3.1.10  Complex  History  Test  No.  11 

This  particular  test  is  very  complex  and  required  two  full  months  to 
complete.  It  was  run  with  a  single  sample  (one  at  a  time)  in  an  Instron  using 
single  6-in.  bar  specimens  of  UTP-3001  and  UTP-19.360B  propellant.  The  samples 
were  pinned  through  the  redwood  end  tabs  and  clamped  onto  the  pins  to  avoid 
crushing  the  redwood  (i.e.,  that  would  generate  a  compression  load  on  the 
sample)  during  attachment.  All  coupling  joints  were  heavy  and  tightly  pinned  so 
that  the  sample  would  be  put  into  compression  when  returned  to  the  zero  strain 
position . 

Only  the  test  on  UTP-3001  is  selected  for  presentation  here,  but  it  is  very 
similar  to  the  UTP-19,360B  test.  Since  the  test  is  so  complex,  it  has  been 
divided  into  three  parts.  Part  1  is  described  in  Table  1.  Since  an  actual 
Instron  trac<=  of  the  load-time  curve  was  obtained,  the  peaks  and  minimum 
(compression)  stresses  were  selected  data  reduction  points.  Part  1  of  the  test 
is  expanded  in  time  scale  to  show  some  detail  of  the  process  (Figures  16  through 

19).  Part  2  is  described  in  Table  2.  Test  sequence  for  Part  2  is  shown  in 

Figures  20  through  22.  Part  3  (selected  cycle  maximum  and  minimum)  are  given  in 

Table  9.  "Hie  last  figure  of  this  part  (Figure  23)  is  some  of  the  cyclic  loading 

at  the  end  of  the  test.  The  chart  speed  was  set  such  that  good  definition  of 
the  cycle  could  be  recorded. 

The  latter  part  of  the  cycling  represented  only  by  the  maximum  and  minimum 
stress-strain  points.  Part  3  of  this  test  (the  balance  of  cycling  to  failure) 
is  recorded  in  Table  3  as  maxima  and  minima  for  selected  cycles  sufficiently 
close  to  describe  the  upper  and  lower  bounds.  A  plot  of  the  data  would  be 
similar  to  Figure  22.  The  strain  values  in  Table  3  are  stable,  while  the 
maximum  stress  shows  a  continual  decay,  and  compressive  (negative  values) 
stresses  are  less  compressive. 


TABLE  1.  TEST  NO.  11,  PART  1  -  QUINLAN  COMPLEX  HISTORY  FOR  UTP-3001 


Test 

No. 

Cycle 

Rate , 
in. /min 

Remarks 

1  to 

Load 

2 

Approximately  15-min  rest  after  cycle 

7 

Unload 

2 

Approximately  15-min  rest  after  cycle 

8 

Load 

1 

Approximately  15-min  rest  after  cycle 

Unload 

1 

9 

Load 

5 

Approximately  15-min  rest  after  cycle 

Unload 

5 

10 

Load 

0.5 

Approximately  15-min  rest  after  cycle 

Unload 

0.5 

11 

Load 

10 

Approximately  15-min  rest  after  cycle 

Unload 

10 

12 

Load 

0.2 

Approximately  15-min  rest  after  cycle 

Unload 

0.2 

13 

Load 

2 

Relax 

1/2  hr 

Unload 

2 

Approximately  30-min  rest  after  cycle 

14 

Load 

2 

Relax 

1  hr 

Unload 

2 

4 -day  rest  after  cycle 

15 

Load 

2 

Relax 

1  hr 

Unload 

2 

7-day  rest  after  cycle 

3.2  TWO-DIMENSIONAL  AND  VARIABLE  TEMPERATURE  INVESTIGATION 

The  biaxial  and  nonisothermal  testing  was  conducted  on  specimens  of 
UTP-3001  and  UTP-19,360B  propellants  as  detailed  in  Figure  23.  The  biaxial 
samples  were  cast  into  prelined  redwood  boxes  with  a  1.25-in.  gage  length 
by  6-in.  wide  and  machined  flat  to  a  0.25-in.  thickness.  The  response  prope»'t  i< 
rather  than  failure  properties  were  of  interest  so  the  discontinuity  at  the 
redwood  interface  did  not  affect  the  desired  behavior.  The  1/.'-  x  1/2-  x  b-m. 
specimens  were  used  for  st.rn  ininp.-cool  ing  and  cyclic  strain  tests.  Shear 


Nominal  temperature 


Figure  Id.  Test  No.  11,  Part  1  -  Stress  While  Cycling  for  UTP-3001-750/7768 


w 


mu 


Test 

No. 

Cycle 

Rate, 
in. /min 

Remarks 

16 

Load 

Unload 

■  1  B 

Approximately  30-min  rest  after  cycle 

17 

Load 

0.02 

Unload 

0.02 

Approximately  30-min  rest  after  cycle 

18 

Load 

0.05 

Relax 

3  hr 

Unload 

0.05 

2-weeks  rest  after  cycle 

19 

Load 

0.02 

Unload 

0.02 

Approximately  30-min  rest  after  cycle 

20 

Load 

0.02 

Unload 

0.02 

Approximately  30-min  rest  after  cycle 

21 

Load 

0.05 

Relax 

3  hr 

Unload 

0.05 

1 -month  rest  after  cycle 

22  to 

Cycling 

5 

Several  cycles  monitored,  followed 

42 

by  several  with  only  maximum  and 

minimum  recorded 

relaxation  tests  were  run  with  1-  x  1-  x  3-in.  specimens  bonded  directly  to 
steel  anvils.  Details  are  given  in  later  sections,  except  for  test  21  in  Fig¬ 
ure  23  which  was  uniaxial.  The  six-channel  test  equipment  discussed  in  subsec¬ 
tion  3.1  was  used  for  this  testing;  however,  only  three  biaxial  specimens  could 
be  tested  at  once  because  of  space  limitations  in  the  conditioning  boxes. 


The  biaxial  specimens  used  in  this  part  of  the  program  were  cast  into 
redwood  boxes  similar  to  that  shown  in  Figure  5.  The  space  between  redwood 
blocks  was  1.25-in.  instead  of  the  6-in.  for  the  uniaxial  bars.  The  mill  fin¬ 
ished  geometry  is  shown  in  Figure  24.  The  propellant  was  left  flat  rather 
than  necking  it  down  as  is  done  with  standard  JANNAF  biaxial  specimens.  The 
gage  length  was  designated  as  the  wood-to-wood  distance  for  strain  evaluation. 


Nominal  temperature  =  70  F 


Stress  While  Cvcline  for  UTP-7001-750/7768 


Time 


Strain 


Stress 


Remarks 


74388.0479 

74390.3479 
74390.4539 

74391 .3819 
74391.4239 

74391.9432 

74391.9859 

74392.0799 

74392.1219 

74393.0099 

74393.0519 

74394.7799 

74394.8219 

74395.7279 

74395.7479 

74418.3479 
74418.3979 

74451 .5979 
74451.8779 

74474.2779 

74474.3179 

74698.0179 


74928.0179 


75158.0179 


75388.0179 


mm 

-6.37 

End  of  cycle  20 

93.23 

Peak  of  cycle  30 

mmi 

-5.98 

End  of  cycle  30 

12.17 

92.03 

Peak  of  cycle  34 

2.40 

-5.58 

End  of  cycle  34 

12.17 

91.24 

Peak  of  cycle  40 

2.40 

-5.58 

End  of  cycle  40 

12.17 

89.64 

Peak  of  cycle  50 

2.40 

-5.18 

End  of  cycle  50 

12.17 

88 . 45 

Peak  of  cycle  60 

2.40 

-5.18 

End  of  cycle  60 

12.17 

84.86 

Peak  of  cycle  80 

2.40 

-4.78 

End  of  cycle  80 

12.17 

80.10 

Peak  of  cycle  102 

2.40 

-4.78 

End  of  cycle  102 

12.17 

76.49 

Peak  of  cycle  164 

2.40 

-4.78 

End  of  cycle  164 

12.17 

71.71 

Peak  of  cycle  330 

2.40 

-4.38 

End  of  cycle  330 

12.17 

68.92 

Peak  of  cycle  442 

2.40 

-3.98 

End  of  cycle  442 

12.17 

62.95 

Peak  of  cycle  980 

2.40 

-3.59 

End  of  cycle  980 

12.17 

60.16 

Peak  of  cycle  1980 

2.40 

-3.59 

End  of  cycle  1980 

12.17 

56.97 

Peak  of  cycle  2980 

2.40 

-3.59 

End  of  cycle  2980 

12.17 

56.18 

Peak  of  cycle  3980 

2.40 

-3.59 

End  of  cycle  3980 

TABLE  3.  TEST  NO.  11,  PART  3  -  UTP-300 1-750/7768  1/2-IN.  BAR  STRESS 

WHILE  CYCLING 

(Sheet  2  of  2)  T8715 


Time 

Strain 

Stress 

Remarks 

75618.0179 

12.17 

53.39 

Peak  of  cycle  4980 

2.40 

-3.59 

End  of  cycle  4980 

75695.5179 

12.17 

51.39 

Peak  of  cycle  5317 

2.40 

-3.59 

End  of  cycle  5317* 

*  Sample  broke 

Ten 

Test 

Description 

Oamaga  Cyclt/Ttsl 

1 

Strain  Cycle 

14 

Biaxial 

constant  rate 

Biaxial  samples  ot  UTP-3001  and  UTM9.360B  were  ramp  loaded  to 
failure  at  rates  ol  2.  0  2  and  0  02  in  /min  at  temperalu'es  oi 

41.  70.  and  1 20  *F 

1 

( 

1 

1 

1b 

16 

Biaxial 

slraining> 

cooling 

Biaxial 

relaxation 

Biaxial  samples  ol  UTP-3001  and  Utr'  H.360B  were  simultaneously 
stram  and  cooled  from  120  to  40®F  c  e  a  40  hr  period 

Riaxial  samples  ol  UTP-3001  and  UTP-19.360B  were  run  m  stress 
u-'axation  tests  at  40.  70.  and  12Q°F 

gga 

SB 

17 

18 

Shear 

relaxation 

6-in.  bar 
straining-  cooling 

Shear  samples  ot  UTP-3001  and  UTP-19.360B  we'e  run  m  sness 
relaxation  tests  at  70  °F 

6-m  bars  ol  UTP-3001  and  UTP-  19.360B  were  simultaneously  strain 
and  cooled  Irom  120  to  40  °F  at  three  slow  rates 

1 

19 

'Biaxial  Quintan 
complex  history 

Biaxial  samples  ol  UTP-3001  and  UTP-  19.360B  were  cycled  lor  the 
Quinlan  complex  history  test  at  70 'F 

c 

MM 

] 

20 

6-m  bar  cyclic 
lest 

6-m  bars  ol  UTP-3001  and  UTP-19.360B  were  run  in  cytln  strain  tests 
at  0  1  in. /min  and  70°F 

■■FTTT'~TTTflH 

!!,™riiui^rrrTr''n 

21 

Biaxial 

thermal 

similitude 

Biaxial  samples  ol  UTP-3001  and  UTP-19.360B  were  run  in  ramp 
relaxation- ramp  tests  with  simultaneous  cooling  oi  heating  n  e  lor 
reverse  ramp) 

81 

1 

Note.  Nominally  three  samples  were  run  for  each  test  and  condition 

legeno 

T  -  temiwralim- 
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Figure  23.  Biaxial  and  Nonisothernal  Phase  III  Testing 
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Figure  24.  Finished  Biaxial  Specimen 
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Response  properties  rather  than  failure  properties  were  of  interest,  so  the 
boundary  perturbation  was  neglected. 

"The  shear  samples  (test  No.  17)  were  1-  x  1-  x  3-in.  blocks  of  propellant 
•that  were  bonded  to  the  test  fixture  (shown  in  Figure  25)  after  being  machined. 
The  pull  rods  were  attached  to  the  offset  plates  so  that  the  load  was  trans¬ 
mitted  through  the  center  of  the  sample  as  shown.  Since  strain  was  limited 
to  5%  for  the  shear  relaxation  test,  the  sample  was  assumed  to  be  in  simple 
3hear.  The  shear  strain  (7)  was  calculated  as  the  tangent  of  the  displacement 
angle  orAL/G.L.  The  shear  stress  (t)  was  calculated  as  force/area  (arcs 
=  3  in.2). 


The  data  modification  to  insert  peak  and  minimum  stress  points  previously 
discussed  for  uniaxial  cyclic  tests  was  also  used  for  the  biaxial  and 
nonl30thermal  tests. 
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a.  Propellant  and  Loading  Fixture 


b.  Propellant 


Figure  25.  Shear  Sample  and  Test  Attachment 
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3.2.1  Biaxial  Constant  Rate  Test  No.  14 

The  biaxial  constant  rate  tests  to  failure  were  run  with  the 
1/4-  x  1-1/4-  x  6-in.  specimens  of  UTP-3001  and  UTP-19.360B.  The  tests  were 
conducted  at  40,  70,  and  120  F  at  crosshead  rates  of  2,  0.2,  and  0.02  in. /min 
The  typical  stress-strain  response  is  shown  in  Figure  26  for  UTP-19,360B 


at  71  F  and  2  in. /min.  Because  of  the  fixtures  and  more  difficulty  in  adjust 


a  f^vws^ vi  a  r, 


Figure  ?6.  Test  No.  14  -  Stress  for  UTP-19 ,  ^>0n-'ino/1777 


linkage  than  with  the  6-in.  bar  specimens,  the  three  samples  did  not  load 
simultaneously.  Consequently,  it  was  necessary  to  adjust  the  data  for  each 
sample  to  have  a  common  starting  point. 


3.2.2  Biaxial  Straining-Cooling  Test  No.  15 

Biaxial  specimens  of  UTP-3001  and  UTP-19,360B  propellants  were 
simultaneously  strained  and  cooled  from  115  to  40  F  at  a  crosshead  rate  of 
approximately  3  x  10~5  in. /min  over  a  40-hr  period.  Typical  test  results  are 
shown  in  Figure  27  for  UTP-3001.  The  stress-time  traces  for  all  three  samples 
appear  to  start  together  but  spread  out  as  the  test  progressed. 

3.2.3  Biaxial  Stress  Relaxation  Test  No.  16 

Biaxial  stress  relaxation  tests  were  run  with  the  1/4-  x  1-1/4-  x  6-in. 
specimens  of  UTP-3001  and  UTP-19,360B  propellants  at  a  nominal  3%  strain  and  teim 
peratures  of  40,  70,  and  120  F.  The  loading  ramp  rate  was  0.2  in. /min.  Typical 
relaxation  data  for  UTP-19,360B  at  40  F  are  shown  in  Figure  28. 

3.2.4  Shear  Relaxation  Test  No.  17 

Shear  relaxation  tests  were  run  on  1-  x  1-  x  3-in.  samples  of  UTP-3001  and 
UTP-19.360B  prope1  ..ants.  The  samples  were  post-bonded  to  steel  plates  as  shown 
in  Figure  25  and  run  one  at  a  tin  by  loading  them  at  0.2  in. /min  and  ambient 
temperature  with  offset  fixtures  so  the  load  was  transmitted  through  the 
centerline  of  the  sample.  The  three  samples  for  each  propellant  were  hand- 
reduced  and  digitized  for  computer  storage  and  printout.  Typical  test  data  are 
presented  in  Figure  29  for  UTP-3001. 

Peak  stresses  and  strain  were  very  close,  as  was  the  1  hr  relaxation 
stress  on  each  propellant,  even  though  the  samples  were  run  separately. 

3.2.5  Uniaxial  Straining-Cooling  Multiple  Rates  Test  No.  18 

The  rate  effect  on  the  straining-cooling  response  was  determined  on  UTP- 
3001  and  UTP-19.360B  propellants.  The  1/2-  x  1/2-  x  6-in.  uniaxial  bar  sample 
wsa  used  so  that  testing  could  be  completed  in  the  shortest  time  possible.  The 
rate  effect  for  the  uniaxial  specimens  was  then  applied  to  biaxial  test  No.  15. 
Cooling  was  from  110  to  40  F  at  the  crosshead  rates  of  0.002,  0.0002,  and 
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temperature 


Stress  While  Step  Straining  for  UTP-19 , 360B-,)00/ 177 


Nominal  temperature 


0.0004  in. /min.  Typical  data  a~e  shown  in  Figures  30  for  UTP-3001.  Good 
reproducibility  is  shown  within  the  set  of  three  samples. 

3.2.6  Biaxial  Quinlan  Complex  History  Test  No.  19 

The  1/4-  x  1-1/4-  x  6-in.  biaxial  samples  of  UTP-3001  and  UTP-19.360B  were 
subjected  to  the  complex  cycling  and  relaxation  history  indicated  in  Figure  23. 
When  the  tests  were  run,  the  linkage  (misalignment,  etc.)  was  such  that  the 
samples  were  not  loaded  with  an  equivalent  amount.  The  first  sample  to  be 
loaded  had  the  correct  strain  determined  but  strain  on  the  other  samples  had  to 
be  adjusted  to  the  time  the  ramp  started  on  each.  Each  sample  was  separately 
reduced  and  data  were  modified  to  pick  up  the  peak  and  minimum  stress  points. 
Sample  1  for  UTP-3001  is  shown  in  Figures  31  through  33  where  the  complex  test 
has  been  divided  into  segments  on  an  expanded  time  scale  to  show  test  details. 
The  first  cycle  in  Figure  31  shows  no  load  and  the  second  cycle  shows  very 
little.  By  contrast  sample  3  (not  included  here)  had  a  first  peak  stress-strain 
of  37  psi,  2.26*  and  a  second  peak  of  76  psi  and  4.97*.  During  the  unload  part 
of  the  cycle  after  the  stress  reached  zero,  strain  decay  was  estimated  from 
other  tests  in  which  the  strain  was  mea'  ured  using  a  cathetometer . 

3.2.7  Biaxial  Ramp-Relax-Ramp  Test  No.  21 

Ramp-relaxation-ramp  tests  We..  e  run  on  UTP-3001  and  UTP-19,360B  propellants 
with  the  1/4-  x  1-1/4-  x  6-in.  biaxial  specimens.  The  first  test  was  ramped  at 
0.0005  in. /min  to  6?  strain  and  simultaneously  cooled  from  120  to  70F.  It  was 
held  at  6*  strain  nearly  23  hr,  then  ramp-loaded  to  failure,  while  cooling 
towards  40  F.  Data  for  UTP-3001  are  shown  in  Figure  34;  those  for  UTP-19.360B 
were  similar.  The  cooling  cycle  did  not  end  at  the  peak  strain;  consequently, 
the  relaxation  of  stress  was  not  the  normal  type  behavior.  The  continued 
cooling  increased  the  propellant  stress  so  that  the  normal  relaxation  behavior 
did  not  start  until  the  propellant  temperature  stabilized. 

This  test  was  repeated  starting  at  110  F  and  taken  to  6*  strain  with  a  peak 
stress  of  70  psi  compared  to  30  psi  for  the  above  test.  The  longer  ramp  time 
allowed  the  cooling  to  reach  40  F  at  the  peak  stress.  The  samples  were  allowed 


Fieure  70.  Test  No.  18  -  Stress  While  Straining  and  Cooling  for  UTP-3001-750/7768 
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Figure  3?.  Teat  No.  19,  Part  2  -  Stress  for  UTP-3001-750/7768 
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Figure  34.  Teat  No.  2 1  -  Stress  While  Complex  Straining  and  Cooling 

for  UTP-300 1-750/7768 


to  relax  overnight  and  then  unloaded  to  3%  while  warming  the  samples  to  room 
temperature . 

3.2.8  Uniaxial,  Biaxial,  and  Shear  Comparison 

Comparisons  of  the  different  samples  were  made  in  order  to  show  that  the  pro¬ 
pellant  used  in  each  of  the  tests  was  of  the  same  family.  These  comparisons 
were  done  at  ambient  temperatures  for  selected  rates  and  were  limited  to  strain 
rates  that  were  close  to  each  other.  Such  testing  minimized  the  time- 
temperature  equivalence  shifts  to  small  changes  for  negligible  data  input 
errors.  The  adjustments  made  for  strain  levels  are  given  in  Table 

The  comparisons  between  uniaxial  and  biaxial  in  the  table  are  close  to  the 
theoretical  ratio  of  75  to  80$.  The  shear  to  uniaxial  ratios  of  0.28  and  0.39 
bracket  the  nominal  theoretical  value  of  1/3.  For  the  comparisons  ma^o  in 
Tab_e  M,  the  UTP-3001  and  UTP~19,360B  propellants  have  to  be  considered  part  of 
the  same  family.  Any  minor  differences  can  be  attributed  to  carton-to-carton 
variations. 

Uniaxial  and  biaxial  stress  relaxation  modulus  data  for  UTP-19.360B  are 
compared  in  Figure  35.  The  biaxial  modulus  was  shifted  using  the  uniaxial  AT 
shift  factor.  Figure  35  also  shows  that  excellent  agreement  is  obtained  when 
the  uniaxial  modulus  is  multiplied  by  the  theoretical  value  of  *4/3  to  correct  it 
to  a  2:1  principal  state  of  stress. 


4.0  THEORETICAL  DEVELOPMENT 


4.1  INTRODUCTION  AND  PRELIMINARY  STUDY 

Generally  speaking,  solid  propellants  may  be  considered  as  lightly 
crosslinked  long-chain  polymers,  highly  filled  with  coarse  solid  particles. 

They  respond  viscoelastically  to  the  action  of  external  stimuli.  But  certain 
aspects  of  their  behavior  cannot  be  reproduced  by  classical  linear  or  nonlinear 
theories  of  fading-memory  materials.  Thus,  in  recent  years,  much  work  has  been 
concerned  with  the  development  of  appropriate  models  to  predict  the  mechanical 
response  of  solid  propellants.  A  current  trend  is  to  express  the  observed 
response  in  terms  of  some  measure  of  ’'damage"  at  the  continuum  level  where 
damage  is  described  as  the  difference  between  the  observed  response  and  that 
predicted  by  a  fading- theory,  such  as  linear  viscoelasticity.  There  is  now 
sufficient  experimental  evidence  to  show  that  damage  per  se  is  a  microscopic 
phenomenon  associated  with  the  initiation  and  growth  of  flaws,  debonding  between 
matrix  and  solid  filler  particles,  and  molecular  chain  scission.  Although  it  is 
largely  irreversible,  damage  is  partially  recoverable  shortly  after  removal  of 
the  loading  system.  This  recovery  from  damage  is  termed  as  "healing."  It  is 
clear  that,  depending  on  the  propellant  and  service  requirements,  it  may  also 
have  to  be  accounted  for  in  a  constitutive  theory  for  solid  propellants. 

In  the  present  program,  two  approaches  to  characterizing  damage  have  been 
followed.  In  the  first  one,  damage  is  treated  as  the  algebraic  difference 
between  the  measured  stress  and  that  predicted  by  linear  viscoelasticity,  so 
that: 


Oq  (t)  =  a(t)  -  og  (t) 


(1) 


in  which  og  and  oc  are  the  linear-viscoelastic  and  correction  terms,  respectively, 
with  a ,  the  measured  stress.  In  the  second  approach,  the  difference  between 
measured  and  fading-memory  type  stresses  is  handled  through  a  stress-correction 
function  in  the  following  form: 


I 


The  softening  function  (C)  is  made  to  depend  on  the  past  maximum  strain  or 
stress  and  a^(t)  represents  an  appropriate  function  of  the  fading-memory  type 
stresses . 

Broadly  speaking,  the  models  of  R.  Farris,  M.  Gurtin  and  M.  Quinlan  are  of 
the  type  presented  in  equation  (1)  above,  while  those  of  R.  Schapery,  W.  Hufferd 
and  Swanson  are  of  the  form  given  by  equation  (2). 

The  following  section  presents  some  experimental  evidence  on  the  types  of 
nonlinearities  exhibited  by  solid  propellants,  and  briefly  discusses  the 
pioneering  work  of  Mullins  and  Tobin  in  treating  the  large  hysteresis  observed 
in  tire  rubbers.  During  the  preliminary  phase  a  brief  evaluation  of  the 
potential  of  basing  constitutive  models  on  molecular  dynamics  considerations  was 
reviewed  by  Dr.  S.  M.  Breitling,  subcontractor  to  CSD,  and  Dr.  R.  P.  Wool  of  the 
University  of  Illinois.  Their  conclusions  and  recommendations  are  discussed 
next.  Subsequently,  the  nonlinear  theory  of  Farris 5  is  presented.  This 
theory  was  used  during  the  first  pha_ .  of  the  program  to  predict  the 
response  of  TP-H1011  and  to  compare  wit.*  preliminary  results  of  the  other  five 
constitutive  laws  investigated. 

Next,  the  theory  of  linear  viscoelasticity  is  applied  to  predict  the 
response  of  UTP-19.360B  and  UTP-300'i  under  various  strain  histories.  The 
ensuing  results  are  the  basis  for  comparing  the  propellant  response  as  predicted 
by  each  of  the  constitutive  laws.  This  comparison  is  most  meaningful  because 
each  of  the  theories  considered  evolved  from  a  set  of  modifications  to  linear 
viscoelasticity. 

Finally,  a  detailed  description  is  given  of  each  of  these  five  stress-strain 
relations.  These  Include  the  original  concept  of  the  models,  their  current 
versions,  comparisons  of  predicted  and  measured  stresses  for  a  variety  of  strain 
histories,  and  some  pertinent  guidelines  for  characterizing  solid  propellant 
according  to  each  theory. 

4,1.1  Experimental  Background 

The  complex  behavior  of  solid  propellants,  as  well  as  some  attempts  at 
developing  usable  stress-strain  laws  for  these  materials,  are  well  documented  in 
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numerous  studies. 3-31  These  studies  show  that  a  given  deformation  process 
causes  a  change  in  the  response  properties  of  solid  propellants,  for  instance  a 
drop  in  the  relaxation  modulus.  An  stated  before,  this  deviation  from  some 
expected  response  is  what  has  been  called  damage.  It  is  evidenced  as 
phenomenological  microscopic  changes  that  are  caused  by  undefined,  but  real, 
irreversible  or  partially  reversible  microscopic  changes.  Polymer  bend 
breakage,  vacuole  formation  in  the  polymer  matrix,  dewetting  between  the  polymer 
matrix  and  solid  filler  particles,  microcracking,  etc.,  are  among  the  possible 
microscopic  causes  of  observed  permanent-memory  effects  in  propellants. 

Studies  on  uniaxial  solid  propellant  samples  have  indicated  that  these  mate¬ 
rials  exhibit  large  hysteresis  even  at  small  strains.  These  studies  have  also 
revealed  that  the  state  of  damage  in  solid  propellants  is  determined  primarily 
by  the  maximum  strain  or  stress  undergone  during  the  loading  histories. 

The  typical  nonlinear  hysteresis  and  permanent-memory  effects  exhibited  by 
solid  propellants  are  illustrated  in  Figure  36.  A  series  of  finite-duration, 
variable-strain-level  ramp  pulses  were  used  to  obtain  the  propellant  response 
subsequent  to  a  given  damage  history. ^2>  1 3  All  ramps  had  the  same  initial 
moderate  rise  rate,  with  two  exceptions  to  be  noted  later,  and  all  ramps 
had  the  same  very  slow  decline  rate. 

Observations  of  the  load  on  the  specimen  after  returning  to  its  original 
length  (zero  strain)  showed  that  it  took  approximately  30  min  for  the  stress  to 
relax  to  zero. 

A  series  of  tests  was  run  on  a  1/4-  x  1/4-  x  4-in.  tab-end  sample.  The 
virgin  specimen  was  initially  strained  to  a  level  of  7.04J  and  allowed  to  relax 
to  achieve  a  rest-state  condition.  The  first  part  of  the  testing  is  shown  in 
the  lower  half  of  Figure  36  (curves  A-H)  and  the  last  part  in  the  uppe-  halr 
(curves  H-M) . 

Curve  A  shows  the  load  response  to  the  first  pulse.  The  specimen  was  then 
subjected  to  four  successive  ramp  strain  pulses  ranging  from  e°  =  3.8?|  t  t°  = 
6.34$.  There  was  a  rest  period  allowed  between  each  pulse.  The  resjl*  - 
shown  in  Figure  36  as  curves  B  through  E. 
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Figure  36.  Relaxation  After  Damage 
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Two  aspects  of  the  propellant’s  behavior  are  worth  noting.  First,  during 
the  constant  strain  rate  portion  of  the  ramp,  each  successive  load-time  curve  is 
essentially  identical.  This  phenomenon  indicates  that  the  "new  material"  has 
the  sane  nonlinear  rate-dependency  under  repeated  strain  conditions  as  long  as 
the  strain  levels  are  below  the  previous  maximum  strain  of  c0=  7.0UJ.  Second, 


the  relaxation  portions  of  the  curves  are  essentially  homologous,  indicating 
that  a  viscoelastic  relaxation  process  is  taking  place. 

Curves  F  and  G  present  the  results  of  two  additional  tests  at  two 
successively  lower  strain  rates  where  the  sample  was  loaded  to  5.64$  strain  each 
time.  A  strong  rate  dependency  is  observed  during  the  rise  portion  of  the  ramp; 
however,  curves  F  and  G  rapidly  rejoin  curve  D  indicating  that  the  material  is 
behaving  in  a  viscoelastic  fading-memory  fashion. 

The  specimen  was  next  subjected  to  a  ramp  strain  pulse  reaching  a  higher 
strain  level  (e°  =  9.86$)  than  the  maximum  7.04$  previously  experienced  (Fig¬ 
ure  36,  curve  H).  The  first  part  of  curve  H  repeats  the  loading  ramp  portion  of 
curves  B-E  to  indicate  the  same  "new  material"  rest-state.  Note  that  the  load¬ 
time  curve  returns  to  the  initial  or  virgin  constant  strain  rate  curve  once  the 
previous  maximum  strain  (7.04$)  has  been  passed. 

Subsequently,  the  specimen  was  strained  with  the  ramp  pulse  to  four 
different  strain  levels  less  than  9.86$  (f°  =  4.26$,  5.64$,  7.04$,  and  9.01$), 
as  shown  in  curves  I  through  L.  The  results  show  that  a  new  rate-dependency  has 
developed  (compare  the  rising  portions  of  I  through  L  with  the  rising  portion  of 
H).  Thus,  another  "new-material"  rest-state  has  been  produced  as  a  result  of 
the  second  maximum  strain  level  of  9.86$.  Lastly,  the  specimen  was  strained  to 
another  new  maximum  of  e 0  =  12.68$  as  shown  in  curve  M.  It  returned  again  to 
the  virgin  undamaged  curve  once  the  9.86$  strain  level  was  exceeded. 

The  above  experimental  evidence  suggests  that  the  form  of  the  constitutive 
equation  should  remain  unchanged  with  respect  to  the  material's  current  rest- 
state.  This  condition  should  remain  as  long  as  the  damage  is  unchanged  (i.e., 
the  e  max  is  unchanged  during  its  subsequent  strain  histories) . 

Figure  37  shows  a  replot  (curves  N  and  0)  of  some  of  the  results  just 
discussed.  After  an  initial  maximum  strain  (7.04$)  the  specimen  was  allowed 
to  return  to  a  rest-state  and  then  strained  to  a  value  of  e°  =  5.64$,  with 
the  result  shown  as  curve  P.  These  three  identical  strain  history  tests 
of  three  different  material  states  indicate  that  the  higher  the  state  of  damage 
(primarily  enax),  the  softer  the  material  response  upon  subsequent  testing. 
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Figure  37.  Relaxation  after  Damage 
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In  addition,  other  experimental  studies  have  pointed  out  the  importance  of 
healing  effects,  load  duration,  and  initial  strain  rate.12'  ^  Finally, 

It  Is  important  to  note  that  the  behavior  of  solid  propellant,  depicted  in 
Figures  36  and  37,  is  similar  to  that  exhibited  by  rubber.  The  nonlinear 
uniaxial  stress  response  of  rubber,  with  or  without  carbon  black  filler  and 
1q  the  absence  of  time  effects,  was  well  characterized  by  Mullins  and  Tobin^7 
with  equation  (3). 


e  *  eu  F 


where : 


c  =  engineering  strain.  The  Mu  11 ins -Tobin  model  is  not  limited  to  small 
strains . 

eu=  <  u(o)  =  strain  as  a  function  of  engineering  stress  for  the  polymer 
without  filler  and  without  damage.  The  characteristic  shape  of  this 
function  is  shown  in  Figure  38. 

F=  F(omax,  N)  =  damage  or  softening  function  which  depends  on  the  maximum 
stress  experienced  by  the  rubber  and  the  number  N  of  loading  and 
unloading  cycles.  F  is  not  very  sensitive  to  N,  but  depends  strongly  on 
any  hard  filler  particles  that  may  be  present. 

A  large  amount  of  rubber  data  can  be  predicted  by  means  of  this  equation 
when  the  samples  are  not  allowed  to  rest  between  cycles.  Recovery  or  healing 
occurs  as  a  function  of  the  rest  time.  Therefore,  healing  would  have  to  be 
considered  in  an  accurate  characterization  of  rubber. 

Introducing  the  inverse  of  £u  =  eu(o),  equation  (3)  may  be  put  in  the  form: 

(4) 

c  =  f  (e/F) 


which  shows  that  F  (where  F<1 )  is  a  strain-magnification  factor.  The  ratio  e/F 
is  interpreted  by  Mullins  and  Tobin  to  be  the  average  strain  in  the  rubber  phase 


Figure  38.  Stress-Strain  Curve  for  Rubber 
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of  a  hard  particle-filled  rubber.  Without  damage  in  a  highly-filled  rubber,  F<C1 , 
A a  the  rubber  is  cycled  between  the  strains  e  =  0  and  e  =  e  the  ratio  e/F  at 
any  strain  decreases,  and  therefore  the  stress  decreases.  The  shape  of  the 
stress-strain  curve  is  still  as  shown  in  Figure  38.  It  is  similar  to  that 
for  solid  propellant  after  first-time  loading.  This  fact  and  the  ability 
of  the  model  presented  in  equation  (3)  to  reproduce  a  large  amount  of  rubber 
data  explain  the  great  influence  of  the  Mullins-Tobin  approach  on  the  development 
of  nonlinear  constitutive  theories  for  solid  propellants. 

4.1.2  Potential  Physical  Parameters  for  Damage  and  Healing 

A  fundamental  problem  of  continuum  mechanics  is  the  rationalization  of  mathe¬ 
matical  models  with  a  body  of  experimental  data.  As  the  complexity  of  the 
physics  increases  it  becomes  progressively  difficult  to  identify  unambiguously 
those  theories  and  models  which  properly  treat  the  critical  mechanisms  operating 
in  the  physical  process.  Nowhere  is  this  more  apparent  than  in  viscoelastic 
materials  undergoing  damage.  Frequently,  a  theory  is  a  product  of  years 
of  evaluation  with  layer  upon  layer  °  parametric  refinements  that  lead  to 
better  numerical  fits  of  given  data  set-  which  may  not  elucidate  the  underlying 
physics.  It  becomes  very  difficult  to  transfer  the  acquired  experience  with 
a  given  model  to  a  new  material  system  or  a  markedly  different  application 
in  the  same  system. 

Preoccupation  with  chemical-physics  is  equally  fascinating  and  unproductive. 
The  detailed  investigation  of  molecular  and  chemical  properties  has  been  the 
domain  of  many  investigators  for  long  periods  of  time  with  little  information 
that  is  applicable  to  real  problems.  However,  several  trends  have  been 
developing  in  the  last  decade  that  may  prove  productive  if  prudently  implemented 
and  coordinated  with  an  applied  effort.  Boyd  and  associates  at  the  University 
of  Utah  have  successfully  treated  cooperative  molecular  motions  in  crystals  and 
in  the  pseudo-amorphous  state.  The  results  of  this  work  has  not  been  applied  to 
the  development  of  continuum  models. 

A  relevant  question  may  be  asked:  "Is  it  possible  to  apply  molecular 
dynamic  methods  directly  to  continuum  mechanics?"  The  answer  is  clearly  "no". 
Continuum  theories  do  not  consider  the  discrete  molecular  processes,  Just  as  the 
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molecular  dynamics  neglect  microscopic  deformations  and  stresses.  However,  the 
point  of  commonality  between  the  two  approaches  is  the  treatment  of  energy. 

Both  methods  must  consider  energy  as  a  function  of  movement  and  dissipative 
processes.  Thus,  it  may  be  possible  to  use  the  molecular  dynamic  work  done  to 
date  to  guide  in  the  selection  of  mathematical  continuum  models  by  evaluation  of 
model  functional  forms  and  parameters  based  on  molecular  energetics.  We  were 
unable  to  accomplish  this  goal  within  the  resources  available  to  this  program. 
Still,  this  is  a  laudable  objective  and  should  be  continued  in  future  constitu¬ 
tive  theory  development  programs. 

The  evaluation  would  not  be  of  a  formulation  of  models  based  on  first  prin¬ 
ciples  nor  a  calculation  of  model  parameters.  Rather,  both  functional 
form  and  parametric  values  could  be  judged  for  physical  appropriateness  using 
the  molecular  dynamic  information.  Conversely,  the  experimental  continuum  data 
could  eliminate  many  molecular  processes  that  do  not  contribute  significantly  to 
the  macroscopic  viscoelastic-damage  mechanism.  The  interplay  of  a  more  fundamen¬ 
tal  molecular  viewpoint  with  the  mathematical  models  and  an  expanded  experimen¬ 
tal  data  base  could  prove  extremely  useful  not  only  in  the  rapid  identification  of 
the  fruitful,  productive  mathematical  formulations,  but  in  the  elimination  of  poor 
models. 

Included  in  the  area  of  molecular  dynamics  is  rate  theory  applications  to 
macromolecular  mechanical  and  dielectric  loss  experiments,  cooperative  polymer 
motion  energetics,  fracture  and  electron  microscope  fracture  studies.  These 
studies  do  not  represent  a  coherent  single  field  of  study  but  rather  a  wide 
range  of  studies  conducted  on  many  types  of  materials,  most  frequently  unfilled. 
The  initial  effort  in  this  program  was  to  collect  and  review  some  material 
likely  to  be  of  most  interest  and  relevance  to  the  propellant  behavior.  An 
evaluation  should  lead  to  the  identification  of  a  few  molecular  modes  of  motion 
and  associated  energetics  applicable  to  the  binaer.  This  information  can  then  be 
interpreted  in  a  form  (e.g.,  energy  density,  surface  energy)  helpful  to  the 
mathematical  continuum  model.  Implicit  in  this  interpretation  is  the  descrip¬ 
tive  physical  model  of  molecular  behavior  that  gives  rise  to  the  macroscopic 
observations. 
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A  small  cooperative  effort  between  molecular  and  continuum  approaches  was 
carried  out  in  this  program  to  identify  the  compatibility  of  the  mathematics  and 
physics.  Only  minimal  molecular  dynamics  experimentation  was  conducted.  The 
objective  of  this  initial  effort  was  to  capitalize  on  the  molecular  information 
presently  available  and  evaluate  the  productiveness  of  the  initial  approach. 

Table  5  briefly  outlines  potential  areas  where  molecular  dynamics  can  be 
expected  to  make  an  impact  on  damage-healing  phenomena. 

Processes  occurring  at  the  molecular  and  microstructural  levels  control  the 
macroscopic  response  of  materials.  Therefore,  a  knowledge  of  these  events, 
particularly  to  the  extent  of  controlling  damage,  would  provide  guidelines  for 
constructing  the  constitutive  relations  for  the  material.  The  nature  of 
mechanically  induced  damage  can  be  investigated  in  experiments  which  monitor 
both  the  occurrence  of  the  damage  and  disappearance  of  damage  (healing).  Damage 
may  express  itself  mechanically  as  a  stress-softening  effect,  similar  to  the 
"Mullin' s  effect"  in  filled  elastomer  ' ,  and  may  or  may  not  involve  cavitation, 
crazing,  or  bond  rupture.  Damage  can  of‘ an  occur  in  a  time  scale  and  manner  which 
crjeate  difficulty  in  evaluation  or  detection.  However,  healing  studies  provide 
a  method  of  reversing  the  process,  or  "turning  back  the  clock"  in  a  controllable 
(via  temperature)  manner.  Thus,  the  damage  reversibility  concept  provides  a 
method  of  evaluating  damage  mechanisms.  This  consideration  was  investigated 
briefly  and  at  a  preliminary  level  by  Dr.  Richard  Wool  of  the  University  of 
Illinois  during  the  initial  preliminary  studies  of  the  program.  Available 
methods  for  evaluating  crack  or  damage  healing  are  summarized  in  Table  6. 

Damage  healing  studies  were  conducted  with  the  intent  of  learning  the  mecha¬ 
nism  of  material  damage  at  the  molecular  and  microstructural  levels  because 
processes  occurring  at  these  material  levels  control  the  macroscopic  response  of 
the  material.  The  approach  is  to  study  the  nature  of  mechanically  induced  damage 
through  experiments  which  monitor  both  the  occurrence  and  disappearance  (i.e., 
healing)  of  damage.  Dr.  Wool  at  the  University  of  Illinois,  Urbana,  has  con¬ 
ducted  recent  healing  studies  cn  the  polymeric  materials  shown  In  Table  7. 
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TABLE  6.  METHODS  OF  EVALUATING  CRACK  HEALING 


T998? 


Mechanical 

Spectroscopic 

Tensile  testing  (MTS,  Instron) 

Infrared  with  stress  and  tem- 

Dynamic  mechanical  (Rheovibron) 

perature  dispersive  and 

Fourier  transform  (FTIR) 

Izod  impact  fracture 

X-ray,  SAXD,  WAXD 

NMR,  broadline,  pulse,  magic  angle 

Optical 

Eximer  flourescence 

Bright  field  plus  stress  plus 

Electron  microscopy 

temperature 

Dark  field  plus  stress  plus 

Transmission  (TEM)  with  stress 

temperature 

Scanning  (SEM)  with  stress 

Birefringence  plus  stress  plus 

Scanning-transmission  (STEM) 

temperature 

Videorecording 

Other 

Photometric 

Picometry  (density  gradient) 

Small  and  wide  angle  light 

scattering 

Dilatometry 

Visual 

Neutron  scattering 

. -  — . —  — . . — 

TABLE  7.  RECENT  HEALING  STUDIES 
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Materials  Studied 

Damage  Investigated 

Crystalline  polymers 

Microvoids  <1  micron 

Amorphous  polymers 

Crazes  > 1  micron 

Glassy  polymers 

Izod  fracture 

Block  copolymers 

Microvoids 

Filled  composites 

Stress  softening 

Studies  of  these  materials  have  reached  the  following  conclusions: 


•  In  each  material,  complete  healing  was  usually  observed;  i.e.,  all  cracks 
crazes,  and  microvoids  disappear  as  heat  in  the  mechanical  sense. 

•  In  each  case,  the  virgin  mechanical  strength  was  restored. 

•  The  kinetics  of  healing  could  be  excellently  modeled  by  the  following 

equations:  _  _dD_  _  krp  (5) 

dt 

k,  n  constants 


D  =  Damage 


M  +  Kt)a 
K  =  «o  exp  -  E/kT 


K,  aconstants 


(6) 

(7) 


•  There  exists  a  master  healing  curve  (equation  6)  for  each  material 
which  is  obtained  by  similar  WLF  approach  to  viscoelastic  materials. 

The  activation  energy,  E,  is  obtained  from  a  plot  of  log  (shift  factor) 
versus  1/T.  The  master  healing  curve  is  sigmoidal  with  respect  to  log 
t,  with  a  lower  plateau  at  R  =  R0  and  an  upper  plateau  at  R  =  1,  at 
long  times. 


From  these  studies  it  has  been  concluded  that  the  general  mechanisms  of 
healing  consist  of  four  stages: 

1.  Surface  Rearrangement  -  Molecular  motion  on  newly  formed  crack  surfaces 
may  cause  new  configurations  which  affect  the  following  stages.  The 
concept  of  damage  history  is  therefore  important.  If  significant 
rearrangement  occurs,  healing  can  be  prevented. 

2.  Wetting  -  The  surfaces  of  the  crack  can  wet  and  the  crack  disappears 
optically,  but  only  a  minor  fraction  of  the  local  mechanical  strength 
is  restored. 


3.  Diffusion  -  The  chains  on  contacting  surfaces  interdiffuse  and  estab¬ 
lish  the  virgin  configuration  of  interpenetrated  chains  with  a 
subsequent  increase  in  mechanical  strength.  This  step  is  the  slowest 
(rate  determining)  and  is  responsible  for  the  major  fraction  of  mechani¬ 
cal  strength  recovery. 

Randomization  -  The  final  stage  is  subtle  and  involves  long  time  molecu¬ 
lar  motion  which  completely  erases  the  damage  memory. 

From  these  observations,  it  may  be  postulated  that  with  mechanical  loading  the 
mechanisms  of  damage  consist  of: 

•  Derandomization 

•  Dlffusional  demixing 

•  Dewetting 

•  Surface  rearrangement. 

Voids  or  cracks  appear  at  step  3*  The  majority  of  stress  softening  damage, 
however ,  occurs  at  step  2.  Therefore,  considerable  damage  can  occur  under 
conditions  approximating  constant  volume. 

I 

In  the  present  limited  stud> ,  the  objectives  were  to: 

•  Obtain  a  master  healing  curve 

•  Study  matrix  molecular  mechanisms 

•  Determine  voiding  processes 

•  Investigate  molecular  theories  which  may  have  an  impact  on  constitutive 
applications. 

Two  propellants  were  investigated,  an  inert  (NaCl-fill)Thiokol  propellant 
and  a  "live"  PBAN  propellant  supplied  to  Dr.  Wool  by  CSD.  The  strain  tempera¬ 
ture  history  which  was  investigated  is  shown  in  Figure  39. 

Reoovery  data  were  fitted  to  Wool's  empirical  kinetic  theory: 

1  -Ro 

(1  ♦  Kt  )a 
62 


R  =  1  - 


(8) 


K  =  K0  exp  (-Ea/KT) 


(9) 


Experimental  parameters  are  compared  in  Table  8. 


One  interesting  observation  was  that  the  rate  parameter,  K,  in  equation  (9) 
was  constant  for  the  inert  propellant,  whereas  it  was  not  for  CSD  live  propellant 


The  experimental  results  are  shown  in  Figures  40  to  46.  It  can  be  seen 


Figure  39.  Strain  Temperature  Healing  History 
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TABLE  8.  HEALING  RESULTS 
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Based  upon  this  limited  study,  a  new  theory  of  healing  as  a  function 
of  stress  level  and  time  was  proposed  in  the  following  form: 


A  (o.t)  -f  t 

•*  r=- oo  . 


with  several  forms  possible  for  <P( t'  or  0(H): 


-  at 
e 

5  (t)  (11) 

Ks  U(t) 

K0  exp(-k-|(t-  r0)2) 

Analytical  solutions  are  obtainable  for  these  forms;  however,  the  "new"  constitu¬ 
tive  theory,  which  would  be  based  on  a  physical  molecular  model,  is  even  more 
complicated  at  the  present  time  than  available  phenomenological  continuum 
models. 


4.1.3  Farris  Nonlinear  Theory  for  Solid  Propellants 

The  WQrk  of  R.  Farris^  was  a  major  attempt  at  predicting  the  nonlinear 
response  of  solid  propellants  in  rocket  motor  analyses. 
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Elongation,  in. 


Log  lime,  min 

Healing  Master  Curve  for  Inert  Propellant 


20  C  lest 


First  load  • 


Second  load  ■ 


First  unload 


Second  unload 


Strain  (  e  )  — ■ 

Figure  43.  Load-Unload  Curves  for  PBAN  Propellant  with  Zero  Healing  Time 
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First  load  - 


■Second  load 


■First  and  second  unload 


Strain  (  e  ] 


Figure  44.  Load-Unload  Curves  for  PBAN  Propellant  with  1000-min  Healing  Time 
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Figure  45.  PBAN  Propellant  Recovery 


Experience  with  Farris’  stress-strain  law  at  CSD  showed  that  this  theory 
could  not  predict  the  response  of  solid  propellants  under  strain  histories 
that  were  not  included  in  the  set  used  for  material  characterization.  For 
loading  histories  of  the  types  included  in  such  a  set,  however,  the  predictions 
were  generally  acceptable. 

Farris'  nonlinear  theory,  in  the  form  presented  next,  was  considered  in 
the  first  phase  of  the  program  to  compare  it  with  the  other  five  approaches 
originally  proposed. 

Based  upon  previous  work  on  rubber  elasticity,  Farris  postulated  a  model  to 
account  for  the  permanent-memory  effects  exhibited  by  many  solid  propellants 
under  uniaxial  loading.  This  model  presumes  the  existence  of  inhomogeneities  in 
the  local  strain  field  between  filler  particles,  a  distribution  of  polymer 
chain  lengths  between  filler  particles,  and  a  uniform  failure  strain  for  each 
polymer  chain.  This  model  has  been  successful  in  predicting  the  nonlinear  perma¬ 
nent  memory  response  of  solid  propellants  before  dewetting  when  the  models  predict 
the  same  response  in  compression  as  in  tension.  This  prediction  is  in  agreement 
with  experimental  observations,  although  the  molecular  mechanisms  contributing 
to  the  permanent  memory  response  in  compression  are  clearly  different  from  those 
in  tension. 

Once  dewetting  occurs,  the  model  is  modified  to  account  for  vacuole  forma¬ 
tion  and  different  results  in  compression  and  tension  are  expected.  Farris 
presented  the  constitutive  equation  as  the  sum  of  an  essentially  time  independent 
bulk  stress  and  a  time-dependent  deviatoric  stress  o^d  ,  so  that  in  general 

d  (12 ) 

»ij(t)  =0B  5iJ  +aij  (t) 


where  6^  is  the  Kronecker  delta  equal  to  unity  if  i  =  J  and  zer  otherwise.  The 
form  developed  for  the  deviatoric  stress  is 
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d  d  /  I>  \  °2  d 

aiJ  (t)  =  e-BI  d/I^  Ai  ejj  (t)  A2  I -  J  ejj  (t) 

\  II  II  p2  / 


I* 

/.d  /  T7  \  m4  > 

A3  (t  -  |)  ejj  (?)  d?  ♦  [  \  /  *n  (t 

VM>J  J 


-  ?)  *ejj  (?)  d? 


where: 


Ij  =  volume  dilatation  =  *  ep2  ♦  633  for  small  strains 

1  r  2  2 

I-y  =  octahedral  shear  strain  =  —  j ( e -j -j  -  ep2)  ♦  (e22  -  @33) 

3  2-.1/2 

+  (e33  ~  e11)  J 

d 

ejj  a  deviatoric  strain  tensor 

B»  *i ,  A2,  A3,  Ajj,  M2,  M4,  P2,  P4  ,  a  constants 


r  1  p  i1/Pi 

Mpt =  /'  I7(f)  I  Xdf 


Equation  (13)  has  been  applied  to  reasonably  complex  deformation  histories  using 
unpressurized  and  pressurized  uniaxial  and  biaxial  test  specimens.  The 
agreement  was  not  as  good  as  would  have  been  desirable,  but  it  was  still  better 
than  with  linear  viscoelasticity.  Time-temperature  superposition  was  included 
in  equation  (13)  by  introducing  a  time-temperature  shift  factor,  A^,  and 
redefining  the  Lp  norm  of  equation  (14)  in  the  form: 


i  M. 


)l7(f)|pi  \i/Pi 

- df  j 

aT(f)  / 


with  |,|  denoting  the  absolute  value.  Experimental  data  for  simultaneous 
cooling  and  straining  have  been  fit  using  equation  (13)  with  the  introduction 
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of  a  time- temperature  shift  function  AT  through  equation  (15).  The  justification 
for  introducing  an  A?  in  the  above  manner,  however,  is  not  immediately  obvious 
nor  adequately  explained  in  the  available  literature. 


To  represent  the  response  to  interrupted  and  cyclic  constant  strain  rate 
tests,  equation  (13)  was  modified  by  setting  P4  =  oo  and  A3  =  -A4  so  that: 


The  multiplier  for  the  hereditary  integral  in  equation  (16)  vanishes  whenever 
the  current  value  of  1 7  is  at  its  largest  and  is  non-zero  for  all  other  values. 
This  representation  allows  for  viscoelastic  (fading  memory)  response  on 
unloading. 

The  bulk  stress,  Og  in  equation  (12),  was  taken  to  be  essentially  time- 
independent,  although  there  is  coupling  between  distortion  and  dilatation  as 
indicated  in  the  exponential  multiplier1  in  equations  (13)  and  (16). 


The  first  attempt  to  represent  the  bulk  stress  took  the  form  of  a  series: 


"  i  1 

i  Aij  ±(j  I7  ?  Aqo  =  0 


This  equation  adequately  predicts  the  bulk  stress  as  long  as  it  does  not 
vary  greatly.  When  a  hydrostatic  pressure  is  superimposed,  however,  very 
poor  results  are  obtained.  In  an  attempt  to  overcome  this  difficulty,  Farris 
modeled  the  compressibility  of  the  gas  voids  caused  by  vacuole  dilatation 
by  treating  them  as  spherical  voids  contained  in  an  elastic  medium.  Assuming 
that  the  voids  themselves  offer  no  resistance  so  that  all  void  dilatation 
is  caused  by  distortion  of  the  surrounding  elastic  material,  the  void  content 
at  zero  pressure  may  be  represented  as  a  power  law  in  terms  of  the  octahedral 
shear  strain  ly,  and  assuming  that  the  bulk  behavior  varies  linearly  with 
hydrostatic  pressure,  P,  yields: 


Id  =  C-i  P 


l2 


08) 


for  the  dilatation,  where  C^,  C2,  and  n  are  constants,  and  G  is  the  shear  modulus 
of  the  elastic  matrix  material. 


The  Farris  approach  to  nonlinear  characterization  was  only  superficially 
investigated  at  the  start  of  this  program,  since  our  previous  experience  on  the 
Failure  Mechanism  Program  had  demo..jtrated  that  this  approach  was  not  useful. 
However,  it  was  felt  that  the  Farris  code  could  give  a  meaningful  baseline  for 
comparison  to  the  Quinlan  theory,  at  least  the  predictions  made  in  AFRPL-TR-78- 
37.  Consequently,  the  same  12  cycle  test  data  that  were  reported  in  AFRPL-TR- 
78-37  were  input  into  the  Farris  nonlinear  code.  Results  of  this  computation 
are  shown  for  the  first  two  cycles  in  Figures  47  and  48.  Note  that  during  the 
first  cycle  (Figure  47),  the  Farris  fit  gives  much  better  correlation  on  both 
the  loading  and  unloading  parts  of  the  cycle.  The  second  cycle  shows  less  clear 
agreement  on  the  loading  part  of  the  cycle,  similar  agreement  at  the  peak 
stress,  and  better  agreement  with  the  measured  stress  on  the  unloading  cycle. 
Figure  49  shows  the  error  in  peak  stress  calculation  for  each  cycle  for  both  the 
Quinlan  fit  and  the  Farris  fit.  This  measure  of  error  also  shows  that  the 
Farris  theory  fit  the  data  quite  well.  The  eleventh  cycle  data  is  questionable. 
This  is  the  highest  strain  rate  test  and  in  fact  could  be  erroneous  as  the  test 
rate  is  specified  as  4.0  in. /min,  which  is  not  common  nor  even  possible  on  some 
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Figure  47.  Quinlan  N.C.  Calculation  with  TPH  1011  Propellant  (First  Cycle) 


Stress  (o ),  psi 


vs 


TABLE  9.  FARRIS  COMPUTED  N-L  PARAMETERS 
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Constants 

Farris  ASPC  Final  Report 
TPH  1011  (54  Tests) 

Quinlan  Test  j 

Run  21 

Run  22 

Run  23 

Run  24 

B  .  1 

60 

60* 

94.0 

60.0* 

93.8 

nNT  +  ' 
°NC 

-0.150 

-0.150* 

-0.150* 

-0.45 

-0.237 

BNC  +  1 

-1.69 

3.94 

-3.32 

8.15 

-3.92 

®NC  ♦  2 

-1.23 

-1.23 

-1.23 

-1.23 

-1.23 

BNC  ♦  3 

10.5 

12.5 

8.1 

18.4 

10.8 

BNC  4 

2.00 

19.9 

19.6 

19.1 

18. 1 

BNC  ♦  5 

102.0 

-44.0 

15.4 

129.0 

128.0 

bNC  6 

239.0 

151 

426 

48.0 

352.0 

bNC  7 

56.5 

167 

90.9 

17.8 

32.0 

Average  deviation 

mmm 

-14.2 

■BBB 

-9.0 

Standard  deviation 

EEH 

24.1 

■Sal 

26.3 

•  Denotes  held  fixed  j 

Instron  machines,  whereas  all  other  cycles  were  at  the  usual  intervals  of  1,  2, 
and  5  in. /min.  Table  9  shows  the  valu< s  for  the  Farris  parameters  from  several 
solutions  to  the  Quinlan  data.  Shown  for  comparison  in  the  first  column  are 
data  generated  by  Aerojet  during  the  development  of  the  Farris  constitutive 
theory  on  the  same  propellant  (circa  1972),  TPH  1011.  Note  the  range  in  values 
from  the  various  runs.  This  is  reminiscent  of  the  CSD  experiences  with  the  fail 
ure  mechanisms  program. 

It  might  also  be  noted  in  passing  that  each  of  these  solutions  cost  about 
$1.50  in  computer  time  as  compared  to  about  $100.00  for  the  CSD  program  to 
evaluate  the  Quinlan  theory.  Hence,  computation  costs  may  be  highly  dependent 
on  the  theory  chosen  and  the  sophistication  of  the  techniques  used  to  evaluate 
the  constitutive  expressions. 

One  additional  test  was  run  at  the  request  of  Dr.  J.  E.  Fitzgerald,  a  con¬ 
sultant  to  CSD  on  this  program.  He  requested  a  prediction  for  a  test  in  which 
a  24-hr  rest  period  was  included  in  the  test  history.  The  results  are  shown  in 
Figure  50  for  the  parameters  selected  from  Table  10.  Contrary  to  what  was 
expected,  the  predicted  response  did  not  show  anomalous  behavior  on  the 
reloading  following  the  rest  period. 


96 


Fitzgerald  2M-hr  "Dip"  Test 


TABLE  10.  PARAMETERS  FOR  FARRIS  FIT 


T9991 


Contstants 

QEX12 

GEX3A 

GEX3B 

0NT 

18 

18 

18 

^C 

-0.49 

-0.27 

-0.61 

bNC  +  3 

5.44 

2.90 

3.00 

bNC  ♦  4 

2.20 

13.0 

19.0 

bNC  +  5 

152 

30 

140 

bNC  ♦  6 

102 

360 

225 

bNC  ♦  7 

15 

22 

11 

4.2  LINEAR  VISCOELASTIC  CONSTITUTIVE  EQUATION 
4.2.1  Isothermal  Tests 

The  one-dimensional  stress-strain  law  for  a  theroorheologically  simple 
linear  viscoelastic  solid  may  be  expressed  as: 

i 

•<«'/  E  <St-S,)$(r)  dr  (, 

0 

where: 

os  stress 
e*  strain 

E(t)  s  relaxation  modulus  (FRONT  series  representation 
using  a  matrix  solution  for  curve  fitting  data; 

GSD  Data  Analysis  Procedure  No.  7.3) 

$£  -Sr  ■  temperature-shifted  time,  given  by: 


(20) 


(21) 


in  which  TR  is  the  shift  reference  temperature,  and  both  Ta  and  m  are  material 
parameters.  The  material  parameters  were  obtained  using  CSD  Data  Analysis  Pro¬ 
cedure  No.  7.4,  which  is  curve  fit  routine  using  Powell's  algorithms. 

The  linear  viscoelastic  model  was  used  to  predict  the  response  of  UTP- 
19,360B  and  UTP-3001  under  several  strain  histories.  The  corresponding  results 
are  included  here  as  a  basis  against  which  to  compare  the  stress  predictions 
obtained  using  the  nonlinear  stress-strain  laws  considered  in  the  program. 

4.2.2  Isothermal  Stress  Predictions 

The  measured  response  is  compared  against  that  predicted  by  linear 
viscoelasticity  for  UTP-19.360B  in  the  following  order  (Figures  51  through  76). 

The  results  for  the  lowest  and  highest  constant-rate  tests  (Test  No.  1) 
appear  in  Figures  51  and  52.  Those  for  the  dual-rate  tests  (high-to-low  and 
low-to-high,  Test  No.  3)  are  shown  in  Figures  53  and  54.  Figure  55  contains  the 
comparisons  for  a  saw-tooth  strain  history  (Test  No.  5)  with  increasing  strain 
peaks  and  rest  periods  between  cycles.  The  results  corresponding  to  complex  mul 
tiple  loading,  24-hour  relaxation,  long-duration  similitude,  three- step 
relaxation,  and  predamage  relaxation  are  presented  in  Figures  56  to  60.  In  addi 
tion,  the  time- temperature  superposition  principle  is  put  to  use  with  constant 
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(Text  continued  on  pg.  126) 


Figure  54.  Linear  Viscoelastic  Stress  Predictions  for  UTP-19,360B-400/1777  Two-Rate  Test  History 

(Code  No.  3) 


UOO/ 1777  Multiple  Loading  Test  History 


TEKRERRTURE  HISTORTi  7?  F 
3TRRIH  HISTORTi  }B  X 


Linear  Viscoelastic  Stress  Predictions  for  UTP-19, 360B-400/1777  24-hr  Relaxation 


/m/ 


UHPEAB 
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Figure  60.  Linear  Viscoelastic  Stress  Predictions  for  UTP-19 , BSOB-^OO/I??? 

Predamap;ed  Relaxation  (Test  No.  9) 


Figure  61.  Linear  Viscoelastic  Stress  Predictions  for  UTP-19 , 360B-400/1777  Constant-Rate  History 


Linear  Viscoelastic  Stress  Predictions  for  UTP-19 , 360B-H00/1777  Constant-Rate  Test  History 


Figure  64.  Linear  Viscoelastic  Stress  Predictions  Tor  UTP- 3001 -750/7768  Constant -Bate 

(Code  No.  1) 


/7768 


in  * 
co 


Figure  66.  Linear  Viscoelastic  Stress  Predictions  for  UTP-3001 -750/7768  (Sample  No,  2) 

Multiple  Rate:  High  to  Low  (Test  No.  3) 


3001-750/7768 


MEASURED  STRAIN  .  .  TEMPERATURE  HJST0RT:  74  F 

STRAIN  HISTORY:  12.8.12.4.8.0.4.0 


Complex  Multiple  Loading  (Test  No.  10) 


Figure  70.  Linear  Viscoelastic  Stress  Predictions  for  UTP-3001-750/7768 

2’’ -hr  Relaxation  (Test  No.  8) 


Linear  Viscoelastic  Stress  Predictions  for  UTP- 3001-750/7768 
Similitude  Test  History  (Code  No.  12) 


HERSUREO  STRAIN 


Figure  72.  Linear  Viscoelastic  Stress  Predictions  for  UTP-3001-750/7768 
Short  and  Long  Time  Relaxation  (Test  No.  13) 


re  73.  Linear  Viscoelastic  Stress  Predictions  for  UTP-3001-750-7768 
Predama^e  Relaxation  (Test  No.  9) 


Figure  7^.  Linear  Viscoelastic  Stress  Predictions  for  (JTP-300 1-750/7768 
Constant-Rate  Test  History  (Code  No.  1) 


Linear  Viscoelastic  Stress  Predictions  for  UTP-3001-750/7768 
Constant-Rate  Test  History  (Code  No.  1) 


rate  testa  (Test  No.  1)  at  70,  40,  and  123  F,  as  shown  in  Figures  61  to  65, 
respectively.  Stress  predictions  for  the  same  tests  on  UTP-3001  are  shown  in 
Figures  66  through  76. 

Predictions  for  biaxial  samples  were  also  made  using  equation  (19).  For 
OTP-19, 360B,  it  was  assumed  that  Eblaxial  =  4/3  EUniaxial*  This  assumption  is 
supported  by  the  limited  biaxial  modulus  tests  conducted  on  the  program.  Predic¬ 
tions  for  UTP-3001  were  made  by  substituting  the  biaxial  relaxation  modulus  for 
that  of  the  uniaxial. 

The  results  for  UTP-19,360B  are  given  as  follow.  The  lowest  and  highest 
constant  rate  tests  (Test  No.  14)  are  shown  in  Figures  77  and  78.  Figure  79 
contains  the  comparisons  for  the  stress  relaxation  history  (Test  No.  16).  The 
results  for  UTP-3001  for  the  same  tests  are  given  in  Figures  80  to  82. 

4.2.3  Linear  Thermoviscoelasticity  Analyses  of  UTP-19.360B  Propellant 

Simultaneous  cooling  and  straining  tests  of  uniaxial  and  biaxial  specimens 
and  the  complex  biaxial  ramp-relaxatic..  ramp-test  while  cooling  and  heating  have 
been  analyzed  according  to  linear  viscoelasticity  theory  in  order  to  provide  a 
baseline  for  comparison  with  the  nonlinear  theoretical  predictions  and  to 
evaluate  the  error  between  linear  predictions  and  measured  experimental  results. 

The  one-dimensional  equations  are  as  shown  in  equation  (19),  with  the 
exception  that  e  is  replaced  by  e  -  oAT  and  AT  =  T0-T(t).  In  order  to  simplify 
numerical  calculations,  the  time- temperature  shift  function  was  represented  by 
the  power  law  given  in  equation  (21)  . 

Calculations  have  been  made  for  prony  series  and  power- law  representations 
of  the  relaxation  modulus.  Both  types  of  computations  give  comparable  results 
and  the  computer  codes  have  been  verified  with  handbook  type  calculations. 


Figure  78.  Linear  Viscoelastic  Stress  Predictions  for  UTP-19, S^OB-^OO/I??? 

Biaxial  Constant-Rate  (Test  No.  1*0 


Figure  79.  Linear  Viscoelastic  Stress  Predictions  for  UTP-19,360B-1400/1777 
Biaxial  Stress  Relaxation  (Test  No.  16) 


Using  Prony  series  representation  for  the  modulus  and  assuming  a  con¬ 
stant  coefficient  of  linear  thermal  expansion,  equation  (19)  car  be  written 
in  the  form 


oft) 


f  e  -  a  AT) 


where  it  is  also  assumed  that  e  and  AT  are  constant.  Substituting  equation  (21) 
into  equation  (20)  and  performing  the  indicated  integrations  leads  to 


oft) 


(e  -a AT) 


-au1+m 

e 


(22) 


where: 


a 


1 

fl  +  m) 


(23) 


Equation  (22)  was  evaluated  numerically  calculating  the  constant  e  and  T  at  each 
data  point. 

Assuming  a  power  law  representation  for  the  modulus  of  the  form 

Eft)  =  Ee  +  E]  t'"  (24) 


leads  to 


oft)  =  (  c  -  oAT) 


Ee  + 


E-!  f  t/AT)-n 


-It 


(25) 
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where 


1-n  /*1  p-  +  n-2\ 

/  (l-u)-m-uWl  'du  (26) 

J(H  +  ATn> 

Figures  83  through  85  compare  the  measured  and  predicted  stresses  for  the 
uniaxial  tests  and  Figure  86  compares  stresses  for  the  biaxial  test.  The  good 
agreement  with  the  biaxial  straining  and  cooling  test  is  believed  to  be  fortui¬ 
tous.  In  fact,  the  test  data  is  questionable,  since  the  linear  predictions  did 
not  agree  even  during  the  initial  cooldown  of  the  biaxial  ramp-relaxation  ramp 
test  as  shown  in  Figure  87. 

Handbook  type  calculations  were  also  made  to  determine  a  thermomechanical 
coupling  coefficient,  Ap,  with  the  results  shown  in  Figures  88  through  90.  The 
results  obtained  using  a  constant  Ap  agree  reasonably  well  with  experiment, 
with  the  exception  of  the  biaxial  tests.  Using  an  Ap  calculated  from  uniaxial 
tests,  the  predictions  shown  in  Figures  91  and  92  were  made.  Observe  that  the 
initial  cooldown  of  the  biaxial  ramp-relaxation-ramp  test  is  matched  well,  but 
the  cooling-only  test  is  not,  giving  further  credence  to  our  contention  that 
these  test  results  are  invalid. 

4.3  M.  QUINLAN'S  THEORY  OF  MATERIALS  WITH  VARIABLE  BONDING 
4.3.1  Original  Model 

In  developing  a  mathematical  framework  for  his  stress-strain  law,  Qunilan,1* 
reasoned  that  since  propellants  consist  of  minute  rigid  particles  embedded 
in  a  polymer  matrix,  such  materials  would  respond  to  a  deformation  process 
with  a  change  in  the  amount  of  species  to  species  bonding.  He  thus  proposed 
to  correct  the  deficiencies  of  fading-memory  type  theories  by  introducing  a 
correction  term  that  accounted  for  the  changes  in  the  state  of  bonding  that 

(Test  continued  on  pg.  145) 
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Linear  Viscoelastic  Stress  Predictions  for  UTP-19, 36013-400/1777 
6-in.  Bar  Straining  -  Cooling  (Test  No.  18) 


Linear  Viscoelastic  Stress  Predictions  for  UTP-19 , 3(>0B-^00/M71 
Biaxial  Ramp-Relxation-Ramp  (Test  No.  21) 
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are  induced  by  a  deformation  process.  His  constitutive  model  then  took  the 
f  orm: 

o  =of  +  ob  (27) 

in  which: 

o  =  current  stress 
=  fading-memory  type  stress 

=  stress  correction  due  to  change  in  the  state  of  bonding. 

Motivated  to  some  extent  by  reaction-rate  theory,  Quinlan  modeled  the  evolution 
of  the  bonding  state  through  the  following  ordinary  differential  equation: 


7r=a  -  qv  ( <t>- 

subject  to  the  initial  condition: 

ir(0)  =  1 


(28) 


(29) 


in  which  x represents  the  state  of  bonding;  a  ,  m  ,  and  v  are  material  parameters, 
and 

<p  =  1  (30) 


is  the  stretch  ratio;  with  e,  the  strain. 

The  unique  solution  of  equation  (27)  may  be  readily  obtained  for  piecewise 

ji 

linear  stretch  histories. 

Taking  a  linear  viscoelastic  relation  for  and  considering  the  stress 
correction  term,  ob,  as  proportional  to  the  state  of  bonding,  Quinlan  arrived  at 
the  following  stress-strain  law: 
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a  (t) 


■/ 


G(t  -  t)  4>  (r)  dr  +  B0  rr(t) 


with: 


G(t)  =  E0 


n  (t)  =a  1  M 


*  -  n  [l  -  e* (*-»)]  1 


jr  (t  =  0)  =  1 


The  six  parameters  (EQ,  n,  BQf  a , n  and  v  )  needed  in  this  theory  to  characterize 
a  propellant  may  be  obtained  by  fitting  the  model  to  the  observed  response  of 
the  material  when  subjected  to  a  saw-tooth  strain  history  that  has  increasing 
peak  strains  and  sufficiently  long  rest  periods  between  cycles.  Alternatively, 
the  studies  reported  in  the  literature  on  the  effects  of  employing  different 
data  bases  for  characterization  show  that  the  test  history  should  primarily 
include  the  maximum  expected  strain  level,  the  expected  range  of  strain  rates, 
as  well  as  rest  and  relaxation  periods. 

The  model  was  used  to  predict  the  response  of  TP— HI  Oil  under  several 
loading  histories,  and  it  reproduced,  somewhat  accurately,  the  general  trend  of 
observed  solid  propellant  behavior.  CSD  developed  a  computer  code  to  automate 
Quinlan’s  constitutive  equations. 


HJ6 


The  parameters  entering  Qunilan's  law  were  determined  by  the  code  through 
minimization  of  the  root-mean-square  (rms)  error  between  the  observed  stresses 
and  those  predicted  via  the  foregoing  equations.  This  procedure  may  be  formally 
stated  as  follows: 

Determine  E,  n,  B,  a,  v,  v  so  that 


acquire  a  minimum  value,  subject  to 
E,  B,o,P,r>0;  and  -1<n<0 

where  N,  ojL,  and  oci  stand,  respectively,  for  the  number  of  data  points,  the  ith 
measured  stress,  and  the  corresponding  stress  computed  from  the  proposed  law. 

Such  an  optimization  process  is  performed  in  the  code  with  the  Powell-BOTM 
sequential-search  algorithm. 

4.3.1 . 1  Verification  of  Published  Experiment 

In  reproducing  Dr.  Quinlan's  experiment,  CSD  took  as  a  basic  data  set  the 
reported  values  of  time,  measured  stress,  and  stretch  rate.  The  corresponding 
saw-tooth  stretch  history  was  internally  built  by  the  computer  code.  This 
approach  introduced  a  difference  in  the  degree  of  accuracy  of  representation 
between  the  values  contained  in  Dr.  Quinlan’s  report  and  those  used  by  CSD 
(Table  11).  Thus,  small  discrepancies  were  observed  in  verifying  Quinlan's 
work.  An  overall  measure  of  such  deviations  is  given  by  the  rms  error  per 
cycle.  While  Dr.  Quinlan's  reported  value  is  2.7  psi,  that  obtained  at  CSD 
was  2.98  psi. 

4. 3. 1.2  Verification  of  Uniqueness  of  Solution 

Establishing  the  uniqueness  of  parameter  values  for  use  with  the  nonlinear 
law  is  of  capital  importance.  To  this  end,  test  history  number  4  of  Table  12  was 


TABLE  11.  SUMMARY  OF  QUINLAN'S  AND  CSD'S  DATA 
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12-rain  rest  at  zero  strain  between  each  cycle 


TEST  HISTORIES 


employed  as  the  basic  input  data  for  CSD's  optimization  and  stress  prediction 
code,  with  two  distinct  sets  of  initial  guesses  for  the  parameters.  As  shown  in 
Table  13,  the  end  results  of  both  computer  runs  are  essentially  the  same. 


4. 3. 1.3  Modifications  to  Damage  Equation 

Under  the  hypothesis  of  a  piecewise  linear  stretch  history,  say, 

0(t)  =  <t>{ tQ)  +  R(t-tQ);  tQ<  t  <  tf 


the  damage-evolution  equation  provides  the  solution: 


evU-.)=c  j,  .  e-fU-to)  [,  .  (V]}'1;  (37) 

-1 

;  whenever  c  =  0 

for  tQ^t<tf,  where  the  following  auxiliary  parameters  have  been  used: 


T  (38) 

f(t0)5  0(to)  -7T[  <t0)] 

Now,  continuous  piecewise  linear  stretch  histories  have  discontinuous 

stretch-rate  histories,  as  schematically  shown  in  Figure  93«  It  does  make  a 

difference  whether  the  right-  or  the  left-hand  derivative  stretch  is  employed  in 

computing  the  evolution  of  damage.  If  the  right-hand  derivative  is  used,  a 

Jump  discontinuity  in  the  response  is  obtained  which  propagates  through  the 
»> 

cycle  because  0  also  affects  the  value  of  ev(  <P-n)  through  the  auxiliary  parameter 
If,  on  the  other  hand,  the  left-hand  derivative  is  consistently  used,  no 
jump  phenomenon  is  observed.  Additionally,  since  in  solving  the  damage-evolution 
equation,  the  values  of  0  and  ir  at  the  end  of  the  previous  ramp  constitute 


or: 


’(0  -  »)  =  a 


A<  v 


TABLE  13-  UNI QUENESS - OF -S OLUT ION  TEST 


the  initial  condition  for  the  present  ramp,  another  interpretation  is  yet 
possible  of  the  0  value  to  be  used,  namely,  to  compute  £(tQ)  with  the  left- 
hand  derivative  of  the  stretch  (i.e.,  this  corresponds  to  the  final  value 
of  the  previous  ramp),  and  evaluate  ir  (tQ)  using  the  right-hand  derivative 
of  the  stretch.  This  alternative  yields  a  jump  only  at  the  start  of  a  cycle 
and  nowhere  else. 

From  the  foregoing  discussion  it  is  apparent  that  the  damage  evolution 
equation  may  be  unambiguously  defined  as: 

!(♦,)  -  .{i'-*  -N  [1-  ev(‘-)]}  (39) 

where  0^  represents  the  left-hand  derivative. 

CSD  proposed  a  modified  version  of  tue  damage-evolution  equation  in  the  form 

*(0t)  -  a  -  1)  -  v  [l  -  eV^  |  (40) 

which  would  allegedly  recover  the  equilibrium  component  of  the  response  and  at 
the  same  time  would  eliminate  the  jump  phenomenon.  The  reduced  equation  (in 
terms  of  {  =  0  -  w  )  corresponding  to  this  version  constitutes  a  differential 
equation  with  variable  coefficients  which  may  be  numerically  solved.  The 
required  subroutines  to  accomplish  this  have  already  been  developed  at  CSD  and 
the  performance  of  the  modified  damage  equation  (either  explicitly,  as  in  equa 
tlon  (39),  or  implicitly  through  0  ,  as  in  equation  (40))  was  tested  modifying 
the  damage  equation  according  to: 

;  -  o*  -  1] 


(41) 


m 


Use  of  equation  (41)  seriously  affected  time  to  convergence;  furthermore, 
stress  predictions  based  on  this  equation  yielded  larger  errors  than  those 
based  on  equation  (39). 

Finally,  it  is  interesting  to  note  that  even  when  an  inconsistent  fora  of 
equation  (40)  was  used,  whereby  i  was  computed  from  equation  (40)  but  *  itself 
wets  determined  from  equations  (37)  and  (38),  better  stress  predictions  were 
obtained  than  with  Quinlan's  original  form.  In  this  case,  the  first  two 
cycles  of  Quinlan's  experiment  were  fitted  with  these  two  approaches  starting 
with  the  same  initial  values  for  the  parameters.  After  the  second  iteration 
of  the  optimization  routine,  the  error  corresponding  to  the  inconsistent 
form  was  always  smaller  than  that  for  Quinlan's  original  equation. 

The  previous  comments  seem  to  Indicate  that  equation  (40)  constitutes 
an  appropriate  definition  of  the  damage  evolution.  It  must  be  tested,  however, 
under  different  types  of  stretch  histories  before  it  can  be  either  accepted 
or  discarded. 

4. 3. 1.4  Effect  of  Data  Base 

Employing  the  stretch  histories  listed  in  Table  12,  a  preliminary 
evaluation  was  carried  out  of  the  features  a  test  should  include  to  properly 
characterize  the  material  response.  The  tests  described  below  were  performed  to 
gain  some  insight  into  the  effects  of  data  base  on  stress  predictions  and  are  by 
no  means  exhaustive  for  the  purpose  of  completely  ascertaining  such  effects. 

4. 3. 1.5  Reducing  the  Data  Base 

History  No.  1  (H.1)  of  Table  12  was  curve-fitted  and  the  corresponding 
parameters  were  used  to  predict  material  response  under  H.2.  Errors  in  the 
predicted  response  seem  to  depend  on  whether  or  not  the  stretch  maximum  and/or 
the  stretch  rates  are  included  in  the  data  base  used  to  determine  the  constitu¬ 
tive  parameters.  Typical  plots  of  these  results  appear  in  Figures  94  to  97. 

Based  on  the  above  observations,  H.4  was  next  used  to  determine  the 
parameters,  and  the  response  to  H.5  was  then  predicted  with  virtually  the  same 
ras  error  per  cycle.  It  is  pointed  out  here  that  the  base  history  (H.4) 


153 


Figure  94.  Effect  of  Jump  and  Strain  Maximum 
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included  all  the  stretch  rates  of  the  history  employed  for  stress  predictions 
(H.5),  but  did  not  contain  the  stretch  maximum.  The  last  fact  introduced  sig¬ 
nificant  errors  in  the  response  prediction  along  the  ramp  leading  to  the  maximum 
stretch,  as  shown  in  Figure  98. 
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Strain  (  e  ).  % 

Figure  95.  Effect  of  Jump  and  Strain  Maximum 
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It  is  interesting  to  note  that  the  measured  response  during  the  third 
stretch  plateau  of  H.5  corresponds  to  a  damaged  recovery  relaxation;  and  even 
though  the  base  history  ends  before  this  phenomenon  takes  place,  the  predicted 
response  to  H.5  actually  reproduces  this  type  of  behavior,  as  depicted  in 
Figure  99. 


4. 3. 1.6  Extrapolating  to  a  Different  Test 

History  No.  2  (Table  12),  with  Quinlan's  reported  values  for  the 
parameters,  was  used  to  predict  response  to  H.3  without  success.  Material  aging 
and/or  differences  in  stretch  maxima  as  well  as  in  stretch  rates  might  have  been 
responsible  for  these  bad  predictions. 


»vv 


Strain  (  e  ) .  % 

t 

Figure  96.  Effect  of  Strain  Maximum  and  Strain  Rate  (0  =  0.0029) 
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Next,  the  parameters  corresponding  to  H.H  were  used  to  predict  response 
under  H.3;  the  samples  used  in  both  tests  came  from  the  same  propellant  batch. 
The  rms  error  per  vcle  of  the  predicted  response  to  H.3  with  parameters 
corresponding  to  H.M  was  approximately  13-75  psi,  whereas  that  obtained  with 
parameters  corresponding  to  H.3  was  about  13.2. 


Alternatively,  when  parameters  obtained  by  curve-fitting  H.3  were  used  to 
predict  material  response  under  H.5,  large  errors  resulted.  These  errors  may  be 


attributed  to  differences  in  stretch  maxima  and  stretch  rates  present  in  the  two 

tests  - 

Fran  these  preliminary  observations,  it  appears  that  parameters  adequate 
for  stress  prediction  purposes  should  be  determined  from  tests  which  cover  the 
expected  range  of  stretches  and  stretch  rates. 
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Note 

(1)  Ninth  cycle  Loading  Ramp 

(2)  4>  =  0.01 15  in. /in. /sec; 

=  6  5% 


Symbol 


Approach 


Measured 

Predicted  from  first  two 
cycles  only  (  i  =  0.0058 


Strain  (  e  ) ,  % 

Figure  97.  Effect  of  Strain  Maximum  and  Strain  Rate  ($  =  0.0115) 
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4. 3. 1.7  Temperature  Effects 

As  a  first  step  to  ascertaining  the  effect  of  temperature  on  the  constitu¬ 
tive  parameters  entering  Quinlan’s  constitutive  law,  isothermal  tests  were 
performed  at  three  different  temperatures.  A  saw-tooth  strain  history  with 
successively  higher  peaks  (2,  4,  and  6$  nominal  strain)  and  without  rest  periods 
between  cycles  was  used  in  the  tests.  The  corresponding  computer  determined 
parameter  values  are  summarized  in  Table  14. 

These  preliminary  results  seem  to  indicate  monotonically  decreasing 
behavior  of  the  parameters  with  increasing  temperature  except  for  the  coeffi¬ 
cients  a <and  u,  which  appear  to  have,  respectively,  a  bell-shaped  and  a  monoto¬ 
nically  increasing  type  of  dependence  on  temperature. 


Additional  tests  at  different  temperatures  are  needed  to  more  precisely 
define  the  functional  dependences  previously  referred  to  and  to  hopefully  remove 
part  of  the  hidden  nature  of  the  parameters  involved. 


3  4  5  6 

Strain  ( e  ),  % 

Figure  98.  Effect  of  Strain  Maximum  (e 
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Reproduction  of  Modulus  Recovery 


TABLE  14.  EFFECT  OF  TEMPERATURE  ON  CONSTITUTIVE  PARAMETERS 


T999I4 


Parameters 


E 

n 

B 

a 


Temperature,  F 

2 

70 

120 

9019 

1543 

965 

-0.327 

-0.141 

-0.142 

547,250 

40,210 

15,458 

0.00236 

0.0675 

0.00427 

0.229  x  10‘S 

0.0735 

1.415 

321 

6.73 

4.91 

rms  error/ 
cycle,  psi 


4. 3. 1.8  Desirable  Test  Features 

Test  data  used  to  determine  the  constitutive  parameters  have  a  definite 
bearing  on  response  predictions  and  may  seriously  limit  the  applicability  of  the 

corresponding  results. 

From  a  practical  standpoint,  the  preliminary  studies  on  data  base  effects 
for  stress  prediction  purposes  indicate  that  it  is  not  imperative  to  duplicate 
the  expected  loading  history;  rather,  the  test  history  should  primarily  include 
the  maximum  expected  strain  level,  the  expected  range  of  strain  rates,  and,  if 
applicable,  other  characteristics  such  as  rest,  or  relaxation  periods.  On  the 
basis  of  these  observations,  it  appears  that  a  test  may  be  so  designed  to 
include,  within  a  reasonable  time  span,  the  minimum  amount  of  information 
necessary  to  properly  characterize  the  propellant  for  the  intended  use.  Two 
possible  stretch  histories  are  shown  in  Figure  100. 

4.3.2  Current  Model 

The  theory  developed  by  Quinlan  has  undergone  several  changes,  mainly  in 
the  expression  defining  the  evolution  of  damage.  With  the  theory  as  developed, 
Quinlan  has  not  been  able  to  satisfactorily  fit  stress  relaxation  data.  One 
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possible  modification  here  is  to  allow  the  relaxation  modulus  to  depend  cr 
damage  parameter,  ir .  This  modification  is  reasonable,  as  tests  generally  sh'w 
a  lowering  of  modulus  in  stress  relaxation.  In  particular,  Figure  101  showr 
four  stress-relaxation  tests,  each  corresponding  to  a  different  prestrair.  r 
the  three  upper  curves,  the  initial  strain  in  the  test  is  greater  than  or  e 
to  the  prestrain;  the  three  curves  approximately  coalesce.  For  the  lower  curvf 
however,  the  3%  strain  is  less  than  the  5%  prestrain  and  there  is  consider  v’: 
modulus  reduction.  This  type  of  strain- sensitive  damage  phenomenon  has  h 
reported  throughout  the  solid  propellant  literature.1®’  12 ’ 

The  most  effort  however,  has  been  spent  on  modifying  the  bonding  law  v. 
reduce  the  jump  phenomena,  to  improve  the  curve  fit  to  data,  and  to  ir.cl  ■ 
healing.  The  general  form 


o(t)  =  j  H(t-r)  <j>  (t)  dr  +  Cff(t) 
o 


(4. 


has  been  retained  in  all  versions  of  Quinlan's  theory,  but  the  rate  mechar 
underlying  damage 


ff=  p(ir,  0,  0) 


(43) 


for  which 


tt(0  =  1)  =  1 


was  assumed  to  contain  a  neutral  rate,  ?  ,  at  which  damage  remain  o- 


P(jt,  0,  {■)  =  0 
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Strain  (£  ).  V 


Figure  101.  Typical  Stress  Relaxation  Behavior  for  Undamaged 

and  Predamaged  Samples 


21967 


This  concept  allows  the  concept  that  at  rates  higher  than  f  ,  bond  breakage  would 
take  place,  while  at  rates  slower  than  the  neutral  rate,  bond  formation  would 
occur . 


Equation  (45)  for  the  neutral  rate,  ?  ,  may  be  rewritten  as: 


f  =  Q  (7T,  0  ) 

$  =  0  if  and  only  if  ir  =  0 


(46) 


which,  upon  expansion  in  Taylor  series,  becomes: 


(=Q(7r,7r  -  0  ) 


+  9Q  ( 7T  ,  v  -  0  ) 
0=0  a  (n  -  0  ) 


U  -  0)  +  Q  (ltr-0  2) 


0=  0 


The  first  term  on  the  right-hand  side  of  this  equation  vanishes  by  virtue  of 
equation  (46),  so  that,  neglecting  the  higher  order  terms,  and  defining: 


def  3Q  (»  ,  »  -  0) 

M  =  - 

d(n  -  0)  0=0 


leads  to  the  following  first-order  expression  for  the  neutral  rate: 


f  =  M(?r  -  0) 


where  for  bond  breakage 


0  -?>0 


while  for  bond  formation 


0  -  r<o 


In  addition,  equation  (43)  may  be  cast  in  the  following  form: 


jr  =  P  (jt,  0  ,  0 )  =  R(s,  u) 


where : 


s  =  s  (»,  0) 


U  =  0  —  f 


with  f  given  by  equation  (48) 
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Now,  under  the  assumption  that  the  process  of  bond  formation  is  slower  than 
that  of  bond  breakage,  the  function  R  defining  the  evolution  of  damage  must  be 
of  the  form  shown  in  the  sketch  below,  in  which  the  parameter  a  would  represent 
the  maximum  rate  of  bond  formation. 


u 


Hence,  R  may  be  defined  through  the  following  differential  equation: 

dR(s,  u)  r  .. 

=  v  I  R(s,  u)  -  a  (54) 


where  a  and  v  are  positive  constants. 


Integration  of  equation  (54)  yields: 


R(  s,  u)  =  a(ei;u  -  1 ) 


(55) 


Finally,  putting  equations  (48),  (53),  and  (55)  into  (51)  results  in: 


7i  =  a  |  exp  {<t>  -  *)]  -  1  | 


(56) 


Equations  (42)  and  (56)  subject  to  (44)  were  used  by  Quinlan  in  several 
ways  to  characterize  the  response  of  UTP-19,360B.  One  such  stress-strain  law 
took  the  following  form: 


in  which  Eo,  ,  E2,  n,  and  C  are  constants  and  €  is  the  strain.  Although  some 
aspects  of  propellant  behavior  were  better  modeled  than  with  the  original  version 
of  the  theory,  others  were  not,  and  further  revisions  were  necessary.  In  the 
latest  version  of  his  constitutive  law,  Quinlan  used  a  Prony  series  to  represent 
the  relaxation  function  and  changed  strain  for  stretch  in  the  original  equation 
of  evolution  for  damage,  so  that,  in  summary,  the  current  model  looks  as 
follows : 


0  (t)  =  ov(t)  + 

t 

ov(t)  =  Ge  e(t)  J  G(t  -  r)  e(r)  dr 

(58) 

n 

V  ”Tifc 
o(t)  =  2  ci e 

i  =  1 

ob(t)  =  Be7(e  “**  e(t) 

•  V(e  -  it) 

ir  (t)  =  a  e  '  ;  -  1  ;  tt  (0)  =  0 

This  final  version  of  the  theory  has  been  employed  to  characterize  UTP-19,360B, 
but  has  not  been  used  by  Quinlan  to  predict  the  response  of  the  propellant 
under  any  loading  history  other  than  a  single  characterization  test.  Its 
generalization  to  transient  thermal  loadings  and  two-  and  three-dimensional 
states  is  expected  to  require  significant  modifications,  possibly  including 
a  oomplete  reformulation. 


166 


w  V ^ V* ■'vIP’'r-T-1  *■  v-  '-*.  T. '  v-' .  ■.»  ’  ,  V  V’> .' 


j-n 


r. 


M 


.v 


K 


K 


>• 


4.4  R.  SCHAPERY’ s  NONLINEAR  STRESS-STRAIN  LAW 
4.4.1  Original  Model 

Dr.  R.  Schapery  used  the  constitutive  theory  advanced  by  him  for  viscoelastic 

17  1ft 

materia  Is  with  mi crocracking  '<  lo  was  taken  by  Dr.  R.  Schapery  as  the  starting 
point  for  predicting  the  response  of  solid  propellants  under  general  loading  conditions, 
The  one-dimensional  version  of  this  law  takes  the  following  simple  form: 


Ap 

0  z:  -  n 

x 


(59) 


where  is  the  linear  viscoelastic  stress  for  a  thermorheologically  simple  mate¬ 
rial  : 


°t  = 


/ 


de. 


E  (£  -  £')  — a-  dr  , 
d  r 


(60) 


with 


t/i  =  t 


u  (T  -  T0)  =  strain  due  to  mechanically  applied 
stress 


$  =  j  dtVAT  [t  ( t ' )  ] 


s  HO 

E(0  =  linear  viscoelastic  relaxation  modulus 


temperature  at  t  =  0 


Ap  =  Ap  (T)  =  temperature-dependent  material  function 
A  =  A  (S^):  softening  function  in  which  the  damage  parameter: 

q 


St  = 


aF 


(61) 
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depends  only  on  the  strain  and  temperature  histories,  and: 


1 


f  =  l  (e/<M  Y 


0 


(e2/el) 

for  constant  threshold  strains  e-j  and  e 


for  0  ^  e  i 

for  e  i  <  £  <  e  2 
for  e  >  e2 

2,  with  P>0. 


The  function  F  =  F(e0)  and  the  positive,  constant  exponent,  q,  originate 
with  the  equation  for  microcrack  speed, 

dA  q 

—  =  MOCr/f) 

d!  (62) 

where  M  is  a  positive  constant  and 


d£  =  dt/Ac 


(63) 


in  which  AQ  =  Ac(T)  is  the  shift  factor  for  microcrack  growth  rate. 

The  functional  form  of  the  softening  function,  A  =  A(Sg),  depends  on  the  type 
of  behavior  that  needs  to  be  reproduced.  The  following  special  case  was  used: 

*  =  [i*=se]p/q 

where  c  and  p  are  positive  constants.  Note  that  when  Sj,  =  0  or  c  =  0,  a  linear 
viscoelastic  stress-strain  equation  is  recovered  from  equation  (59). 

Taking  Ap  =  1 ,  several  sets  of  numerical  values  for  the  constitutive  parame¬ 
ters  corresponding  to  TP-H1011  were  tried  without  success.  This  theory  was  also 
used  to  predict  the  response  of  UTP-19,360  and  UTP-3001.  Having  failed  to 
perform  better  than  linear  viscoelasticity  in  many  cases,  it  has  undergone  sev¬ 
eral  changes  since. 

M.4.2  Current  Model 

The  essential  form  of  the  modified  uniaxial  stress-strain  relation  is  given 

by: 
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i 


o  = 


(64) 


where: 


f  (e°,  <  J  ,  S) 


o  =  engineering  stress 
f °  :  psuedo  strain 
t 

1 

«.  o 


■T,  / 


dt 

E  (t  -  t) - dr 

dr 


(65) 


t°  =  maximum  value  of  K  0  I  up  to  the  current  time 


S  =  damage  parameter 


S  = 


( /Vi 


1/q 


dt 


(66) 


Er  =  arbitrarily  selected  reference  modulus, 
E(t)  =  linear  viscoelastic  relaxation  modulus, 

=  Ee  +  E2  t-n  ■  E2  (Et  +  t"n) , 

ET  -  Ee/E2, 


(67) 


i. 

* 

V; 

i; 

E 


V 


V 

V 


and 


q  =  positive  constant. 

The  functional  form  of  f  in  equation  (64)  depends  on  the  material  considered. 
Studies  on  solid  propellant  to  date  indicate  it  may  be  taken  as  follows  for  some 
solid  propellants: 


f  =  Yi  Y2  Y3  P 15  sign  (f°) 


(68) 
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in  which 


1 

if 

e°  >  0 

sign  (e°)  =  - 

0 

if 

eo  =  0 

-1 

if 

eo  <  0 

and  P15  is  used  to  normalize  function  Yg  to  unity,  at  a  reference  point 
Yj^'s  are  the  following  functions  of  damage  and  pseudo  strain: 


Y-j  =  Y,  (S) 


1  +  A^S  +  for  S  ^  S0 

A4  Sa5  for  S  >S0 


*2 


A2  S0.63-Sx  (cg)  (0.463-Mx-Lx) 


Y3  =  CQ  +  C^x  +  C2x2  +  C$x3  +  C4X14  +Cgx5 


where: 


in  which  xr  is  the  only  root  of  the  equation: 

max(Sr)  =  Y3  (xr) 


.  The 

(69) 

(70) 

(71) 

(72) 

(73) 


with  tux  (Sr)  representing  the  maximum  value  of  Sr  up  to  the  current  time,  and: 


while  X  is  a  factor  that  accounts  for  relatively  small  higher  order  effects 
possibly  due  to  rehealing  and  particle  interaction. 

The  constants  entering  the  definitions  of  Y-j,  Y2,  and  Yj  depend  on  the  mate 
rial.  For  UTP-19,360B  they  are: 


42 

0.637 

-0.387 

0.85 


and  the  factor  X  is  given  by: 


X  =  Kx  C 


x 


The  resulting  form  of  equation  (64)  for  UTP-19.360B  is  thus: 


o  =  P15  Ag  Yi  Y3  f°  sign  (e°) 


Clearly,  if  Lx  =  1,  equation  (77)  may  be  written  as: 


0  =  Ap  |  E  ( t  -  t)  -  dr 

/  dr 


in  which  Ap  =  Ap(e°,e°m,  S)  plays  the  role  of  a  softening  function,  remininscent 
of  the  Mullins-Tobin  approach. 
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4.4.3  Stress  Predictions 

The  current  version  of  the  nonlinear  model  developed  by  R,  Schapery  may  be 
used  to  predict  the  response  of  solid  propellants  with  a  rather  remarkable 
degree  of  accuracy,  as  may  be  seen  in  Figures  102  to  109,  which  are  sample  cases 
of  the  isothermal  tests  considered  in  the  program.  The  first  two  plots  (Fig¬ 
ures  102  and  103)  correspond  to  the  highest  and  lowest  constant  rate  tests  (Test 
No.  1)  available  in  the  data  for  which  the  difference  between  theory  and  test  is 
greatest.  Figure  104  shows  the  saw-tooth  test  (Test  No.  5)  at  constant  rate 
with  increasing  strain  peaks.  Figures  105  and  106  pertain  to  the  dual-rate 
tests  (Test  No.  3)-  Results  for  the  short-  and  long-duration  similitude  tests 
(Test  No.  12)  are  given  in  Figures  107  and  108,  and  Figure  109  includes  a 
three-step  relaxation  test. 

Finally,  it  is  important  to  mention  that  a  complete  characterization  of 
UTP-19.360B  was  also  carried  out  using  Lx  =  1  (the  value  leading  to  equation 
(78)),  and  the  ensuing  response  predictions  were  very  close  to  those  obtained 
with  Lx  =  0.85;  only  the  iow-to-high  dual-rate  test  of  Figure  8  was  predicted 
somewhat  better  with  Lx  =  0.85. 

4.4.4  Material  Ch;  acterization 

In  evaluating  the  material  cor  tants  and  property  functions,  the  following 
observations  may  be  valuable: 

•  Y1(  its  variation  being  brought  about  by  vacuole  formation,  appears  to 

be  a  decreasing  and  concave  down  function  of  damage,  as  presented  in 
the  following  figure: 
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The  function  Sr,  which  provides  a  certain  measure  of  damage,  increases 
as  a  direct  result  of  a  reduction  in  the  number  of  polymer  chains 
supporting  the  internal  stresses;  the  larger  the  Sr,  the  higher  the 
stress  on  each  chain. 

For  UTP-19,3oOB,  the  state  of  damage  is  essentially  constant  during 
unloading  and  reloading,  and  the  shape  of  the  so-called  backbone  curve 
resembles  the  stress-strain  curve  for  rubber,  which  is  of  the  form: 


(€0)1, 


resembles  a  secant  modulus  and  that,  for  most  tests  of  UTP-19.360B  Lx 
=  0.85  produced  a  finite  limiting  value  as  e°  approached  zero,  it  is 
suggested  that  Lx  be  found  in  this  fashion  for  other  propellants. 

•  For  constant-rate  tests,  one  has: 


6°  =  e° 
m 

X  =  1 


and  thus,  from  equation  (72): 


X  =  Xr 


Sr  =  max  (Sr) 


from  which 


Y3  «  Sr 


Equation  (77)  then  reduces  to: 


0  =  1.861  Yt  S°*637  (eO)0.463 


with  Yi  =  Yi  (S)  given  by  equation  (69), 


For  very  small  damage: 


Y1  (S)  =  1 


so  that  equation  (79)  becomes: 


o  -1.861  S0-637  (eo)0.463 

m 

in  which  the  stress  increases  with  damage,  probably  because  of  molecular 

chain  stiffening  due  to  an  increase  in  stress  per  chain. 

With  the  foregoing  observations  in  mind,  determination  of  the  material 
properties  can  be  accomplished  as  follows: 

(1)  The  exponent,  n,  appearing  in  the  relaxation  function,  is  obtained 
from  relaxation-modulus  data. 

(2)  The  normalized  coefficient,  Er,  entering  the  relaxation  modulus,  is 
determined  to  make  unloading  curve  2  in  the  figure  above  pass  through 
the  origin. 

(3)  The  exponent,  q  (present  in  the  definition  of  the  damage  parameter), 
is  evaluated  using  equation  (80)  and  two  constant-rate  tests  at  small 
values  of  damage. 

(4)  The  function  Y1  is  obtained  by  curve-fitting  equation  (79)  to 
constant-rate  tests  over  all  strains  out  to  failure. 

(5)  Experience  to  date  indicates  that  the  function  I2  independent  of  S 
and  em,  and  therefore  equation  (73)  nay  be  used  instead  of  the  more 
general  form  of  equation  (80). 

(6)  The  backbone  curve  Y3  is  determined  using  unloading  and  reloading  data 
like  those  available  in  a  cyclic  test  whose  first  peak  strain  is  the 
largest. 

(7)  Finally  the  correction  factor,  X,  can  be  ascertained  from  a  relaxation 
test  at  a  large  strain  level. 

4.4.5  Multiaxial  Generalization 

A  micromechanics  model  has  been  developed  which  predicts  the  form  of 
equation  (73),  but  its  use  was  not  pursued  on  this  program. 


4.5  M.  GURTIN'S  THEORIES  FOR  NONLINEAR  VISCOELASTIC  MATERIALS 

Four  essentially  different  approaches  have  been  followed  by  M.  Gur'  in 
in  trying  to  predict  the  response  of  solid  propellants  that  exhibit  damage. 
The  stress-softening  theory  appears  to  be  the  most  accurate  of  the  four  laws 
as  will  be  pointed  out. 

4.5.1  Original  Model 

The  one-dimensional  stress-strain  law  for  materials  undergoing  internal 
damage  was  based  on  the  hypothesis  that  the  state  of  damage  at  any  time  is 
completely  characterized  by  the  maximum  strain,  em,  that  the  material  has 
experienced: 


em  (t)  =  max  e(s) 

0  <  s  <  t 


(82) 


The  stress,  a,  is  given  by  a  constitutive  equation  of  the  form: 

o  (t)  *  g[e(t),  em(t)]  (83) 

and  it  depends  only  on  the  current  values  of  strain  and  damage.  Such  an 
equation  is,  of  course,  rate-independent. 

In  this  theory,  if  the  maximum  strain  occurs  at  the  present  time,  then 

em  (t)  =  e(t),  (84) 

and  equation  (82)  reduces  to: 

a  =  G(em)  =  g(ero,  era)  (85) 
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The  stress-strain  curve 

(86) 

o  =  G(em) 

is  called  the  virgin  curve  and  is  traced  out  in  an  experiment  with  monotonically 
increasing  strain. 

Using  the  virgin  curve,  equation  (83)  may  be  rewritten  in  the  form: 

(87) 

o  =  F(*,em)  G(em) 


with: 


*  =  — 
e  _ 


the  relative  strain,  and: 


F(*,*m)  = 


g(  cm*  cm^ 
G(em) 


The  function  F($,  e  m)  is  called  the  damage  curve  at  the  damage  level  em,  and  is 


such  that: 


F  ( 1  >  e  m )  =  ^ 


In  some  situations  of  interest  F(£,  em)  is  independent  of 


F(*,e_)  =  F<*> 
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When  this  is  so,  F(£)  is  referred  to  as  the  master  damage  curve,  and  equation 
(87)  reduces  to 


o=  F(f)G(em> 


(92) 


As  pointed  out  previously,  this  is  a  rate-dependent  theory,  and  as  such,  cannot 
be  used  for  loading  rates  that  differ  much  from  that  used  to  determine  the 
damage  function.  This  situation  was  remedied  by  changing  the  stress-strain  law 
to  the  one  described  next. 

4.5.2  Nonlinear  Model  Based  on  Stress  Softening 

To  develop  a  simple  theory  of  stress  softening  which  allows  for  rate 
effects  and  which  returns  to  Mullin' s  original  ideas  of  using  the  past  stress 
maximum  as  the  damage  parameter,  two  fundamental  ingredients  are  considered. 

The  first  is  the  virgin  stress,  S,  which  represents  the  stress  the  material 
would  experience  in  the  absence  of  softening.  This  stress  is  assumed  governed 
by  a  constitutive  equation  of  the  typu  encountered  in  linear  viscoelasticity. 
The  second  ingredient  is  a  damage  function,  F,  which  gives  the  true  stress,  o, 
when  the  virgin  s'.ress,  S,  and  its  past  maximum,  Sm,  are  known. 

The  one-dimensional  form  of  the  constitutive  law  for  a  classical  linear 
viscoelastic  material  is  given  by: 

t 

a  (t)  =  Jg  ( t 

—OO 

in  which  o(t)  is  the  stress;  e(t),  the  strain;  and  G(t),  the  relaxation 
function.  It  is  further  assumed  that  f(t)  =  0,  prior  to  t  =  0. 


The  generalization  of  equation  (93)  is  begun  by  defining  the  quantity: 


S(t)  =  G(t  -  7)  t(r)  dr 


which  is  called  the  virgin  stress  and  which  represents  the  stress  that  would  be 
present  in  the  absence  of  softening.  It  is  assumed  that  the  extent  of  softening 
is  governed  by  a  constitutive  equation  giving  the  true  stress,  o(t),  when  S(t) 


and  its  past  maximum  are  known: 


Sm(t)  =  max  S  (7) 


0  <  r  <  st 


Without  loss  of  generality,  this  constitutive  equation  is  written  in  the  form: 


0  ~  ^m  F  (  »  ^m 

\  Sm 


and  it  is  assumed  that  the  damage  function,  F,  satisfies  the  following 
conditions: 


F(1,  Sm)  =  1 

F(x,  Sm)<  x  for  x<1 


These  restrictions  imply  that: 


0  (t)  <  SCt) , 


also  that: 


am  (tO  =  ( t) , 


(99) 


and,  that  the  following  conditions  are  equivalent: 


i)  o(t)  =  S(t) 

ii)  S(t)  =  SmCt) 

iii)  o(t)  =  om  (t) 


(100) 


where  om  is  the  past  stress  maximum,  defined  analogically  to  Sm.  The  inequality 
equation  (98)  asserts  that  the  material  actually  softens,  while  equation  (99) 
indicates  that  this  softening  occurs  when  and  only  when  S(t)<Sm(t)  lor 
equivalently  a(t)<om(t).  The  results  of  equation  (98)  and  (99)  show  that 
one  may  equally  well  use  the  true  stress,  a(t),  as  the  damage  parameter. 

Equation  (97)  and  the  fact  that  l  .e  first  relation  of  equation  (100)  implies 
the  third  are  direct  consequences  of  the  hypotheses  laid  down  in  equation  (101). 
To  verify  equation  (98),  note  that  if  the  maximum  of  S  on  the  interval  o<T<t, 
occurs  at  r  =  {  ,  then: 


which,  by  virtue  of  (100)  and  the  definition  of  Sffi,  implies  that: 

n(()  =  S($)  =  Sm  (t)>S  (>)><’  (A);  0<  A  «■'  t 
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proving  equations  (97)  and  (98)  and  the  first  two  equations  of  equation  (100) 
To  establish  the  third  equation  (100),  note  that  if: 

Sm(t)  =  S(t) 

which  implies  that: 


Sm(t)  =  o(t) 


because  of  equation  (101);  then,  since: 


Sni  =0m 


one  would  have  that: 


om  (t)  =  o(t). 


Conversely,  if 


a(t)  =  om(t), 


then: 


=  om(t)  =  Sn,(t) 


so  that: 


S(t)  =  ^(t). 


A  _\  A  *.  *.  *. 


(104) 


A  V* 


Returning  to  the  constitutive  equation  (95),  it  is  interesting  to  consider 
the  special  case  in  which  the  damage  function  depends  only  on  S/Sm: 


(105) 


which  is  a  master  damage  curve  of  the  type  considered  in  the  rate-independent  model 
discussed  previously. 


When  the  virgin  stress  obeys  an  elastic  stress-strain  relation: 


S  =  Ec 


(106) 


then: 


-"iii  "  E^m 


(107) 


in  which  Cjj  is  the  past  strain-maximum,  and  equation  (95)  yields: 


0  =  Eem  *  Eem) 


(108) 


so  that,  defining: 


F*  (e  »  e“)sF(e  ’  Ee®) 

'  m  '  '  m  ' 


leads  to  the  starting  assumption  of  Gurtin  and  Franc is^ 


(109) 


* =  '") *■ 


(110) 


presented  earlier  as  the  rate- independent  model. 


Although  implicit  in  equation  (95)  is  the  assumption  that  the  functional 
form  of  F  would  be  the  same  for  unloading  conditions  as  for  reloading,  it 
was  found  experimentally  that  a  different  damage  function  is  needed  for  each 
of  these  processes.  Actually,  there  is  more  than  one  way  of  obtaining  the 
same  damage  function.  For  TP-H1011,  for  instance,  the  following  procedure 
was  employed. 

Considering  the  strain  history  shown  in  Figure  110,  on  the  loading  portion 
we  have: 


G(t)  At  for  t<  T 


(111) 


hence,  S(t)  increases  monotonically  and,  by  equation  (98): 


Sm(t)  =  S(t)  =  o(t) 


(112) 


and,  upon  unloading,  the  past  maximum  of  S  is  the  true  stress: 


o  m  =  c(T) ,  T  <  t  <  2T 


(113) 


Figure  110.  Strain  History  Used  to  Characterize  the  Damage  Function 
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Further,  by  equation  (93): 


=  u,|/ 


G(t  -  T)  d T 


v» 

•/ 


G(t  -  r)  dr 


(114) 


or,  equivalently: 


S(T  +  t)  =  G  (t)  -  o(t) 


(115) 


t  +  T 


G(t)  s 


*'/ 


G  (r)dr 


(116) 


Letting  o.,(t)  and  o2(t)  denote  the  true  stress  during  loading  and  unloading, 
respectively,  with  t  in  02(t)  me^-ured  from  the  time  T  at  which  unloading 
begins,  equations  (95)  and  (110)  yield  the  simple  formula: 


02(t)  j  G(t)  -  o^t) 


(117) 


Thus,  summarizing,  the  stress-softening  approach  to  damage  1s  described  through 
the  following  constitutive  equation: 


(t)  =  SB  F^,  S.) 


(118) 


O  1^0  O  fc'fv*  O  L^V'  V*  «-*  O.  •/  «c*  I.  *  Sit-  *  •  •  *  «  '  V*  »_•  •  •  c‘,  %_  -  V  •J'V.W"  0*0  O. ' 


1V1VWI 


where 


S(t)  = 


G(t  -  t)  —  (r)  dr 
dr 


and  in  which  the  damage  function,  F,  may  be  determined  from  sawtooth  tests  with 
increasing  peak  strains  and  with  sufficiently  long  rest  periods  between  pycles. 
For  conditions  of  reloading,  F  is  given  by: 


a- 1  (t) 

—  =  F  -  ,  °m 

0  m  L°m  J 


(119) 


while  for  unloading  the  following  form  is  employed: 


a?  fG(t)  -o-i(t) 

-:F  - ,  a 

am  _  m 

L  °  m 


(120) 


Typical  curves  for  unloading  and  reloading  damage  functions  of  TP-H1011  are 
included  as  Figures  111  and  112,  respectively. 

Finally,  we  point  out  that  the  use  of  this  approach  to  predict  the  response 
of  TP-H1011  yielded  results  that  were  far  more  accurate  than  those  obtained 
with  any  of  the  other  theories  in  their  original  form. 

4.5.3  Nonlinear  Models  Based  on  Maximum  Strain 

A  series  of  constitutive  relations  based  on  the  past  maximum  strain  have 
been  proposed  by  M.  Gurtin.  The  precursor  of  these  relations  took  the  form: 


(121) 


G(t  -  r)  —  (r)  dr 
dr 


0  0.1  0.2  0.3  0  4  0.5  0  6  0  7  0  8  0  9  1  0 

G(t)-  c,(t)/em 

Figure  111.  Damage  Function  for  Unloading  ( TP  — H 1011) 
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where  G  represents  the  relaxation  modulus,  and  the  function  e  was  expressed  as 
a  product  of  the  reduced  damage  function,  F,  and  the  virgin-response  function, 
g,  In  the  following  way: 


e  =  F 


»  emj  8  (fm) 


U22) 


with: 


F ( 1 » =  ) 


(123) 


Figure  112.  Damage  Function  for  Reloading  (TP-H1011) 


and 


em  =  max  e(r) 
0<  t<  t 


so  that,  during  virgin  response,  since: 


one  had: 


o(t) 


t 


/ 


G(t  -  T) 


dg(e) 

de 


dc(r) 

-  d  t 

d  r 


(126) 


and,  by  taking 


k  =  K 

g(e)  =  ^  Ak  ek 

k  =  1 


(127) 


the  best  values  for  the  a^'s  could  be  determined  using  least  squares  and  the 
data  of  all  constant-rate  tests. 

To  characterize  the  reduced  damage  function,  F,  involved  the  determination 
of  a  creep  function,  J,  solution  o; • 


/ 


dJ(r) 

G  ft-  T)  -  dr  =  1 

d  r 


and,  such  that: 


e(t) 


d a  (r) 

t)  -  d  r 

d  r 


(128) 


(129) 


Thus,  taking  the  reduced  damage  function,  F,  in  the  form: 

F  (x,  y)  =  F,  (x)  F2  (x,y)  (130) 


M  -  1 

Ft  (x)  =  xM+  '2  d®  (x“  "  xM) 
o=1 


(13D 


F2  (x,y) 


p  -  Q  -  1 

1  +  bp  yP  [x  -  xQ  ♦  ^  cq  (  q  ”  Q>  ]  (132) 

p  =  1  q  -  2 


and  equating  equations  (114)  and  (118),  the  coefficients  entering  F  were  to  be 
determined  using  least  squares  and  all  the  saw-tooth  data  with  increasing  strain 
peaks . 


When  this  constitutive  law  was  applied  to  UTP-19.360B  data,  it  was  deemed 
necessary  to  change  the  form  of  the  function  e,  because  of  the  large  errors 
observed  in  the  predicted  response. 


The  last  of  a  sequence  of  modifications  yielded  the  following  stress-strain 

law: 


o(t) 


(133) 


where,  as  before,  G  was  the  relaxation  modulus,  and: 


F  (1,eB)  =  1 


(134) 


em  =  max  e(r) 


v 


0<7<t 

In  the  present  case,  the  virgin  response  was  given  by: 


\ 

o(t)  =  | 


de 


m 


G(t  -  r)  K  (em,  cm)  -  dr 

dr 


while,  the  damage  response,  for  which  em  remains  constant,  took  the  form 


t  ( 

f  3F 

o(t)  =  K  (€m,  0)  }  G(t  -  r)  —  (x,y)  dr 

j  3  x 

o 

in  which 

K(f ,e)  =  1  +  At  (e  -e0)  +  A2  (*  -  «0>2  +  A3  (e  "  eo>3  +  ( 

+  c  (B1  e  ♦  B2«  2  +  B3  e  3) 


F(x,y) 


D7  y2) 


4 


a(x)  s  x5  +  Y  Dm  (xm  -  x5) 


(141) 


x  s  e/e. 


y  h  em 


(142) 


A  set  of  stress  predictions  obtained  for  UTP-19,360B,  with  the  resulting 
version  of  the  theory,  is  included  in  Figures  113  through  119.  Figure  113  is 
for  the  high-to-low  dual  rate  test,  Figures  114  through  116  are  segments  of  the 
long-duration  similitude  test,  and  Figures  117  through  119  are  segments  of  the 
three  step  relaxation  test. 

Since  the  dependence  of  the  function  K  on  the  strain  rate  was  felt  to  be 
artificial,  the  treatment  of  damage  was  revised  in  the  manner  explained  next. 

4.5.4  Current  Model 

The  latest  version  of  M.  Gurtin’s  nonlinear  stress-strain  law  is  based  on. 
a  strain-dependent  relaxation  function  and  has  the  form: 


o(t)  =  J  G  j  e  (r),  r  Je(t  -  t  ) 


(143) 


where 


G(e  ,  t)  =  Gr  (t)  Gc  (e  ,  t) 


G  =  relaxation  modulus 


(144) 

(145) 

(146) 


Gc  =  correction  modulus,  defined  as: 


(e,  t)  =  ^  An  (e)  (< 


‘n(0=  X 


^np  e 


For  this  material,  the  virgin  curve  (e  =  em  =  e )  is: 


o  =  <?r  ♦  ac 


with  o  the  linear  viscoelastic  stress 


■I 


or(t)  =  I  Gt(t)  c  (t  -  t)  dr 


and  the  correcticn  stress,  oc,  given  by: 


(t)  =  J"  G c[6(t)>t]  e(t-r)dr 


0  49; 


Hence,  to  characterize  the  virgin  response,  only  constant  rate  tests 
need  be  employed.  In  this  instance ,  o and  ap  are  known,  so  that  form  equation 
(131)i  °c  nay  be  computed,  and  equated  to  equation  (133)  using  the  fact  that 
€  is  a  constant;  i.e.: 


o(t)-or(t)=oc(t)  s  j  G  |f  (X),  X  |  d  X 


(  1r  0  ) 
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Relaxation-Unload-Reload  of  6-in.  Bar  of  UTP-19 , 360B-400/1777 


Figure  117.  Three-Step  Relaxation  of  6-in.  Bar  of  UTP-19 , 360B-400/1777 


Relaxation  of  6-in.  Bar  of  UTP-19 , 360B-K00/1777 


which,  upon  recalling  equations  (134)  and  (135)  becomes: 


N 

o(t)  -  0r  (t)  =  ^ n  (t) 

n  =  1 


where: 


*n  <t) 


P  t 

X  Anpe*  J  *P<X>  e 

p  =  1  o 


Furthermore,  since  for  a  constant-rate  test 


e  (X)  =  cX 


it  follows  that: 


*»  («  =  ^ 


P+1 


Anpe 


P  *  1 


and  after  integrating  by  parts: 


with 


fo  (x)  =  i  _  e-x 


(156) 


fp  (x)  =  -xR  e~x  +  pf  1  (x)  ;  p  =  1,  ...I 


Clearly,  equations  (132)  and  (135)  to  (137)  may  be  used  to  determine  the 
coefficients  anp  appearing  in  the  definition  of  the  correction  modulus.  The  pro¬ 
cedure  suggested  by  M.  Gurtin  to  accomplish  this  is  as  follows: 

1.  Take  N  tests  with  constant  rates  e2,  . ..  eN;  and  set: 


T i  =  “  5  i  =  1 i  N 

e  . 


(157) 


2.  Select  the  degree,  P,  of  the  series  expansion  of  the  correction  modu¬ 
lus,  as  it  appears  in  equation  (130). 

3.  Use  the  test  and  the  approximation: 


oe  (t)  s  (€1t  t) 


(158) 


to  find  the  a-|p. 


4.  Use  the  e 2  test  and  the  approximation 


oc  (t)  -  (e2,  t)  =  *2  <e1> 


(159) 


to  find  the  ajp. 

5.  Use  the  e  3  test  and  the  approximation 


ac(t)  -  ^  (f3,  t)  -*2  (e 3,  t)  =^3  (t3,  t) 


(160) 


to  find  the  a3p,  and  so  on  for  the  s^p^.a^ 


V.n"SS 


6.  Iterate  this  procedure  if  necessary;  that  is,  define 


N 


*(e,  t) 


2  *n  («. 


t) 


n  =  1 


so  that  4*  is  known.  For  each  ramp  test,  define 


°c  *  °c  ~ 


(161) 


(162) 


and  repeat  the  above  procedure  using  oQ  to  find  constants  anp.  The  new 
values  of  the  a^  are 

,  .  -  (163) 

'anp'new  =  anp  +  anp 


7.  Repeat  the  process  if  nec  ".‘ary. 

It  is  important  to  point  out  here  that  numerical  difficulties  may  be 
encountered  in  applying  this  technique  to  characterizing  the  virgin  response  of 
solid  propellants.  In  fact,  some  convergence  problems  were  faced  in  connection 
with  the  OTP-19, 360B  data.  Moreover,  characterization  of  the  damaged  response 
calls  for  a  large  number  of  cyclic  tests  over  a  wide  range  of  rates.  This 
increases  the  convergence  difficulties.  The  model  was  employed  with  the 
constant-rate  tests  only  for  this  reason.  Figures  120  to  123  show  the  results 
of  the  stress  predictions  obtained  with  the  current  version  of  the  model.  The 
first  two  plots  correspond,  respectively,  to  the  lowest  and  highest  rates 
available  at  ambient  temperature. 

So  far  in  the  program,  none  of  the  models  developed  by  Gurtin  have  taken 
into  account  the  effects  of  temperature  on  propellant  response.  However,  the 
tlae-temperature  superposition  principle  was  tested  with  the  current  version  of 
the  theory.  The  results  appear  in  Figures  122  and  123  for  the  thermal  tests  at 
123  and  40  F,  respectively.  The  use  of  the  superposition  principle  breaks  down  at 
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Nonlinear  Viscoelastic  Stress  Predictions  for  UTP-19 , 360B-U00/1777 
at  70  F  (M.  Gurtin's  Theory) 


Nonlinear  Viscoelastic  Stress  Predictions  for  UTP-19,3&0B-1J00/1777 
at  123  F  (M.  Gurtin’s  Theory) 


a  strain  of  about  39X.  This  was  apparently  due  to  the  ambient  temperature  data 
base  being  limited  to  a  low  strain  level.  Also,  It  minht  be  necessary  to  allow 
the  relaxation  fund  ion  (G(, )  to  depend  on  the  glass  fr.inst  '  on  temperature  (Tg) 

Ge  =  Gc  (e,  1^,  t) 

H.6  RUSSIAN  APPROACH  TO  PHYSICALLY  NONLINEAR  VISCOELASTIC  SOLIDS 

The  Russians  have  explored  two  approaches  for  characterizing  damage  effects 
in  solid  propellants.^*  ®  They  are  a  general  functional  approach  and  a  kinetic 
equation  of  evolution  for  damage.  Both  approaches  are  based  on  in temalvariable 
concepts,  and  either  approach  appears  general  enough  to  also  incorporate  cumula¬ 
tive  damage  and  propellant  response  under  multiaxial  stress  states.  However, 
the  general  functional  apDroach  may  be  of  little  practical  engineering  value 
because  a  very  extensive  testing  program  may  be  required  to  evaluate  material 
parameters.  This  approach  requires  introduction  of  damage  measures  which  should 
reflect,  rnicrostructural  damage  mechanisms,  and  a  damage  functional  which  charac¬ 
terizes  the  accumulation  of  damage  or  defects.  The  damage  functional  is  then 
expanded  into  a  series  of  multiple  integrals  in  an  analogous  fashion  to  that 
followed  by  Green  and  Rivlin  for  nonlinear  materials  with  fading  memory.  Herein 
lies  the  difficulty.  Even  assuming  isotropy,  four  to  six  different  tyres  of 
multi-axial  tests  are  required  to  evaluate  the  required  material  property  func¬ 
tions  for  a  first-order  theory.  Although  the  approach  has  theoretical  merit  and 
may  even  have  some  practical  application  in  the  future,  its  pursuit  was  aban¬ 
doned  in  favor  of  the  kinetic  approach. 

Tne  essential  feature  of  the  kinetic  approach  is  to  introduce  the  degree  of 
damage  into  the  constitutive  equations  as  a  reduced-time  parameter  in  the  same 
way  that  temperature  is  introduced  as  a  reduced-time  parameter  for  the  therroo- 
rheologically  simple  materials  in  linear  thermoviscoelas'ricity .  Damage  is  then 
defined  in  terms  of  some  strength  parameter  of  the  material,  and  the  degree  of 
damage  is  characterized  through  an  equation  for  damage,  as  explained  subsequently. 

^.6.1  Original  Model 

The  one-dimensional  constitutive  equation,  taken  from  the  Russian  litera¬ 
ture  by  W.  L.  Hufferd  as  a  means  of  predicting  the  response  of  phys  i sally 
nonlinear  viscoelastic  materials,  may  be  expressed  by: 


a  (t)  =  J 


E  (t‘  -  r')  — (r)  dr 
dr 


(164) 


where 


a  =  stress 

e  =  strain  due  to  mechanically  applied  stress 
E(t)  =  relaxation  modulus 
E(t)  =  Ee  ♦  E2fn 


t’  -  r* 


(165) 


aTj  [tj  cr>] 


represents  the  daf ige-reduced  time,  which  is  arrived  at  in  the  manner  described 
next. 

First,  a  normalized  damage  function  u>  =  u>  (t)  is  introduced  through  the 
following  kinetic  equation  of  evolution 


(t)  =  h  (w)  f(t) 


(166) 


in  which  it  is  further  assumed  that: 


f(t)  =  I  F(  t  -  t)  4>  |  o0(  r)  J d  r 


(167) 
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together  with  the  conditions: 


tJ(o)  =  0  (168) 

cj  ( t* )  =  1 

indicating  that  no  damage  exists  at  the  initial  state,  and  that  failure  occurs 
at  time  t*. 

Next,  equation  (166)  is  integrated  with  u>(o)  =  0,  leading  to 


co 


f—  ■  f 

J  h(w)  I 


f (r)  dr 


(169) 


h(u>) 
o  o 


Setting  t  =  t*,  so  that  (t*)  =  1,  and  substituting  equation  ( 1 67 )  for  f(r),  it 
is  obtained  that 


t«  £ 


d?  J  F(£  -  r)  <t>  JoQ  (r)Jdr 


=  1 


(170) 


i 

J 


d  u> 
h(cu) 


If  now  the  function  F(t)  is  assumed  to  have  a  power-law  representation 


Equation  (170)  can  be  written  in  the  form 


j  d£  j  F0(£  -  t  )m^|a0(r )  j  dr 


o  o 


1 


I 


0 


d  cu 
h(w) 


and  integrating  with  respect  to  £,  assuming  that  the  order  of  integration  may  be 
interchanged,  one  arrives  at 


which,  for  the  case  where  o0  and  0  [o0]  =  0[oo]  are  constant,  becomes 


Fo0o^o>  <t*)  2+” 

(1  +  m)  (2  +  m)  1  ”  1  (174) 

d  u> 
h(a>) 

where  t*0  is  the  time  to  failure  under  the  constant  stress  o0.  Thus  equation 
C 1 73 ) »  in  this  case,  may  be  written  as: 


:?  +  m)  j  (t#  - 


f  1  +  m  d  r 

(t!)2+m 


(175) 


If  the  time  to  failure  under  a  constant  stress,  oQ,  has  the  power-law 
representation 


w 

a0a  t0  -  constant  =  (3 


(176) 


then  equation  (175)  can  be  put  in  the  form: 


j  (t*  - 


r)  1+m  o0^2+m)  dr  =  ^ 

2+m 


(177) 


so  that,  motivated  by  equations  (175)  and  (177),  the  degree  of  damage 
accumulation  may  be  introduced  through  the  expression 


C  1+m 

r?(t)  =  (2  +  m)  (t  -  r)  dr 

J  (to)2*" 


(178) 


in  which 


V  (o)  =0 

V  (D  *  1 


(179) 


-'.V-V.V.s 


The  function  tj(t)  can  be  related  to  the  damage  function,  u>,  by 


I 


(180) 


This  means  that  tj  represents  the  relative  damage  in  the  load  history 
for  the  power-law  representation  of  t#Q.  From  equations  (176)  and  (177), 
equation  (178)  may  be  written  as: 


2  +  m 


V  (t)  = 


(t  -  t) 


1  +  m  ^  a(2  +  m) 


(181) 


and  finally,  the  influence  of  damage  is  treated  as  a  reduced  variable  by 
introducing  the  modified  time,  t' ,  defined  by: 


atj  h(t) 


(182) 


on  which  equation  (165)  is  based  and  where  the  shift  function  due  to  damage,  a^, 
depends  on  the  material  at  hand. 

4.6.2  Hufferd's  Modification  to  Il’yushin  Theory 

One  revised  version  of  the  Russian  stress-strain  law  takes  the  form: 


°M- 1  E(ir)7> 


(183) 


where  is  a  damage-related  shift  function,  assumed  to  depend  only  on  the 
current  state  of  strain;  specifically: 


i  • 

=  arj  e(t),  t  (t) 


084) 


Clearly,  if: 


E(t)  =  Ee  +  E2f 


085) 


then  equation  056)  becomes: 


o( t)  =  Ee  e(t)  +  (a^  )"  E2 


(t  -  t)~  — (t)  dr 
dr 


086) 


which  resembles  the  classical  approach  of  the  softening  function  used  as  a 
stress-correction  factor. 

Another  revised  version  of  the  Russian  approach  consists  of  retaining  most 
aspects  of  the  original  law,  but  constant  strain-rate  data  are  employed  to 
express  the  time  to  failure  as 


t-H 


087) 


and  equation  078)  is  changed  to 


*?(t)  = 


(2  +  m)  C 

hi  J 


r)UnU(T)  dr 


088) 
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where 


<P  it)  =  e ( t )  =  <  t 


and 


Mt*)  =  1 


Thus,  evaluation  of  equation  (188)  at  t  =  t#  yields 


(2  +  m)  r  . 

- - —  I  (t*  -  r)  +m0(r)dr=  1 

ftQ  J 

o 


(189) 


which,  through  a  change  of  variables  a- d  after  some  algebraic  manipulations,  may 
be  integrated  to 

(190) 

«  -+ni 

(t0)  =  3  +  m 


The  solution  for  m,  as  a  function  of  strain  rate,  is  easy  to  obtain  using 
equations  (187)  and  (190)  (as  presented  in  Figure  124  for  UTP-300). 

In  much  the  same  way,  integration  of  equation  088)  for  the  relative  damage 
function,  V(t),  leads  to 


»?(t) 


(191  ) 
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Hence,  using  constant  strain-rate  data  to  express  the  time  to  failure 
does  simplify  things,  but  a  major  assumption  would  still  be  needed  regarding 
the  form  of  the  damage  shift  function  as  it  was  in  the  original  set  of  equations, 
In  this  context,  it  is  important  to  note  that  the  linear  expression 


a„ =  1  - n 


(192) 


and  the  exponential  form 


an  * 


(193) 


were  used  for  the  damage  shift  function  without  success.  For  this  reason,  the 
modified  version  used  to  run  the  stress  predictions  included  in  this  report, 
corresponds  to  equation  (186). 

4.6.3  Stress  Predictions 

Figures  125  to  129  show  the  comparison  between  the  observed  response  and 
that  calculated  ir;ng  the  present  theory.  As  may  be  seen,  the  predicted 
response  is  quite  accurate  in  all  cases  considered,  which  include  constant-  and 
dual-rate  tests  as  well  as  a  shcrt-duration  similitude  loading. 

4.6.4  Material  Characterization 

As  may  be  gathered  from  equation  (186),  the  simplest  version  of  this  theory 
requires  the  knowledge  of  only  two  materia  1-property  functions,  to  wit: 

1.  The  relaxation  modulus,  and 

2.  The  damage  shift  function: 


i  def 

( a^'  -  a^  =  a^  <  (t) ,  <  (t) 


(194) 


which  is  determined  in  the  following  way3: 


Figure  125.  Nonlinear  Viscoelast 


idure  128.  Nonlinear  Viscoelastic  Stress  Predictions 


a.  From  constant  strain-rate  tests,  to  correct  the  stress  response  during 
loading; 

b.  From  a  relaxation  test  at  a  large  strain  level,  to  account  for  healing; 
and 

c.  From  a  constant  strain-rate  cycle  carried  out  to  a  large  strain  level, 
to  more  adequately  reproduce  the  hysteretic  behavior  of  the  propellant. 

The  damage  shift  functions  corresponding  to  UTP-19,360B  are  shown  in 
Figures  130  to  134.  The  first  three  of  these  plots  represent  typical  curves  of 
a  for  low,  intermediate  and  high  strain-rate  tests,  while  Figures  133  and  134 
give  the  correction  curves  for  relaxation  and  unloading,  respectively. 

4.6.5  Multiaxial  Formulation 

The  previous  forms  of  the  constitutive  equation  were  specialized  from  a 
three-dimensional  theory.  Returning  now  to  this  consideration  and  continuing  to 
treat  the  influence  of  damage  as  a  reduced  variable  by  introducing  the  modified 
time  t' ,  where 


dt '  = - —  —  and 

a,L’(0] 


f-  I  . 

0  ai?fr(o! 


(195) 


The  constitutive  equations  may  then  be  written  in  the  form 


E(t)  =  I(t)  f  J^t’  -  r\  T)  —  (tr  S’) 
J  O 


(196) 


J2(t'  -  r\  T)  37  S(r)  dr 


S(t)  =  I(t)  /  E-| ( t '  -  T’,  T)  A  E(t)  dr 

Jo  3r 

c  ^ 

+  /  Eo( t '  -  r'  T)  -A  E( r)  dr 


(197) 


modified  shift  function  for  damage,  a 


Strain  (  O.  % 

Figure  130.  Constant-Rate  Test  (0.001  in. /min) 
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In  equations  (196)  and  (197),  Ji  and  are,  respectively,  the  creep  and 
relaxation  functions  in  bulk  and  J2  and  E2  the  creep  and  relaxation  functions  in 
shear . 

Assuming,  for  example,  that  a  n  (tj)  has  the  exponential  representation 


Figure  132.  Constant-Rate  Test  (10  in. /min) 


Softening  function 


!.»* 

Figure  133-  Softening  Function  During  Relaxation 
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Figure  1 3*J •  Softening  Function  During  Unloading 
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:ij(t)  =  A0  Sjj(t)  4  B0  /  (f  -r'fj-jSM  dr 


To  characterize  the  dependence  on  the  magnitude  of  the  previous  loading,  and  to 
simplify  the  previous  representations,  stress  and  strain  intensities  defined  by 


=  (^h  s = 


*  (Z/3  E  ■  e) 


n 

1/2 


(204) 


are  introduced.  Equation  096)  may  then  be  written  in  the  form: 


•l(Ee,  trE)E,j*f  (Se,  trS)  SijU) 

♦  /  T(f  -f)  f^Sg,  trS)  Sij(r)  dr 


(205) 


and 


4>2(trE,  Ee)  trE  =  f2(Se,  trS)  trS 


♦  /  U(t’  -f)  f2(Se,  trS)  trS(r)  dr 


(206) 


Equations  (205)  and  (206)  now  include  coupling  between  deviatoric  and 
dilatational  behavior  as  well  as  time-dependent  bulk  response.  Under  certain 
conditions,  it  is  possible  to  write  equations  (205)  and  (206)  in  the  alternate 
form5’  5 

fj(Se,  trS)  Sjj  =  01  (Ee,  trE)  Ejj 

■  t 

R(t’  -r’)  *i[Ee,  trE]  Ejj(r)  dr  (207) 


/' 
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and 


f2(Se,  trS)  trS  =  $2(Ee’  trE 

-  /  V(f  -f)  trE:)  trE(r)  dr 


(208) 


The  functions  f , ,  f 2 .  <t>  i ,  and  02  may  be  represented  by  a  '’rony  exponential 
series  or  the  simple  power  law  forms: 


f:(s  trS)  =r  A.  S Q  •  1  (trS;‘  (i  =  0,  1.  . 

i  l  c 


0l(Ee,trE)  =  I.  Bj  Ee  '  (trE)»  (i  =  0,  1 . n) 

f  (Se ,trS)  =r  C,  (trS^-  •  Se‘  (i  =  0,  1 . p) 


02(^trE)  =  Zj  D,  w.  )6-  1  Ee*  ( i  =  0 ,  1 . 


where  Ait  Bif  Ct  and  Di  are  constants  and  ur ,  , j ,  0  and  6  are  integers.  Practically 

speaking,  the  first  two  or  three  terms  should  be  sufficient  to  adequately 
characterize  the  damage  response. 

4 .6 .5 - 1  Evaluation  of  Damage  Functions  from  Experimental  Behavior 

Characterization  procedures  for  evaluating  the  damage  functions  described 
in  the  previous  sections  are  discussed  in  the  following  paragraphs. 

The  simplest  experiment  is  creep  in  pure  shear.  In  this  case  there  is  only 
one  non-zero  component  of  stress  (e.g.  S^)  and  one  non-zero  component  of 

strain  <E12).  The  appropriate  measures  of  stress  and  strain  intensity  from 
equation  (204)  are 


Jj 


Since  tr  E  =  tr  S  =  0,  equation  (206)  is  satisfied  identically  and  only  the 
non-vanishing  equations  of  (208)  by  means  of  equation  (209)  has  the  form 


where 


B0  - j  E12Ct)  =  A0/^  Si2\  S120(t) 


0  (t)  =  1  4  f  T  (  t  ’  -  T’  ) 


If  we  normalize  the  response  to  the  jth  test,  it  follows  that 


E12(t) 

TTT- 

Ei2(0 


5ai\*  *  “ 


l-Tfl-  2  E(ij  O  r1  4  f 

♦  to  l4«  *  -,-4, — £ — 


Cnl/2^  l  - 


* 


The  value  of  the  ratio  (1  +a)/(1  + /3  )  can  be  determined  from  a  cross  plot  of 
log  creep  strain  versus  log  stress. 

Consider  next  a  constant  strain  rate  test  conducted  under  superimposed 
hydrostatic  pressure.  In  this  case, 


S11  =  So  *  P 


S22  =  S33  =  -  P 

r  _  r  tr  E 
E11  =  Eo  -  — j” 

E22  =  E33  =  +Avp  * 


Se  =  S0 

Ee=  Ul  +  v)  E0 


tr  S  =  -p  - 

3 

tr  E  =  Avr  +  (1  -  2v)  E0 


Here,  EQ  is  the  applied  axial  strain  and  SQ  is  the  resultant  stress;Avp  is  the 
(negative)  volume  change  associated  with  applied  pressure  in  the  absence  of 
applied  deformations.  Assuming  elastic  volumetric  response  and  a  constant 
strain  rate  test,  E0  =  E0t,  it  follows  that 


£  *  "’(Eo)  2  A„/J-P  \k 

k  v  y 

,,  x  1+?  vfo  /  ^Vo\k 

=  E0  (’-Ir.g-fJ 

-  U  .<•)  *fl2Bk  2  (1.,)  5'  k  /  R(f-r’) 

K  l  J  0 

k^.2,^SYt1.»  dtJ' 


=  (1  ♦  P  )  E„  2 


Av  \  k 


(216) 


The  function  SQ  is  determined  experimentally  for  various  strain  rates  and 
pressures.  The  functions  R(t)  may  be  determined  from  simple  relaxation  tests. 


Equation  (216)  then  represents  a  system  of  linear  equations  for  determining 
the  coefficients  Ak  and  B^. 

More  generally,  it  is  possible  to  evaluate  the  functions  f1  and  <*>i  in 
equation  (207)  from  other  tests,  assuming  time-independent  bulk  response.  If  we 
let  S  =  tr  S  and  E  =  tr  E,  a  suitable  representation  may  be  taken  in  the  form  of 
the  exponential 


f^Sg,  §)  =  A  S“  ekiS 

(217) 

,  -  0  k2E 

01(Ee ,  E)  =  B  Eg  e  1 

(218) 

-  k 

(219) 

=  K(Ee)  E  =  0  E 

1  +  >  Ee 

where  a,  B,  ...,  kQ,  k1  and  kj  are  determined  experimentally.  For  example, 
assuming  linear  behavior  at  sufficiently  small  strains,  then  B/A  =  2G  (i.e., 
twice  the  instantaneous  shear  modulus).  The  functions  R(t)  and  V(t)  represent 
the  creep  and  relaxation  functions  in  shear. 

For  a  creep  test  under  constant  stress  SQ,  then  equation  (207)  gives 


3GI  - 


(1  +v) 


1  .e_  1  .5  _k2E0  (t)  (J-2.)  .  J  k]so/j 


Eo(0 


=  S, 


0(0 


(220) 


when  v  =  constant,  and  for  two  arbitrary  stresses  SQ  and  S 


or 


®  1  +  & 


(1  -  2v)  (En(t)  -  En(t)) 


(22  2) 


(1_L2  \lo9e(f°  ).kl'S°'1 
'  l  -  a  f  e'ci  ’  3d  +  B 


Now,  as  noted  previously,  the  ratio  (1  +a)/(1  + /J )  may  be  d<  ter mined  from  tests 
in  pure  shear.  Then,  from  equation  (120),  k;>(1  +  2r)/(1  +  id  and  k -| / ( 1  +  0)  may 
be  determined  from  a  log-log  plot  of  E0(t)/E0(t)  versus  Sq/S1 0  for  several 
values  of  creep  stress  and  strain.  Finally,  the  function  may  be  evaluated 

from  equation  (221). 

The  constants  kQ  and  1  should  be  determined  from  pressurized  tests  for  vari¬ 
ous  effective  strain  levels,  E  . 


4.6.6  Non-Isothermal  Character: zation 

For  the  situation  of  transien  mperature  loadings,  e.g.,  simultaneous 
cooling  and  straining,  the  uniaxial  o'  .stitutive  equation  was  written 


■I 


a[e(T)  -  aAT  ( r)] 

E(£  -  l') - - - -  dr 


(223) 


-  '/* 

J  T 


(u) ,  e (u) ,  e(u; 


(224) 


The  "damage"  shift  function  was  determined  from  uniaxial  tests,  and  Aj  was 
determined  from  uniaxial  tests  at  different  temperatures.  The  possibility  of 
the  requirement  of  a  thermomechanical  coupling  coefficient  is  permitted  through 
the  shift  function,  Af.  In  practice,  Af  and  Atj  can  be  determined  as  a  single 
coefficient  from  simultaneous  cooling  and  straining  tests. 


A  determination  of  Af  was  not  made,  sir.ro  the  computer  code  used  to  make 
predictions  was  not  modified  to  preform  tho  integration  of  equation  (22?)  for 

transient  temperature  loadings. 


4.7  THE  SWANSON  NONLINEAR  CONSTITUTIVE  L.*W 
4.7.1  Original  Model 


The  framework  for  this  theory  was  established  by  t.n 
typical  behavior  aspects  of  high-elongation  propellants  - 
features  considered  were:  (1)  the  usual  viscoelastic  ‘  : 
on  the  strain  rate,  (2)  the  ability  of  the  solid  propel- : 
strains,  (3)  the  marked  deviation  of  the  solid-propeliar - 
associated  with  Linear  Viscoelasticity,  as  evidenced  by 
exhibited  under  cyclic  loading  of  many  solid  propellants , 
and  (4)  the  dependence  of  the  stress-strain  response  :<r. 

Although  it  is  not  essential  to  have  done  ’  t  t  l'  v 
handling  large  strains  was  incorporated  into  the  cons: 
the  cauchy-stress  tensor  ( O )  as  a  measure  of  the  state 
conjugate,  the  left  Cauchy-Green  deformation  tensor  (3}  ■ 
of  straining.  The  Cauchy  stresses,  defined  in  terms  of 
area,  are  also  called  "true"  stresses.  In  a  principal  c. 
on  the  diagonal  form: 


B  = 


X?  0  0 
0  A§  0 
.  0  0  *3  . 


in  which  the  X^’s  are  simply  the  extension  ratios  in  tn* 

The  remaining  aspects  of  the  observed  response  of 
eled  through  the  use  of  a  softening  function  as  a  stvv- 
major  constitutive  assumption  in  this  theory  relates  .. 
the  deviatoric  stress  and  deformation  tensors  through 


VTv  =  (f)  (g) 
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"C  veto  account  some 
The  principal 
.  .•?  of  the  response 

ii.:  to  sustain  large 
"vrponse  from  that 
V  large  hysteresis 
,  even  at  small  strains, 
..-•cr imposed  pressure. 

the  capability  of 
.•,vr  equations  by  using 

1  i Tress  at  a  point.  Its 
/.as  used  as  the  measure 
force  per  unit  deformed 
o ordinate  system  B  takes 

(225) 


-i'cipal  directions. 

- ..  propellant  were  mod- 
' or r cation  factor.  The 
invariants  of 

(226) 


This  separable  form  has  been  used  previously ^ > 29  and  is  motivated  by  the 
fact  that  the  constant  strain  rate  tensile  curves  are  roughly  similar. 


In  equation  (224),  f  is  the  following  viscoelastic  function: 


f  = 


(22  7) 


with  G  being  the  relaxation  modulus  in  shear,  taken  in  this  theory  as  one  third 

the  tensile  relaxation  modulus,  and 


g  =  softening  function 

aij  s  °ij  “  (akk^3)  6  ij  =  deviatoric  stress  tensor 
B'ij  =  Bij  -  (Bkk/3)5ij  =  de-iatoric  deformation  tensor 


II  a 


[a  11  a 22  +  <*22  a 33  +  a33 


Ml  l 


,v  2 
a12 


:23 


L31 


(228) 
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Second  invariant  of  tensor  a  =  o,  B 

Now,  g  is  a  function  of  deformation  and  pressure  (mean  stress)  and  can  be 
considered  to  be  primarily  a  strain-softening  function.  It  is  defined  as  that 
function  of  the  invariant  VUg’  that  will  force  the  viscoelastic  Cauchy  stress  to 
coincide  with  the  experimental  results;  thus,  unloading  hystersis  as  well  as  the 
effects  of  pressure  may  be  readily  incorporated  into  this  theory,  simply  by 
obtaining  the  corresponding  forms  of  the  softening  function  under  such 
conditions. 

The  softening  function  corresponding  to  virgin  loading  is  obtained  by 
fitting  the  model  to  uniaxial  tensile  tests  at  constant  crosshead  speed.  Under 
these  conditions,  the  deviatoric  stress  invariant  reduces  to 


where,  again  i  is  Cauchy  stress. 


Assuming  incompressibility 


M  ^2  x3  = 


(230) 


and  noting  that 


X  2  =  X  3 


(23D 


the  deformation  invariant  becomes 


(232) 


Taking  the  rate  of  change  of  the  invariant  as  being  approximately  constant 
results  in 


G  (t  -  t)  dr 


(233) 


so  that,  from  equations  (227),  (230),  and  (232)  the  following  is  obtained: 


—  -  gVlx*  f 

VT  J 


G  (t  -  t)  d  r 


(234) 
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from  which  the  softening  function,  g,  may  be  obtained.  The  assumption  leading 
to  equation  (232),  that  the  time-rate  of  change  of -^Ilg'  is  approximately 
constant,  need  be  guarded  against  conditions  of  changing  strain  rate.  For 
example,  as  in  dual-rate  tests,  where  viscoelasticity  does  not  predict  as  fast 
a  response  to  the  rate  change  as  is  experimentally  observed. 

The  modification  to  linear  viscoelasticity  necessary  to  accommodate  this 
behavior  is  as  follows.  The  response  of  the  function  f  in  jouation  (235)  to  a 
constant  time  rate  of  change  of  the  deformation  invariant  is  defined  as  fc.  It 

can  be  expressed  as 


The  modification  to  the  f  function  ic  one  in  an  incremental  manner  through 


and  the  following  incremental  relationship  is  used: 


The  parameter  /3  governs  the  response  of  the  f  function  under  changing  strain 
rates.  As  0>O,  the  response  is  analogous  to  linear  viscoelasticity. 


The  algorithm  developed  by  Herrman  and  Peterson^O  has  been  used  to  impiement 
the  calculation  of  the  convolution  integral  for  f.  In  brief,  let  the  shear 
relaxation  modulus  be  represented  by  a  Prony  series  as 


which  is  directly  analogous  to  linear  viscoelasticity.  The  modification 

proposed  above  can  then  be  implemented  as 


(^n,i) modified  =  ^n,i  +  ^[*cn,i  “  ^n,ij 


(244) 


and  the  I  terms  can  be  incremented  according  to 


zn , i  =  In-1,i  +  (Aln, ^modified 


(245) 


This  has  the  effect  of  changing  each  term  of  the  series  so  that  it  approaches 
the  value  it  would  have  been  if  it  was  always  at  the  new  strain  rate.  Note 
again  that  (as  discussed  in  Reference  30)  varying  temperatures  can  be 
incorporated  into  the  time  scale  as  usual. 


Unloading  tests  required  a  fu  th°r  refinement  of  the  model.  For  lack  of 
more  detailed  information,  the  paramet  r  0  may  be  taken  as  zero  for  unloading 
states  (i.e.f  states  in  which  ^  Hg'  is  decreasing).  The  large  amount  of  hys¬ 
teresis  seen  in  load-unload  cycles  is  tnen  modeled  in  part  by  the  hysteresis 
inherent  within  linear  viscoelasucity  primarily  through  the  g  function.  This 
is  accomplished  by  giving  g  a  different  value  when  the  deformation  invariant 

is  less  than  its  maximum  previously  achieved  during  the  loading  history. 
If  yiIB»max  is  the  current  maximum  value,  the  function 


I'WlIlU  >  '  -  c1  1'  -v«  B'  W%BI  1 


provides  plasticity-like  behavior. 


The  behavior  of  the  g  function  for  unloading  and  reloading  conditions  is 
illustrated  in  Figure  135  (taken  from  Reference  9). 

The  Swanson  approach^  was  used  with  an  only  limited  degree  of  success  to 
predict  the  stress  response  of  TP-H1011  and  UTP-19,360B  under  several  strain 

histories. 
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Figure  135. 


Effect  of  Deformation  and  Pressure  on  the  Strain 
Softening  Function 
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In  the  case  of  TP-H1011,  the  errors  in  the  predictions  were  believed  to  be 
due  to  uncertainties  in  the  value  of  the  changing-rate  coefficient  ((3).  There 
was  no  data  available  to  determined  0  directly  for  this  propellant. 


It  was  possible  to  characterize  UTP-19,360B  in  a  complete  fashion.  The 
corresponding  predictions  were  not  any  better  than  those  obtained  for  TP-H1Q11. 
This  led  to  changing  the  law  as  discussed  below. 


4.7.2  Current  Model 

Analysis  of  the  stress  predictions,  carried  out  for  UTP-19,360B  with  the 
original  Swanson  theory,  revealed  the  importance  of  several  inadequacies  and 
oversimplifications  listed  below. 


1.  The  softening  function  (g)  should  depend  not  only  on  the  strain  and 
pressure  but  also  on  the  strain  rate. 

2.  The  softening  function,  as  defined  by  equations  (227)  and  (235),  should 
be  different  for  unloading  than  for  reloading. 

3.  The  softening  function  for  unloading  or  for  reloading  should  never 
become  zero  for  conditions  of  tensile  straining  only.  A  zero  value 
could  occur  with  the  softening  function  defined  by  equation  (246). 
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4.  The  healing  process  observed  during  relaxation  in  solid  propellants 
like  UTP-19i360B  was  not  taken  into  consideration  by  the  original 
Swanson  theory. 

5*  The  reverse-recovery  observed  in  solid  propellants  during  relaxation  or 
rest  periods  that  follow  an  unloading  process,  only  poorly  modeled  by 
classical  viscoelasticity,  is  not  considered  in  the  approach  by 

Swanson . 

6.  The  changing-rate  coefficient  ( /3)  is  more  a  mathematical  device  than  it 
is  a  material  property.  If  the  softening  function  is  made  to  depend  on 
the  strain  rate  then  0  need  not  be  used. 

7.  The  use  of  a  softening  function  as  a  stress  correction  factor 
eliminates  the  need  of  using  the  Cauchy  stress  ( O )  and  the  nonlinear 
measure  of  stretching  (B). 

All  these  observations  were  incorpcr  led  into  the  original  stress-strain 
law  but  the  general  form  of  the  corresponding  equations  remained  the  same, 

namely 


n  =  ^3~  (g)  J  a  (st  -  sr)dVn^T  dr 

n  a  t 


(247) 


valid  for  one-dimensional  loading,  with: 


(248) 

t 


representing  temperature-t educed  time;  and  where  the  time- temperature  shift 
function  was  taken  in  the  power-law  form 


At 


Tr  -  Ta  V  ® 
T  -  T  a  / 


(249) 
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in  which  Tp  is  the  shift  reference  temperature,  Ta  and  m  are  material 
parameters,  and  T  is  the  current  temperature. 

The  modified  version  of  the  Swanson  theory  was  most  successfully  used  to 
predict  the  response  of  UTP-19 ,360B,  as  explained  next. 

4.7.3  Stress  Predictions 

The  degree  of  accuracy  of  the  predictions  made  with  the  current  version  of 
the  Swanson  approach  may  be  realized  by  examining  Figures  136  through  148.  The 
first  two  figures  correspond  to  the  lowest  and  highest  constant-rate  tests 
available.  Figures  138  and  139  present  the  results  for  the  dual-rate  tests 
while  Figure  140  pertains  to  the  saw-tooth  test  at  constant  rate  and  increasing 
peak  strains.  Figures  141  to  145  show  the  predictions  corresponding  to  complex 
multiple  loading,  twenty-four  hour  relaxation,  long-duration  similitude,  three- 
step  relaxation,  and  predamage  relaxation.  Finally,  figures  146  to  148  show  the 
results  obtained  for  comparable  constant  rate  tests  at  70  F,  123  F,  and  40  F. 

They  testify  to  the  fact  that  the  time-temperature  superimposition  prin¬ 
cipal  may  be  used  without  sacrificing  more  accuracy  than  is  already  lost  in 
fitting  equation  (185)  to  the  very  limited  time-temperature  shift  data. 

4.7.4  Material  Characterization 

According  to  this  theory,  only  the  following  listed  properties  are  needed 
to  characterize  a  solid  propellant  completely: 

1.  The  relaxation  function,  G,  as  defined: 

Erei  (t)  (250) 

G  (t)  =  - — 

3 

where  Erei(t)  is  the  linear  viscoelastic  relaxation  modulus. 


(Text  continued  on  pg.  263) 
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Nonlinear  Viscoelastic  Stress  Predictions  for  Constant-Rate  Test 


Nonlinear  Viscoelastic  Stress  Predictions  for  Two-Rate  Test 
(UTP-19,360B-lJ00/1777) 


Figure  139.  Nonlinear  Viscoelastic  Stress  Predictions  for  Two-Rate  Test 

(UTP-19 , 360B-400/1777 ) 


iRure  140.  Nonlinear  Viscoelastic  Stress  Prediction  for  the  Saw-Tooth  Test 


Figure  141.  Swanson  Theory  Stress  Predictions  for  the  UTP-19,360B-400/1777  Complex  Multiple  Loading 

(Test  No.  10) 


ieure  1H?.  Swanson  Theory  Stress  Predictions  for  U7P-  i  ,  ~!l'  top 

(Test  No.  8) 


Swanson  Theory  Stress  Predictions  for  UTP-19,360B-^00/1777  Predamaged  Relaxation 


Figure  147.  Swanson  Theory  Stress  Predictions  for  UTP-19,360B-400/1777  Constant  Strain  Rate 


Swanson  Theory  Stress  Predictions  for  UTP-19 , 360B-^00/1777  Constant  Strain  Rate 


The  noft.on  i  rip  function,  g ,  defined  end  obtained  an  follows: 

a.  For  loading  conditions: 

gL  =  g  (e,  e)  (25D 

and  it  is  obtained  from  a  sequence  of  constant  rate  tests  with 
at  least  three  different  rates  that  span  the  range  expected  in  the 
applications. 

b.  For  unloading  conditions: 

Su  =  S  (f/fmax)  (252) 

where  emax  represents  the  maximum  strain  previously  achieved  dur¬ 
ing  the  loading  history.  The  g  is  determined  from  the  unloading 
portion  of  a  loading-unloading  cycle  carried  up  to  an  intermediate 
strain  level. 

c.  For  relaxation  conditions: 

gr  =  gr  (t  -  tD)  (253) 

in  which  tQ  is  the  time  at  which  the  relaxation  process  begins 
and  g  is  evaluated  from  a  relaxation  test  at  an  intermediate 
strain  level. 

In  addition,  for  relaxation  after  partial  unloading  or  during  rest  periods 
starting  at  t  =  tQ: 

gr  =  - - -  (25*) 

gr  (t  -  tQ) 

Also,  the  stress-correction  function  for  reloading  is  taken  as  a  linear 
function  of  the  relative  strain.  It  is  a  straight  line  from  the  point  where 
reloading  starts  to  the  point  of  maximum  loading  over  the  past  history. 


4.7.5  Three-Dimensional  Version  of  the  Model 

The  (general)  constitutive  assumption  used  to  relate  the  deviatoric  compo¬ 
nents  of  the  stress  and  deformation  tensors,  takes  the  following  form: 


a 


ij 


2,  3 


(255) 


together  with 


Vila''  =  (g)  (f) 


(256) 


or,  equivalently 


r  >  a  ^IIB*  d  t 
a  t 


(257) 


where 

=  i-J  component  of  the  deviatoric  Cauchy  stress  tensor 
Bj'j  =  i-j  component  of  the  deviatoric  Left  Cauchy-Green  deformation  tensor 
II0 ,  IIg,  =  second  invariants  of  the  deviatoric  stress  and  deformation  tensors 


with 


def 
na'  ■ 


°ij  =  °ij  (011  +  °22  +  °33)6iJ;  i,  j  =  1,  2,  3 
”[®11  °22  +  011  °33  +  °22  °33^  +  (°12)2  +  (°13^2  +  ^°23^2 


(258) 


(259) 
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and  similarly  for  and  Ilg';  and  in  which 


6 


ij 


def 


1  for  i  =  j 
0  for  i  =*  j 


(260) 


also 


g  *  softening  function  that  depends  primarily  on  the  strain  level,  the  strain 
rate,  and  the  applied  pressure. 

and 


G(t) 


def 


E(t)/3 


(261) 


where  E(t)  represents  the  tensile  relaxation  modulus  at  a  small  strain. 

According  to  the  constitutive  assumptions  (255)  and  (257),  the  distortional 
behavior  of  the  material  is  completely  characterized  through  the  softening 
function,  g,  and  the  relaxation  function,  G,  which  may  be  evaluated  from 
one-dimensional  tests,  as  explained  in  the  previous  section.  Indeed,  the 
stress-strain  relations  set  forth  in  equation  (255)  and  (257)  reduce,  as 
they  should,  to  those  employed  in  the  one-dimensional  version  of  the  model. 

To  complete  the  theory,  an  assumption  is  still  needed  about  volumetric 
behavior;  and  although  time-dependent  bulk  response  may  be  important  in  some 
applications,  the  elastic  relation 


~  (a11  +  ° 22  +  °33)  =  K  0^1  X  2  A  3  [  -  1  ) 


may  be  employed;  in  which  K  is  the  bulk  modulus  and  theX^’s  are  the  stretch 

ratios. 

For  an  incompressible  material  (and  solid  propellants  are  nearly 
incompressible) 


>1X2X3  =  1 


(263) 


so  that  equation  (262)  breaks  down,  and  the  stress  tensor  has  to  be  considered 
a  function  of  the  mean  pressure  (  ^  +  a22  +  °3 3 ) / 3 ,  as  well  as  of  the  deformation 

tensor,  leading,  eventually,  to  a  stress-strain  law  of  the  form  given  in 
equation  (255). 

*1.7.6  Application  of  the  Model  to  Two-Dimensional  Problems 

In  order  to  use  the  stress-strain  law  presented  in  the  foregoing  section, 
one  must  have  available  the  deformation  tensor  at  each  point  of  the  continuum 
where  the  stresses  are  desired.  rhis  solution  in  terms  of  deformation  may  be 
arrived  at  numerically  or  analytically  through  finite  elements. 

The  accuracy  with  which  the  present  constitutive  theory  may  predict  the 
two-dimensional  response  of  solid  propellants  may  be  seen  in  Figures  1  49  to  154, 
which  correspond  to  constant  strain-rate  tests  of  strip-biaxial  samples  of  UTP- 
19.360B.  The  first  three  figures  belong  to  tests  performed  at  a  nominal 
crosshead  displacement  rate  of  0.02  in. /min  at  40  F,  70  F,  and  120  F,  respec¬ 
tively;  while  Figures  150  to  154  show  the  results  for  a  crosshead  displacement 
rate  of  0.2  in. /min  at  the  same  low,  intermediate,  and  high  temperatures  of  40, 

70,  and  120  F.  The  plotted  data  refer  to  the  direction  of  applied  loading,  which 
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is  also  the  direction  of  maximum  principal  stress  and  strain.  The  geometry  of 
the  strip-biaxial  sample  used  is  as  presented  in  the  following  sketch: 


a/b  =  1.25/6  2  4.8 

o  22//°11  h  0  =  0.481 

a 

e  22^e  11s  0C  =  0.025 


The  stress-  and  strain-axiality  factors,  0O  and  0e  were  taken  from  Refer¬ 
ence  31,  and  are  valid  at  the  center  of  the  sample  for  small  strains  only. 


The  constitutive  relations  given  in  (255)  and  (257)  yield: 


o 


aVHjT 

3  T 


dr 


(265) 


where 


aYnjj’.  dViigi  3B^ 

9  t  '  3  at  ’ 


i,  J  = 


2,  3 


(266) 


with  summation  implied  over  repeated  indices. 

Now,  under  conditions  of  plane  stress  of  an  incompressible  material,  and 
along  the  principal  directions,  one  has 
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vv 
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(2o,  _o2)  o 

0  (20  2  -  oi) 

0  0 
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0 

(o-,  +  02 ) 


IIB’  =  — B  -j  -j  B22 


•  1  ft 

BH  B33  -  B22  B33 


To  evaluate  (266),  we  first  write  it  in  unabridged  notation,  noting 
case,  if  i^j  then  B^j  =  0;  thus: 


aVTiB-*  3VnB,  ae^  aVTf^TdB^  aVnB,  3b33 

- =  - —  -  +  - . - + - ; -  - 

3  t  3  Bi  1  3  t  3  B22  3 1  3  B33  3  t 


and  using  (269): 


3Viib  .  1 

3  t  ”  2Vi  lB^ 


B33) 


3Blj  + 
3t 


(-B33 


B11 ) 


3B22 
3  t 


+ 


,  »  f  3b33 

+  (-Bn  -  B22)  — — 

0  t 


where,  from  (268)  one  has,  for  instance,  that 


3Bn 

at 


dA1  A  dA2  X  dA3  \ 
dt  ^  dt  J  dt  / 


and  similarly  for  the  derivatives  of  B22  and  B33. 


(267) 


(268) 


(269) 


that  in  this 


(270) 


(271) 


(272) 


27*4 


h<:  previous  derivations,  the  stretch  ratios  are  computed  as 


*  1  =  1  +  f 1  (t) 

^  2  =  1  +  c2  (t)  =  1  +  <t> '£«i  (t) 


(273) 


^3  =  1  /  ( X  1^-2 ) 

in  which  c  -,  ( t )  is  the  strain  history  imposed  on  the  sample  along  coordinate  1, 
and  the  last  expression  of  (273)  follows  from  the  incompressibility  condition, 
equation  (263). 


Hence,  the  first  component-equation  of  the  constitutive  relation  (265) 


yields : 


°11  5 


~  (2CT  1  “  ff2) 


>iiL  f 

W  J 
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or,  in  view  of  (264) 


(2  -  <t>0) 


sL  f 

;iB’  J 


o  -~dr 


and  finally: 


(2  ~  *o) 


L  f 

■B'  J 


3  ^  1 1 Q  i 

r )  - dr 

d  t 


(274) 


Using  equations  (264),  (268),  and  (271)  to  (274),  we  obtained  the  response  of 
the  biaxial  sample  in  the  direction  of  the  applied  loading.  The  plots  included 
in  this  report,  3how  the  engineering  stress: 


O  ^  =  0  1  / 


(275) 


rather  than  the  Cauchy  stress,  o  . 


l».7.7  Application  to  Transient  Temperature  Loadings 

The  model  in  the  previous  section  has  been  modified  to  the  form 


Viiol  -  Ap  g(B.B) Erel  fe-  f>  i  V>v  dT 


(276) 


where 


(277) 


and  the  thermomechanical  coupling  coefficient,  Af,  has  been  introduced  for  the 
transient  straining  and  cooling  loading  history.  The  deformation  gradient,  F, 
and  the  Green  strain  tensor,  B,  incorporate  the  free  thermal  expansion,  a  AT. 

Figure  155  shows  that  Ap,  as  determined  from  uniaxial  straining  and  cooling 
tests,  is  nearly  constant  (Af  =  1.27)  for  this  approach.  Figures  156  through 


Figure  155.  Swanson’s  Linear  Viscoelastic  Theory 

3095? 
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1  ( /  how  the  agreement  with  uniaxial  and  biaxial  straining  and  cooling  tests. 
The  biaxial  ramp-relaxation-ramp  test  was  particularly  well  predicted  by  this 
approach;  much  better  than  by  any  other  approach.  The  large  disparity  shown  in 
Figure  160,  in  contrast  to  the  good  agreement  in  Figure  159,  lends  credence  to 
our  previous  contention  in  section  *1.2.3  that  this  is  simply  a  bad  test  result. 
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5.0  SUBSCALE  MOTOR  ANALYSES  AND  STRESS  PREDICTIONS 


5.1  LINEAR  ELASTIC  ANALYSES 

Linear  elastic  stress  analyses,  following  commonly  accepted  handbook  and 
industry  practices,  were  performed  to  provide  a  basis  for  comparison  with 
the  nonlinear  constitutive  theory  predictions  and  to  illustrate  the  inadequacy 
of  elasticity  analyses  for  thermal  cycling  of  solid  rocket  motors. 

The  linear  elastic  analyses  consisted  of: 

1.  Handbook  Calculations/Design  Curves 

2.  T EXGAP  2D  Finite  Element  Analyses. 

5.1.1  Handbook  Calculations 

The  solid  propellant  structural  integrity  handbooks ^ *35  SUgge3t  the 
use  of  the  equations 

0rr<b)  .  {X2-  l)PrEpORAT  (278) 

or 


for  the  midplane  case/grain  interfacial  pressure  (i.e.,  radial  bond  stress), 
where: 

A  =  2b/2a  =  grain  OD/grain  ID 

F.p  =  propellant  modulus 
F,c  =  case  modulus 
®R  =  ®p  -  2/3  (1  +  »’c)  ac 

ap  =  propellant  coefficient  of  linear  thermal  expansion 


ac  =  case  coefficient  of  linear  thermal  expansion 
h  =  case  thickness 
Sc  =  “c  At 

5P  =  3s  +  «pAT  =  ap  { (T0  -  Tc)  +  AtJ 
T0  =  stress  free  temperature 
Tc  =  motor  cure  temperature 
Pr  =  finite  length  correction  factor. 

Equations  (278)  and  (279)  are  derived  from  the  Lame  equations  for  a  composite 
cylinder  assuming  a  thin  elastic  casing  and  an  incompressible  inner  core.  In 
equation  (278),  AT  is  referenced  to  the  stress-free  temperature,  and  in  equation 
(279),  AT  is  referenced  to  the  motor  cure  temperature.  For  the  situation  where 
bEp/hEc«1,  equation  (279)  reduces  to  equation  (278).  Hence,  equation  (278) 
was  used  to  make  elasticity  stress  predictions  for  the  subscale  motors. 

Figures  16 1  and  162  compare  elastic  stress  predictions  using  equation 
(278)  with  measured  radial  stresses  for  the  0.75-in.  subscale  motor  with 
072-19,3608-400/1777  propellant.  In  unsse  calculations: 

*=  5.37 

«p  =  5.3  x  10-5  in./in./p 
ac  =  6  x  10“5  in./in./F 
Ep  =  151 .4  psi 
Pr  =  0.83 
•'c  =  0-3 

Figure  161  shows  stresses  calculated  from  the  measured  strain-free  tempera¬ 
ture  of  166.5  F,  whereas  Figure  162  is  based  on  the  measured  stress-free  tempera¬ 
ture  of  14 1  F.  Not  all  the  points  are  connected  in  these  figures. 

As  expected,  the  linear  elastic  predictions  underpredict  stresses  during 
transient  conditions,  but  compare  favorably  with  the  longer  term  equilibrium 
behavior  whan  the  measured  stress-free  temperature  is  used  in  the  calculations. 
The  experimental  results  obtained  here,  as  suggested  in  Reference  3,  indicate 
that  it  is  incorrect  to  assume  stress  and  strain-free  temperatures  as  identical, 
even  at  the  completion  of  cure.  Since  the  stiffness  of  the  motor  case  is  so 


284 


lant  Temperature 


O  <j*j]  i 

^  <£)  E>|  ™ 

cxfdto  i 

O  <10 L  ^  I 

<350! 


IB 

14 

12 

10 

8 

6 

4 

55 

50 

45 

40 

35 

30 

25 

20 

15 

10 

5 

0 

1-3  x  (j) 

3HniVU3dW31 

0 

3  x 

(isd) 

SS3tiiS 

Fiaure  161.  Elastic  Stress  Predictions  for  UTP-19 ,360B-400/1777 
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■  tf^n  that  of  the  propellant,  the  state  of  deformation  in  a  propel- 

.  largely  independent  of  the  propellant  modulus.  Thus,  the  strain- 

■  .■  :  ."porature  is  determined  by  the  original  geometry.  In  contrast,  stresses 
arc  proportional  to  the  propellant  modulus,  and  there  is  no  a  priori  justifica¬ 
tion  I  or  a.i-  umihg  that  the  stress-free  temperature  and  strain- free  temperature 
arc  the  same.  The  difference  is  significant,  and  care  should  be  exercised  in 
making  motor  stress  and  strain  predictions. 

It  should  be  noted  that  the  stress  and  strain-free  temperatures  can 
generally  only  be  defined  for  one  components  of  stress  or  strain.  It  is 
probably  not  possible  to  make  all  component  of  stress  or  strain  vanish  at  the 
same  temperature.  For  example,  the  normal  procedure  for  determining  the  strain- 
free  temperature  in  a  circular  port  motor  for  example,  is  to  make  bore  measure¬ 
ments  at  the  motor  midplane  at  several  temperatures  near,  below,  and  possibly 
above  the  cure  temperature,  and  then  extrapolating  to  the  temperature  where  the 
bore  strain  goes  to  zero.  One  generally  find3  that  the  entire  mandrel 
dimensions  are  not  obtained  at  the  so-called  strain-free  temperature.  The 
reasons  for  this  are  severalfold.  For  one  thing,  propellants  are  not  truly 
incompressible  even  though  Poisson's  ratio  is  very  nearly  equal  to  one-half. 

Second,  it  is  not  always  possible  to  eliminate  or  compensate  for  slump.  Also, 
the  propellant  is  cured  or  polymerized  under  varying  degrees  of  constraint, 
depending  upon  the  location  within  the  motor. 

With  regard  to  the  stress-free  temperature,  it  is  normally  obtained  in  a 
manner  similar  to  that  for  the  strain-free  temperature  determination  only  using 
stress  transducer  data  for  the  radial  bond  stress  at  the  midplane  of  the  motor. 

Again  only  one  component  of  stress  is  used,  and  the  reasons  for  all  components 
of  stress  not  vanishing  are  the  same  as  those  given  in  the  previous  paragraph 
with  regard  to  strain. 

5.1.?  TEXGAP  Finite  Element  Analyses 

Linear  elastic  axisymmetric  finite  element  analyses  using  TEXGAP  2D  were  also 
performed  on  the  subscale  motors  subjected  to  thermal  testing.  The  geometry 
analyzed  is  shown  in  Figure  163,  and  the  undeformed  grid  is  shown  in  Figure  1 64 . 
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Propellant 


Figure  1 63 -  Circular  Port  Motor  Geometry 

The  analyses  considered  a  thermal  load  of  100  F  and  used  the 
properties: 


13924 

fol lowing  material 


•'propellant  s 

“propellant  =  5  x  10  ^  in. /in., 
^propellant  =  151 *4  psl 
•'case  =  °*31 

“case  =  ^  x  10-5  in./in./F 
Ecase  =  3  x  10?  pai 


The  radial  component  of  the  calculated  mid  plane  bond  stress  is  also  shown 
Figures  161  and  162.  As  expected,  these  stresses  agree  exactly  with  the 
handbook  calculations. 


5.2  APPROXIMATE  VISCOELASTIC  ANALYSES 

-O' 

Thermal  stress  analyses  of  solid  propellant  rocket  motors  are  routinely  con 
ducted  using  an  approximate  (quasi-elastic)  viscoelastic  analysis.  The  pro¬ 
cedure  involves  the  use  of  an  equivalent  propellant  modulus  which  approximates 
the  true  time  and  temperature  behavior  in  conjunction  with  an  elastic  stress 
analysis.  The  immediate  objectives  of  the  motor  analyses  are  to  determine  the 
sophistication  and  the  extent  of  realism  in  the  modulus  approximation . 


I 
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Figure  164.  TEXGAP  Finite  Element  Grid  for  Circular  Port  Subscale 

Motor  Analysis 


The  effective  modulus  is  determined  assuming  thermorheologically  simple  mate¬ 
rial  (TSM)  behavior  or  thermorheologically  complex  material  (TCM)  behavior.  The 
decision  to  perform  a  TSM  or  TCM  analysis  is  based  on  comparisons  of  predicted 
linear  thermoviscoelastic  stress  response  with  the  measured  stress  response 
during  simultaneous  straining  and  cooling  tests  of  laboratory  specimens  con¬ 
ducted  at  strain  rates,  strain  levels,  and  over  temperature  ranges  representa¬ 
tive  of  the  anticipated  motor  environment.  If  the  linear  predictions  agree  with 
the  measured  stresses,  then  a  TSM  analysis  is  conducted;  if  t~.ey  do  not  agree, 
then  a  thermochemical  coupling  coefficient  is  introduced  ana  a  TCM  analysis  is 
conducted . 

The  stress,  c^j  (Xk,T),  at  any  point  xk  (xk  =  x.y.z,  or  r,  0,  '  in  he 

propellant  grain  at  time  t  are  calculated  by  the  quasi-elastic  method  ur  t,.e 
formula 


(Xk,  t)  =  Kij  (X.,) 


Eeff  (t)  “p  A i a ( t ) 
Efem  “fematfem 


(2501 


while  the  strains,  Cjj  (xk,t)  are  calculated  from 


*1j  (xk>  t)  *  Kij  (Xk)  a 


“p  ATef t) 


(281) 


FEM  atFEM 


The  subscripts  FEM  in  equations  (280)  and  (281)  denote  the  values  of  modulus,  E, 
linear  coefficient  of  thermal  expansion  (LCTE),  a,  and  temperature  increment 
from  the  stress  or  strain- free  temperature,  AT,  used  in  an  elastic  finite  ele¬ 
ment  analysis.  The  constants  Kij(xk)  and  K’jj  (xk)  are,  respectively, 
the  stresses  and  strains  at  the  point  xk  obtained  from  elastic  finite  element 
analysis  using  Ep^,  “pgw  30(1  ATFEM  33  inPut*  The  constants  jq  j  and  K'^ 
incorporate  the  effects  of  motor  geometry  and  the  material  properties  of 
the  case,  insulation,  and  liner.  They  also  depend  or.  the  ratio  of  cane  u 
propellant  modulus.  These  latter  properties,  if  taken  to  be  nn  average  of 
the  properties  over  the  temperature  range  of  interest,  introduce  only  r.ma  11 
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r-'r-s  in  equations  (280)  and  (28 1 )  due  '  o  their  own  viscoelastic  behavior. 
Also,  to  ensure  the  validity  of  equations  (.280)  and  ( 28 1 ) ,  Poisson's  ratio 
of  the  propellar'  ,  vn,  must,  be  time-independent;  that  is,  the  bulk  response 
of  the  propellant  must,  be  nea.  ly  elastic,  which  is  the  typical  assumption 
made  in  stress  analysis  of  propellant  grains.  Equations  (280)  and  ( 28 1 ) 
also  assume  that  the  propellant  LCTE  is  a  constant;  however,  they  are  equally 
valid  if  «p  is  a  function  of  temperature. 

Four  stress  predictions  were  made  for  the  0.75-in.  bore  structural  analog 
involving  four  different  approximations  to  Eeff.  The  first  involved  simply 
replacing  the  elastic  modulus  by  the  relaxation  modulus  at  the  appropriate 
time  and  temperature;  that  is, 

E~*  Erel  <t/AT> 


A  comparison  with  measured  stresses  is  shown  in  Figures  165  and  166  without 
all  points  connected.  As  in  the  case  of  the  elastic  analyses,  the  results  agree 
well  with  the  long  term  relaxation  behavior,  but  are  low  by  a  factor  of  2  to  3 
during  transient  temperatures. 

The  next  set  of  predictions  were  made  using  linear  viscoelasticity  theory: 


■if 


AE({  -  )  dr 


(282) 


■f 


dx 


At  LT(x)J 


(283) 


Results  for  this  prediction  are  shown  in  Figures  1 67  through  170.  Again, 
the  agreement  is  not  too  discouraging  during  relaxation  periods,  and  the 
last  temperature  excursion  but  are  low  by  a  factor  of  two  during  the  early  tran¬ 
sient  thermal  conditions.  The  agreement  during  the  last  cycling  may  be  the 
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Figure  166.  Viscoelastic  Stress  Predictions  for  UTP-19,360B-il00/1777 
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Figure  170.  Linear  Viscoelastic  Stress  Predictions  for  UTP-19, 360B-400/1777 

Subscale  Motor  Thermal  Cycling  (C-510) 


result  of  either  one  or  a  combination  of  two  observations.  First,  the  last  tem¬ 
perature  excursion  was  to  approximately  40  F.  The  motor  had  previously  been 
subjected  to  two  temperature  excursions  below  40  F.  Much  industry  work  has 
shown  that  solid  propellants,  when  reloaded  to  a  lower  strain  level  after  being 
subjected  to  a  prior  loading  history,  display  essentially  linear  viscoelastic 
responses.  Hence,  the  propellant  could  be  in  a  state  of  near-constant  damage, 
and  its  behavior  approximates  that  of  a  linear  viscoelastic  material.  On  the 
other  hand,  the  agreement  could  be  the  fortuitous  result  of  a  shift  in  the 
stress-free  temperature  from  141  F  to  a  temperature  closer  to  105  F.  Imme¬ 
diately  prior  to  the  last  temperature  excursion,  the  analogs  had  been  stored  at 
105  F  for  nearly  70  days.  This  time  is  sufficient  to  cause  a  shift  in  the 
stress-free  temperature  to  the  storage  temperature.  The  result  is  that  stresses 
would  be  lower  on  subsequent  temperature  excursions,  but  linear  viscoelastic  pre 
dictions  would  be  incorrect  since  they  are  based  on  too  high  of  a  stress-free 
temperature.  If  one  corrects  this  probable  event,  then  the  viscoelastic  stress 
predictions  are  decreased  by  about  40? ;  which  would  then  bring  the  error  in  line 
with  the  error  for  the  previous  thermal  excursions. 

Several  predictions  were  made  using  the  simple  modification  to  linear 
viscoelasticity  theory  of  incorporating  a  thermomechanical  coupling  coefficient 
(Af)  into  the  equations  for  linea*  viscoelasticity.  The  general  procedure,  as 
described  in  section  4.2.3  assumes  that 

Af  =  Af(t> 

and 


o(t)  =  f  Erel  <«-{•>£  {ftF  *  (28-0 

0 

Under  simplifying  assumptions,  however,  assuming  that  Af  =  constant,  equation 
(284)  can  be  written 


°(t)  s  E  (e-<*Ar)  +  A P  /  AE(£-£')  (e-aAr)  dr 


t  Oflr,  -v 


■  *  *  or:  !  h'-  >. fi  in#  •  n  *.n  1  re;-.u! t::  d  i  u;ed  in  section  4|  .2 . 3  f  or  the 

''.or  -':  jr,.  '-nin  gr  •  j.OCr':>  in. /min  uniaxial  straining:  and  cooling  tests,  Af  varied 
between  2.0  and  A  constant  Af  =  2.56  was  used  for  both  the  transient 

cooling-st ra ■ ring  portion  of  the  loading  as  well  as  the  constant  temperature 
relaxation.  This  led  to  the  results  shown  in  Figure  171.  Although  the  tran¬ 
sient  behavior  is  predicted  reasonably  well,  the  relaxation  behavior  is  over 
predicted  by  50$  or  more.  Thus,  the  attempt  to  use  a  constant  Af  does  not  work 
well  in  this  case. 

5.3  LINEAR  VISCOELASTICITY  PREDICTION  PROCEDURES 

Linear  viscoelastic  stress  predictions  were  made  using  the  model  given  in 
equation  (19)  in  section  -'4.2.1.  The  material  characterization  for  uniaxial  iso¬ 
thermal  test  history  predictions  involved  fitting  the  master  relaxation  and  time- 
temperature  shift  curves.  The  relaxation  modulus  curve  was  fit  to  the  exponen¬ 
tial  series: 

m 

Ut)  =  Gea  +  £  6i  ( 286 ) 

i=l 

where  Gep.  is  the  equilibrium  modulus  and  and  are  curve  fit  constants.  The 
equilibrium  modulus  and  reduced  time-modulus  pairs  were  picked  from  the  master 
relaxation  modulus  curve.  The  series  exponents,  a^'s,  were  chosen  as 

a.  =  1/t.  (28? ) 

1  1 

The  constants,  E^  were  taken  as  the  solution  to  the  system  of  simultaneous 
linear  equations 

m 

Ed.)  -  Gan  =  Y'  G.  e'^j  ,-i  =  l,m  (2881 

I  eq  Z— /  1  > 

i  =  l 
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Ihe  time- temperature  shift  function,  A^.,  was  fit  to 


(289) 


where  Tp  was  the  reference  temperature  and  Ta  and  a  were  fit  constants.  From  tem¬ 
perature  AT  pairs  taken  from  the  At  curve,  Ta  and  m  were  determined  using 
Powell's  algorithm  which  minimizes  the  error  function 


FID  =  £ 

i=l 


(290) 


Predictions  for  biaxial  samples  were  done  by  applying  a  correction  to 
the  uniaxial  modulus  or  utilizing  the  biaxial  relaxation  modulus.  For  UTP- 
19.360B  propellant,  the  biaxial  modulus  was  taken  as 


'biaxial 


=  4/3 


'uniaxial 


(291 ) 


where  4/3  is  the  ratio  of  the  biaxial  to  uniaxial  modulus.  UTP-3001  propellant 
biaxial  samples  were  predicted  by  fitting  the  master  biaxial  relaxation  curve  to 
the  exponential  series  described  above.  In  both  cases  the  time- temperature 
shift  function  used  was  the  same  previously  utilized  for  uniaxial  predictions. 

For  predictions  involving  transient  temperature  test  histories,  the  linear 
viscoelastic  stress  was  expressed  as 


o(t) 


(AT) 


V-  G.  T'  1  de  dT 

1  HT 


(29 2) 
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where  Af  is  the  thermo- mechanical  coupling  factor.  Af  was  determined  from 
uniaxial  straining  while  cooling  tests  as 


,  ,.T,  <W(t)  '  G«l  '  (t) 

A  (AT)  =  - 


f  £,  *,  •- 

Jo  1=1 


-«i ft'-  T') 


(293) 


The  AT-Af  function  was  represented  in  a  discreet  point  form,  where  values  for  Af 
were  determined  by  interpolation. 

5.4  MOTOR  STRESS  PREDICTIONS  WITH  SWANSON  CONSTITUTIVE  THEORY 

Equation  (265)  was  applied  to  predicting  the  stress  response  of  the  0.75-in. 
bore  analog.  Since  equation  (265)  is  formulated  in  terms  of  the  second  in¬ 
variant  of  the  deviatoric  stress  and  strain  tensors,  and  since  only  the  radial 
component  of  the  bond  stress  is  measured,  it  is  necessary  to  use  analytical  or 
finite  element  analysis  results  either  to  modify  the  measured  stress  or  to  solve 
the  predicted  deviatoric  stress  for  vne-  radial  component  of  stress.  In  this 
case  we  have  chosen  to  do  the  former. 


The  softening  function  g  alsc  requires  special  attention  when  determining 
it  for  the  Induced  motor  strain  history.  In  this  case,  the  procedure  followed 
was  to  use  a  volume  averaged  deviatoric  strain;  that  is, 
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Since  we  treat  propellant  as  nearly  incompressible,  and  since  the  stress 
state  at  the  case  wall  at  the  midplane  of  a  motor  is  nearly  equal  triaxial 
tension,  the  principal  components  of  the  deviatoric  strain  tensor  (i.e.,  B^', 
Bgg* ,  833')  are  zero.  Ilg'  is  also  zero.  To  apply  the  modified  Swanson's  theory 
correctly  in  this  situation  one  must  make  use  of  a  volumetric  or  bulk  response 
constitutive  relation.  In  this  application,  however,  a  simple  ad  hoc  approach 
was  used  in  an  attempt  to  predict  the  radial  bond  stress  in  the  analog  motors. 


Since  the  stress  state  at  the  bore  surface  is  nearly  a  2:1  biaxial  stress 
state,  lip'  was  evaluated  at  the  bore  surface  and  used  as  the  strain  history 
input.  The  relationship  between  op  at  the  case  wall  and  Ilo'  at  the  bore 
surface  was  then  used  to  calculate  the  predicted  radial  bond  stress. 

Predicted  stresses  are  shown  in  Figures  172  and  173.  While  the  pro¬ 
cedures  have  many  approximations,  the  resulting  agreement  is,  perhaps,  better 
than  should  be  expected. 

5.5  SWANSON  THEORY  PREDICTION  TECHNIQUES 

Swanson's  constitutive  theory  stress  predictions  were  made  using  the  model 
given  in  equation  (17).  Characterization  of  the  material  involved  defining  the 
softening  function  and  viscoelastic  function  parameters. 

The  softening  function,  g(e,  e),  for  virgin  loading  has  the  form: 

g(e,e)  =  A  |>]  B  e"Ce  (295) 


where  A,  B,  and  C  are  fit  constants.  From  constant  strain  rate  uniaxial 
test  data,  the  softening  function  curve  for  individual  strain  rates  were 
determined  as 
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For  each  curve,  A,  B,  and  C  were  determined  using  Powell's  algorithm,  which 
minimizes  the  error  function 
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For  predictions  made  at  intermediate  strain  rates,  values  for  A,  B,  and  C  were 
obtained  by  interpolation. 


The  softening  functions  for  relaxation  and  unloading  were  expressed  as 


g  (e,e)  =  7(At)  *  g(e(tQ),e) 


(298) 


and 


g(e,e)  =  ?(erel)  *  g(e(tQ),  e) 


(299) 


respectively.  g(e(tQ))  is  defined  as  the  softening  function  value  at  the  begin 
ning  of  relaxation  or  unloading.  At,  and  are  given  by 


Ai  :  t  *  t. 
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and 
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The  garana  functions  curves  were  determined  from  load-relaxation  and 
load-unload  test  histories  as 
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The  curves  were  expressed  in  a  discrete  point  manner,  where  values  were 
obtained  by  Interpolation. 


For  relaxation  after  partial  unloading  or  during  rest  periods,  the 
softening  function  was  expressed  as 


g(e,c) 


B(e(t0),  e) 
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The  softening  function  for  reloading  was  taken  as  a  linear  function  of 
relative  strain.  The  relative  strain  is  given  by 


where  emin  is  the  minimum  strain  of  unloading  and  is  the  maximum  strain  of 

the  previous  loading.  The  softening  function  is  expressed  as 


g('-')  *  [’unload  +  [«  -  ’unload]  *'rel]  x  «  <'<V  ^  (305) 

where  Tunioad  is  the  gamma  function  at  the  end  of  unloading  and  0  is  the  ratio  of 
reloaded  to  virgin  stress. 

The  viscoelastic  function  is  given  by  equation  (292).  The  shear  relaxation 
modulus,  G(t),  is  defined  as 


G(t)  =  E(t)/3 


(306) 


where  E(t)  is  the  uniaxial  relaxation  modulus.  The  uniaxial  relaxation  modulus 
was  expressed  in  the  exponential  series  form  described  in  section  5.3.  The 
time-temperature  shift  function,  ,  was  expressed  and  determined  in  the  manner 
described  in  the  same  section. 


Where  stress  predictions  involved  transient  temperature  test  histories,  the 
thermomechanical  coupling  factor,  Af,  was  used  in  equation  (292)  as 


For  UTP-19,360  propellant,  A^  was  found  to  be  a  constant  equal  to  2.2b. 

5.5.1  3D  Analysis  Evaluation 

The  only  3D  application  of  the  -ciified  Swanson  theory  on  this  program  was 
an  instrumented  test  vehicle  subjected  co  a  complex  thermal-mechanical  load  his¬ 
tory.  This  test  was  selected  because  it  represented  the  most  difficult  condi¬ 
tion  for  predictive  stress  analysis.  Other  simpler  3D  tests  such  as  triaxial 
poker  chips  or  tension-torsion  of  thin-walled  tubes  may  have  been  easier  to  analy¬ 
tically  model  but  would  not  provide  a  realistic  check  for  rocket  motor  applica¬ 
tions.  The  selected  3D  applications  highlighted  some  specific  applications  prob¬ 
lems  generated  by  the  stress/strain  axiality  and  amplitude  variations  between 
the  bore  and  bond  area  and  the  corresponding  variation  in  strain  and  strain  rate 
history  throughout  the  grain.  These  could  be  automatically  accounted  for  in  a 
nonlinear  finite  element  code  if  the  theory  were  incorporated.  However,  this 
application  had  to  approximate  these  variations  using  an  effective  stress-strain 
axiality,  stress-strain  amplitude,  and  equivalent  loading  history  to  account  for 
the  material  nonlinearities.  Even  with  these  approximations  the  modified 
Swanson  theory  followed  all  of  the  measured  stress  changes  much  better  than 
linear  viscoelasticity.  The  predicted  peak  stresses  were  above  the  measured 
peaks,  but  the  predicted  values  were  in  agreement  with  motor  stresses  during  the 
hold  period. 
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The  measured  stresses  also  reflected  the  actual  thermal  gradients  in  the 
grain  which  had  to  be  neglected  in  the  predictive  analysis.  If  the  grains  could 
have  physically  changed  temperature  more  rapidly,  then  the  measured  peak 
stresses  would  be  higher  and  in  better  agreement  with  the  predictions.  When  the 
analytic  approximations  required  for  the  analysis  are  considered,  the  predictive 
stresses  are  considered  in  reasonable  agreement  with  the  stress  measurements. 
These  analytical  limitations  would  be  eliminated  if  the  modified  Swanson  theory 
were  incorporated  into  a  nonlinear  viscoelastic  finite  element  code. 
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Appendix  A 


MULTISTATION  AUTOMATED  DATA  REDUCTION 


INTRODUCTION 

Automated  handling  of  multistation  tester  data  is  accomplished  with  a 
system  of  interactive  programs  on  the  HP  9825  desk  top  computer  (Figure  A-176). 
These  programs  include  data  acquisition,  stress  relaxation-master  modulus, 
straining  while  cooling  or  heating,  straining  to  faixure,  and  complex  histories. 
The  acquisition  of  data  and  test  control  are  functions  of  the  data  acquisition 
program  which  supply  data  to  the  data  reduction  programs.  The  reduction 
programs  reduce  and  output  data  for  a  particular  type  of  test  history.  In 
addition,  terminal  emulating  software  for  the  HP  9825  provide  a  data  link  for 
the  transfer  of  data  to  the  VAX-11  mainframe  computer.  This  makes  the  data 
directly  available  to  the  nonlinear  constitutive  theory  programs. 


SYSTEM  INSTRUMENT  CONFIGURATION 

The  multistation  data  acquisition  instruments  are  configured  to  provide 
load,  crosshead  position,  temperature  and  elapsed  time  data  to  the  data 
acquisition  program.  The  system  consists  of  a  Hewlett  Packard  9825  desk  top 
computer,  3**55  digital  voltmeter,  3^95  scanner,  98035  programnable  clock,  and 
9885  flexable  disk  drive  (Figure  A-177). 

The  HP  9825  and  data  acquisition  program  act  as  the  system  controller.  The 
controller  processes  incoming  test  data  and  crosshead  information  and  responds 
by  sending  instructions  to  other  instruments  in  the  system  over  an  HP-IB  inter¬ 
face.  Output  signals  from  the  tester’s  load  transducers,  linear  potentiometer, 
and  analog  thermometers  are  input  into  the  scanner's  programmable  relay  cards. 
The  scanner's  relays  under  command  of  the  controller  can  be  opened  independently 
to  route  output  data  signals  individually  to  the  digital  voltmeter  where  they 
are  digitized  and  read  by  the  program.  Crosshead  control  information  from 
output  lines  connected  to  the  tester's  motor-clutch  assembly,  is  supplied  to  the 
program  through  the  scanner  in  the  same  way  as  the  data  output  signals.  These 
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signals  enable  the  controller  to  react  to  changes  in  crosshead  movement  and 
direction  without  relying  on  operator  intervention.  The  programmable  clock 
connects  directly  to  one  of  the  computer’s  I/O  ports.  It  provides  the  program 
with  elapsed  time  data  and  a  program  interrupt  capability  for  controlling  the 
rate  at  which  data  is  taken.  The  flexible  disk  drive  provides  a  mass  storage 
medium  where  data  is  stored  during  testing  for  latter  access  by  one  of  the  data 
reduction  programs. 


PROGRAMS 
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Data  Acquisition: 

As  previously  stated,  the  data 
acquisition  program  is  used  to  collect 
data  from  the  multistation  tester. 
Operation  of  the  program  involves  steps 
to  initialize  the  program,  calibrate 
the  system,  and  collect  and  store  test 
data. 

Initialization  of  the  program  is 
accomplished  by  operator  entered 
information  used  to  identify  the 
particular  test  and  define  samples 
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being  tested.  In  response  to  promps  from  the  computer  display,  the  operator 
inputs  descriptions  on  test  material,  crosshead  rates,  strain  levels,  and 
temperature  levels  of  the  test  history.  Data  input  on  the  test  samples  include 
their  number,  gage  length,  and  individual  cross  sectional  areas,  along  with 
their  channel  locations.  In  addition,  the  operator  enters  pairs  of  crosshead 
rates  and  delta  strains  for  each  test  interval  used  to  compute  sampling  rates. 
The  operator  also  determines  how  data  is  taken  during  relaxation  cycles  by 
specifying  whether  sampling  is  to  be  done  in  a  fixed  or  log  time  interval. 


Calibration  of  the  system  is  done  by  an  operator-interactive  procedure  to 
determine  the  tester's  transducers  and  potentiometer  sensitivities.  This 
involves  the  operator  queuing  the  program  to  take  readings  from  the  transducers 
at  differing  load  conditions.  By  comparing  the  change  in  output  signals  for  a 

known  change  in  load,  the  lb/volt  sensitivity  of  each  transducer  may  be 

determined.  Similarly,  by  moving  the  linear  potentiometer  probe  a  known 

distance  its  in. /volt  sensitivity  is  determined.  The  analog  thermometers  are 

not  calibrated  at  the  time  of  testing.  These  units  output  a  10  mv/F  repre¬ 
sentation  of  the  test  chamber  and  internal  sample  temperature.  Calibration  on 
them  is  done  periodically  by  the  CSD  electronics  laboratory.  For  short  time  and 
isothermal  tests,  calibration  is  done  once  before  testing  begins.  For  tests 
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lasting  over  a  long  time  period,  a  second  calibration  is  done  when  testing  is 
complete  to  enable  compensation  for  drift  in  the  tester’s  electronics.  Since 
the  load  transducers  are  temperature  sensitive,  for  thermal  tests  two  cali¬ 
brations  must  be  performed  at  differing  temperatures  to  determine  the  change 
in  sensitivity  per  degree  change  in  temperature. 

When  calibration  is  complete,  the  program  stops  operation  until  testing  is 
ready  to  begin.  On  a  queue  from  the  operator,  zero  load  and  position  data  is 
taken  and  the  system's  instruments  programmed  to  their  initial  conditions.  The 
clock  interupt  period  is  set  for  a  sampling  rate  determined  from  the  initial 
crosshead  rate  and  delta  strain  information.  Scanner  relays  are  also  arranged 
to  monitor  the  tester's  break  input  voltage. 

The  program  monitors  the  breed:  voltage  until  detecting  the  break  has 
disengaged  which  signifies  erosshead  motion.  The  clock's  counter  and  interupt 
units  are  then  started.  Interupt  signals  are  output  by  the  clock  at  the  set 
sampling  rate  until  changed  by  the  program  at  the  end  of  the  straining  interval. 

When  interupt  instructions  are  received  from  the  clock,  program  operation 
branches  to  a  data  collection  subroutine.  The  voltmeter  and  scanner  are  set  to 
read  output  signals  from  each  of  the  transducer,  potentiometer,  and  thermometer 
channels.  Fifty  milliseconds  are  required  to  read  each  channel.  Elapsed  time, 
read  from  the  clock  counter,  is  taken  as  the  mean  time  over  which  the  data  set 
was  read.  The  test  data  is  retained  in  a  memory  buffer  until  transferred  to  a 
disk  storage  file. 

The  program  continues  to  monitor  the  crosshead  break  and  clock  information 
channels  throughout  the  test.  When  a  change  in  crosshead  motion  is  detected, 
the  clock  interrupt  is  stopped.  From  the  test  description  data  correspon Jing  to 
the  test  interval,  a  sampling  rate  is  determined  and  the  clock  reset.  For  log 
time  interval  samplings,  the  clock  interrupt  unit  is  stopped  after  each  reading 
and  the  rate  doubled. 
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Up  to  600  data  sets  may  be  retained  in  the  computer's  memory  at  one  time. 
Data  is  transferred  to  the  disk  either  between  data  set  samplings,  if  time 
allows,  or  when  testing  is  complete. 

Data  Reduction  Programs: 

The  reduction  programs  reduce  and  output  multistation  data  pertinent  to 
particular  types  of  test  histories.  The  programs  are  stress  relaxation-master 
modulus,  straining  while  cooling  or  heating,  straining  to  failure  and  complex 
histories.  A  description  of  the  strain  and  temperature  histories  relevant  to 
each  is  listed  in  Table  A - 1 5 . 


Test  identification,  calibration,  load,  sample  extension,  thermal,  and  time 
data  are  supplied  to  the  programs  from  the  acquisition  data  files  or  entered 
directly  by  the  operator.  In  addition,  relaxation  cathetometer  strain 
measurements  and  the  thermal  expansion  coefficient  may  be  optionally  entered. 

Each  program  reduces  stress,  strain,  modulus,  temperature,  and  elapsed  time 
data  when  applicable.  The  method  by  which  each  is  determined  depends  on  the 
test  history  and  amount  of  information  available  to  the  program. 

Calibration  sensitivity  (S)  of  the  transducer  and  potentiometer  are 
determined  in  general  by 

S  =  load/(load  output-zero  load  output) 

where  load  is  the  transducer  calibration  weight  or  potentiometer  probe 
displacement.  For  tests  where  multiple  calibrations  were  performed,  the 
sensitivities  at  time  t  are  corrected  for  electrical  drift  and  thermal 
variations  with  the  linear  relationships 


Sit) 


s 


initial 


(Sfinal  "  sinitial) 
- X  (t) 

xfinal  “  xinitial 
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TABLE  A- 15 


Program  Test  Strain  History  Output 


1 .  Stress  relaxation 


2.  Straining  while 
cooling  or  heating 


Tabular  -  time,  modulus 
Graphic  -  modulus  vs  time 


Tabular  -  time,  strain, 
temperature, 
stress 

Graphic  -  stress  vs  time 
and  temperature 


3. 


Straining  to  failure 


Test  History  - 
Tabular  -  time,  strain  and 
stress 

Graphic  -  stress  vs  time  and 
straining 

Mech  Properties 
Tabular  -  initial  modulus, 
maximum  stress 
and  strain  - 
corrected  stress 
and  strain  - 
rupture  strain 


4.  Complex  histories 


Combination  of  strain¬ 
ing,  relaxation  and 
temperature  intervals 


Tabular  -  time,  strain 

temperature  and 
stress 

Graphic  -  stress  vs  time 
and  strain  and 
temperature  vs 
time 


where  X  is  time  for  isothermal  and  temperature  for  nonisothermal  tests. 
Temperature  corrections  are  not  made  on  the  potentiometer  sensitivity  sinoe  it 
is  located  outside  the  environmental  test  chamber. 


Zero  load  outputs  (ZO)  for  the  transducers  are  also  corrected  for 
electriaal  drift  and  thermal  variation  by  the  same  method  as  the  sensitivities. 
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Corrections  on  the  zero  position  output  of  the  potentiometer  cannot  be  made 
since  the  crosshead  can't  be  accurately  returned  to  its  initial  position. 

Sample  stresses  (a)  at  time  t  are  calculated  by 

o(t)  a  [transducer  output  -  ZO  (t)]  S(t) /cross-sectional  area 
Sample  strain  (e)  at  time  t  is  determined  by 

e(t)  =  [pot  output  -  ZO]  S(t)/gage  length 

For  nonisothermal  tests  a  strain  correction  may  be  applied  using  the  thermal 
expansion  coefficient  (a).  In  this  case,  the  total  sample  strain  becomes 

e(t)  =  mechanical  strain  +  [T(t)  -  a/gage  length 

An  additional  correction  for  effective  gage  length  may  be  made  using  cathe- 
tometer  measurements  of  actual  sample  strains.  The  correction  factor  is 
determined  as  the  ratio  of  the  mean  intervals  in  the  test  history.  For 
histories  where  multiple  cathetometer  measurements  were  made,  the  correction 
factors  are  linearized  to  measured  strain  between  them. 

Relaxation  and  secant  modulus  (E)  at  time  t  is  determined  by 

E(t)  =  a(t)  1  +  e(t)  /e(t) 

where  e(t)  is  held  constant  over  relaxation  test  intervals. 

Temperature  is  reduced  from  analog  thermometer  readings  by  converting  the 
millivolt  output  to  volts. 

Elapsed  time  is  calculated  as  the  difference  between  when  the  data  was 
taken  and  when  initial  loading  occurred  (tQ)  since  loading  times  may  vary  from 
sample  to  sample,  t0  is  approximated  by  the  time  of  initial  straining. 

Once  data  is  reduced,  a  tabular  and  graphic  summary  of  the  test  is  output 
by  each  program.  A  description  of  the  outputs  is  listed  in  Table  A-15«  To 
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retain  data  for  future  reference  and  reuse,  identification,  calibration,  and  raw  test 
data  are  stored  on  permanent  diskette  data  files. 

Terminal  Emulator 

The  program  is  used  to  transfer  data  between  the  HP  9825  and  VAX-11,  desk 
top  and  mainframe  computers. 

A  link  is  created  between  the  computer  types  utilizing  the  VAX-1 I's  dial- 
in  lines  and  an  RS-232  interface  which  connects  the  9825  to  an  acoustic  coupler. 

The  emulator  software  then  supplies  the  capability  of  using  the  9825  as  an 
intelligent  terminal  through  which  data  may  be  read  from  the  flexible  disks  and 
sent  over  phone  lines  to  the  VAX. 

Data  transfer  is  accomplished  with  a  VAX  program  which  reads  data  sent  from 
the  terminal  and  retains  it  in  storage  files  for  access  by  the  nonlinear  theory 
programs . 
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SYMBOLS 


F 

Af 

AiJ 


a  a* 
■q’ 


a 


aF 

ak 

anp 

Ali  Agi  A^,  Aj|,  Ag, 

2a 

B 

B 

B 

B* 

Bi 

bnt 

Bo 

b 

2b 

2b 


microcrack  growth  rate  shift  factor 

temperature-dependent  material  function 

thermomechanical  coupling  factor 

expansion  coefficients  of  bulk  stress  in  terns  of 
octahedral  strains 

constant 

initial  area 

temperature  shift  factor  (a^) 

(AHETA)  damage  related  shift  function 
half  sample  width,  grain  ID  radius 
softening  function 
constant 

expansion  coefficients  of  correction  modulus 
constants 

biaxial  sample  width,  grain  ID 
Cauchy-green  deformation  tensor 
bulk  modulus  or  a  constant 
volume  average  deviatoric  strain 
deviatoric  deformation  tensor 
constant 

Farris  constants 

constant 

constant 

biaxial  sample  height  (gage  length) 
grain  OD 


c 


CSD 

Cx 

C,  Cp  C21  c± 

c 

D 

Dm»  d6>  ^7  Di 
d 

E 


Er 

E(t),  Erel(t) 

e 

d 

®1J*  eiJ 
E($) 

F 


F 


softening  function 

Chemical  Systems  Division 

rehealing  parameter 

constants 

constant 

damage 

constants 

constant 

modulus 

relaxation  function  in  bulk  or  constant 

relaxation  function  in  shear  or  constant 

activation  energy 

case  modulus 

applied  strain 

propellant  modulus 

equilibrium  modulus 

effective  modulus 

reference  modulus  and  normalized  coefficient  for 
modulus 

relaxation  modulus 

linear  viscoelastic  relaxation  modulus 
product  of  F  and  virgin  response  function  g 
deviatoric  strain  tensor 
linear  viscoelastic  modulus 
damage  function  or  softening  function 
degrees  Fahrenheit 
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force 


* 


F 

F(t) 

M,  *„) 

FE 


fv  f2 

f 


(f) 

*c 

f(t) 

G 


G.  L. 


G 


c 


rel 

G(t) 


g 

g(  ) 
(g) 

®L 

gu 

H 

h 


h 


strain  magnification  factor 

constant  rate  modulus 

damage  curve  at  damage  level 

finite  element 

function 

material  parameter 
deformation  function 

constant  time  rate  of  change  of  deformation  invariant 

viscoelastic  type  function  in  kinetic 

shear  modulus 

gage  length 

corrected  modulus 

equilibrium  modulus 

constant 

relaxation  modulus 
shear  relaxation  modulus 
shear  relaxation  modulus 
virgin  response  function 
function  of 

strain  softening  function 
softening  function  for  loading 
softening  function  for  unloading 
function  of  damage 
case  thickness 

function  of  damage  in  kinetic  equation  of  evolution 
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HTPB 


hydroxyl-terminated  polybutadiene 


I 

**  *7  **  pi 
J 

J1 

J2 

JANNAF 
K,  KrK10 
K 

*1 

Kij  (xk) 
Kij  (*k^ 

K* 

k 

L 

Lx 

LCTE 

M 

Mx.  M2*  Ml, 
m 

m,  m^  i 

N 

N-L 


volume  dilatation 

contribution  to  stress  at  time  tn 

octahedral  shear  strain 

Lp  norm 

creep  function 

creep  function  in  bulk 

creep  function  in  shear 

Joint  Army  Navy  NASA  Air  Force 

constants 

temperature-dependent  constant 

stress  intensity  factor 

stress  at  point  xk 

strain  at  point  x^ 

rehealing  parameter 

constant 

length 

constant 

linear  coefficient  of  thermal  expansion 

constant 

constants 

constant 

material  parameters 
number  of  cycles 
nonlinear 


NLCT 


nonlinear  constitutive  theory 


n 


constant 


P 

P 

P2  *  P1J 
P15 

pr 

PBAN 

P 

Q 


S 

S 

s 


m 


SfSr 


SF 


terms  of  equation  under  summation 

hydrostatic  pressure 

constants 

used  to  normalize  to  1 
finite  length  correction  factor 
polybutadiene  acrylonitrile 
constant 

terms  of  equation  under  summation 

constant 

function 

healing 

lower  limit  of  R  (healing) 
strain  maximum  sensitivity 
root  mean  square 

chemical  bond  radius  or  rotation 

virgin  stress  and  damage  parameter  constant 

damage  parameter 

damage  parameter 

maximum  damage 

resultant  stress 

constant 

certain  measure  of  damage 
temperature-shifted  time 
constant 
stress  free 
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temperature 


T 

Ta 

Tc 

T6 

Tn»  Th 

To 

To 

tR 

TCM 

TSM 

t 


t* 


to 

U 

UTP 

V 

^r*  xr 
x 

x»  y»  z 

Yi»  Y1>  y2 »  y3 

y 

a 


peak  stress  time 

material  property  and  shift  temperature 

motor  cure  temperature 

glass  transition  temperature 

healing  temperature 

temperature  at  t  =  0 

stress-free  temperature 

reference  temperature 

thermo rheo logically  complex  material 

thermorheo logically  simple  material 

time 

time 

healing  time 
modified  time 

time  to  failure  under  constant  load 
chemical  bond  stretching  rupture 
United  Technologies  Propellant 
volume 
root  of  Y3 

sC/em 

coordinate  axes 
functions  related  to  damage 
=*m,  maximum  strain  to  time  r 
kinetic  constant 


coefficient  of  thermal  expansion 


material  parameter 


a,  «i 

**e 

“P 

P 

P 

P 

7 

7 

7 

T'a 

'a 

AE 

AL 

AT 

AT0 


e1 

*2 

*11*  *22*  *33 


constant 

case  coefficient  of  linear  expansion 
propellant  coefficient  of  linear  expansion 
constant 
parameter 

changing  rate  coefficient 
constant 
shear  strain 
softening  function 
adhesive  failure 
cohesive  failure 
modulus  change 
change  in  length 
temperature  change 
thermal  stress 

negative  volume  change  with  hydrostatic  pressure 

constant 

case  qc  AT 

Kronecker  delta 

propellant  ap  AT 

strain 

strain  rate 

principal  strain 

lateral  strain 

principal  strains 
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£°,  €<>r 

<m*  emax 
«(t) 


ee 

f 

v 

V(.t) 

e 

e 


pseudo  strain 

maximum  strain 

strain  at  time  t 

strain  of  unfilled  polymer 

strain  due  to  mechanical  stress 

bore  strain 

neutral  rate  where  damage  remains  constant 
compressibility  (  17  =  0  is  incompressibility) 
related  to  damage  function 
ratio  of  reload  virgin  stress 
axis 


summation 


a 


engineering  stress  or  stress 


o 


°B 


°b 


oc 

a.  d 
ij 

o* 


°o 

°kk 


Or 

Or  r^ 
o(t) 

of(t) 


T 


r 

0 


0 

01  .  02 
0 


^n 

OJ 

Hq 

VnB* 


Cauchy-stress  tensor 
bulk  stress 

stress  correction  due  to  change  in  state  of  bonding 

stress  correction 

deviatoric  stress 

deviatoric  stress 

linear  viscoelastic  stress 

constant  stress 

bulk  stress 

linear  viscoelastic  stress 
case-grain  radial  bond  stress 
stress  at  time  t 
fading  memory  stress 
stress 

reduced  time 
shear  stress 
shifted  time 
function  of  loading 
extension  ratio 
functions 

rate  of  loading  function 

nth  component  of  stress  correction 

normalized  damage  function 

second  invariant  of  tensor  ( a  =  a ,  B) 

second  invariant  of  deformation  tensor 
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*  • 


second  invariant  of  deviatoric  stress 


denotes  absolute  value 


infinity 


